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ABSTRACT 


A new capability has been added to the general purpose finite element 
program NASTRAN Level 17.7 to conduct forced vibration analysis of tuned 
cyclic structures rotating about their axis of syrametry. The effects of 
Coriolis and centripetal accelerations together with those due to linear 
acceleration of the axis of rotation have been included. 

This report presents the Theoretical, User's, Programmer's and 
Demonstration manuals far this nsv/ capability. The work was conducted 
under Contract NAS3-22533 from NASA Lewis Research Center, Cleveland, Ohio 
with Mr. Richard E. Morris as the Technical Monitor. 
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FORCED VIBRATION ANALYSIS OF ROTATING CYCLIC STRUCTURES 


1 . 1 Introduction 

A new capability has been developed and implemented in NASTRAN Level 17.7 to 
perform forced vibration analysis of cyclic structures rotating about their axis of 
symmetry. Fans, propellers, and bladed shrouded discs of turbomachines are some 
examples of such structures. The capability ihcludes the effects of Coriolis and 
centripetal accelerations on the rotating structure which can be loaded with: 

1) directly applied loads moving with the structure and 

2) inertial loads due to the translational ecceleraticti of the axis of 
rotation ('base' acceleration). 

One rotational ly cyclic sector of the N-sectorad structure is modelled and 
analyzed. Steady-state sinusoidal or general periodic loads are specified to represent 

1) the physical loads on various segments of the complete structure, or 

2) the circumferential harmonic components of the loads in (1). 

The sinusoidal loads are specified as functions of frequency and the general 
periodic loads are specified as functions of time. 

The translational acceleration of the axis of rotation may be specified as a 
function of frequency in an inertial coordinate system. 

The details of the User's, Programmer's and Demonstration manuals are presented 
in Sections 2, 3 and 4, respectively. The fc'*lowing sections present the salient 
points in the theoretical development. 

1 . 2 Theory 

The theoretical development of Reference 1 to conduct forced vibration analysis 
of rotating cycl ic structures in conjunction with tfie theory of rotational cyclic 
symmetry as presented in Section 4.5.1 of the NASTRAN Theoretical Manual (Reference 2) 
is summarized in this section. 

1.2.1 Equations of Motion 

The complete structure consists of N identical sectors. The displacement at 
any grid point in any sector can be expressed in any body-fixed coordinate sysiem 
as a combination of: 

1) the steady displacement due to the steady rotation of Lite structure, and 

2) the vibratory displacement (superposed on the steady displacement) due to 
the vibratory excitation pro/ided by the directly applied loads and base 
acceleration. 
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Tl’.s vibratory resoonsri of rotating cyclic structures may be determined by 
this new capability. 


As shown in Reference 1, the equations of forced response can be written as 

. n = 1. 2 ri. (1) 


+ b"u" + 


For the n^^ cyclic sector, represents the vibratory degrees of freedom; 
b" and k” represent its mass, damping and stiffness matrices respectively; 

P” represents the directly applied loads on u"^, and -M^R represents the inertial 
loads on u'^ duo to base acceleration K. The damping matrix consists of the 
viscous and structural damping, and the contribution due to the Coriolis 
acceleration, i.e.. 


b" = b'’'. 

viscous , 


+ 2rsl! . , . , 

Conol 15 ’ 


( 2 ) 


structural 

with SI as the (constant) rotational speed. The stiffness matrix ccnsists of 
elastic and differential stiffness together with the contribution duo to the 
centripetal acceleration, i.e., 

(3) 


j/ n ^ 1 / n . i/H 

'^elastic ^^differential 


'■ centripetal 


The derivation of the coefficient matrices b!^!. . and 

conoiis centripetal 2 

is given in Refererse 1. 

Equations (1) supplemented by the inter-segment boundary compatibility 
conditions (Section 4.5.1, Reference 2), 

“side 1 " “side 2’ " " 

completely describe the vibratory forced motion of the rotating cyclic 
structure. 


1.2.2 Method of Solutio n 

The method of solution of equations (1) consists of four principal steps: 

1) Transformation of applied loads to frequency-dependent circumferential 
harmonic components. 

2) Application of circumferential harmonic-dependent inter-segment 
compatibility constraints. 

3) Solution of frequency-dependent circumferential hani'onic components of 
displacements. 
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4) Recovery of frequency-dependent response (displacements, stresses, loads, 
etc.) in various segments of the total structure. 

An overall floivchart outlining the solution algorithm is shown in Figure 1. 

1. Transformation of Aoolied Loads 


a 


u 

ty 

11 


(3 

M 


I 


I 
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The transformation to frequency-dependent circumferential harmonic com- 
ponents depends on the form in which the excitation is specified by the user. 
The following options are made available in the present capability to specify 
the form of excitation due to the directly applied loads and base acceleration 
loads: 

Directly applied loads specified as; 

- periodic functions of time on various segments 

- periodic functions of time for various circumferential harmonic indices 

- functions of frequency on various segments 

- functions of frequency for various circumferential harmonic indices. 

Base acceleration specified os: 

- function of frequency for circumferential harmonic indices 0 (axial) 
and 1 (lateral ). 

Details of each of the above five loading conditions are as follows. 

Directly applied loads (segment-depen d ent and periodic in time) 

If represents a general periodic load on sector n specified as a func- 
tion of time at M equally spaced instances of time per period (Figure 2), 
the load at time instant can be written as 


m 


-0 


■’l 

p'T = + Y 

£=1 


-£c 


■£s 


P cos(m-l J,b) + P sin(m-l;.b) 


-M/2 


+ (-D'^'V” 


(5) 


m = 1 , 2, 


. M 


I! 


a 


where b = 27;/M, Z, - (M-l)/2 for odd M, f,, = (M-2)/2 for even M. The last 

*• •• ^ 1 1 3^ I I 

term in equation (5) exists only when M is even. The coefficients pn'' 

("jl" = 0; Z.C, Zs, t- 1 , 2, ..., M/2) in equation (3) are inoependent of 

time, and are defined by the relations 
-0 m 


P '1 = 1 y pi 
m~l 


(z - 0) 


part of (G) 
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P “ u z. P -t- ^ (i-i ' i '-1> ' :i 

‘ m=l 


= i I ' ‘T ( -•■• • -t' ) , and 

" m=1 


{(.= 1 , 2 , .... ilj^) 


(6 Contd.) 


p" "" M ^ (“’!)" ‘ (H even only) (£=M/2). 

ni“l 

Each of the coefficient vectors P*' on the left hand sides of equations 
(6) can further be expanded in a circumferential (truncated) Fourier series 

-*‘2,** ' "'*5-** 

pii _ pO ^ j pi< (rPfka) + sin(n-lka) + (-1 )'^'^ ^ 

k=l L J 


where n = 1, 2, .... N , ' 

= 0; ilc, «,s.<l = 1, 2, .... P.^; M/2 
a = 2n/N 

k^ = (H-l)/2 for N odd [ 

l<j=(M-2)/2 for N even. 

«< 

The last term in equation (7) exists only when N is even. The Fourier 


coefficients P ^ ("k" = 0; kc, ks, k = 1, 2, .... k^_; N/2) in equation (7) 

do not vary from sector to sector, and are defined by 


= 1 y p*^ 
N 1 , 


(k = 0) 


okc 2 " r,n ; — , 

P ■ TI > P cos(n-lka) 


~ sin (n-lka), and 


(k = 1, 2, 


... k^) 


" il } f" (N even only) (k = N/2) . 


\ 



v.T.'G'm :A1- ‘1 

Or POOn 


The terms P (‘'£'' = 0; Px, Ss, 7, = 1, 2, M/2 and "k" = 0; kc, ks, 

k - 1 , 2, ...j k^^; N/2) are tha transf ormed frequency-dependent circumferential 
harmonic components of the directly uppliGci loads pi"! (m = 1, 2, .... M and n = 1, 
2, N). 


li recti V applied loads (Circumferontial han’xnic-dependent and periodic in 


time). 


Such loads can be represented as 
n; -0 l. 

(ll,!! ULIi » H|^t: 

nl' _n^ ^ ^ ^ ^ 


“ilS 

^ J - 


-M/2 

>m-lK"k" 


P '' = P ^ ^ T L + (-1) P . ( 

where m = 1, .2, M represent the time instances at which harmonic components 

"k" = 0; kc, ks, k = 1 , 2, . . . , k, ; N/2 of directly applied loads are specified. 

-"i" 

The coefficients p-"k" on the right hand side of equation (10) are obtained 
using equations (6) with sector number n replaced by harmonic number "k". 


Directly appl i< 


/-and segment-dependent) 


This type of loads can be represented as 


p'^ = P^ + P^*^ cosCiv^ka) P^^ sin(rrTka) + (-1)'^"^ , (11) 

k=l L 

where "t" (=1, 2, F) now represents the frequencies at which excitation is 

specified. The transformed frequency-dependent c’rcumferential harmonic components 
-"X." 

•p k _ Q. k = 1, 2, ..., k^; N/2) are obtained using equations (9) 

with "X." as defined above. 

Directly applied loads (frequency- and circumferential harmonic-dependent) 

These loads are the transformed frequency-dependent circumferential harmonic 
components ("k" - (1; kc, ks, k = 1, 2, ..., k^; N/2) with "X," (=1, 2, ..., F) 

representing the various frequencies at which the directly applied loads are 
specified. 

Base acceleration (freque ncy- a nd circumferential harr.’ionic-dependent) 

In Reference 1, it is shown that the components of the translational base 
acceleration contr'bute sc inertial loads on the rotating structure in the follow- 


ing manner-; 


f? 

L? 
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1. Axial component contribiites to P ' vjhiU'e "k" = 0, and ”Z“ represents the 
specified excitation frequencios._,i^,i 

M|.<l 

2. Lateral components contribute to P where "k" = Ic and Is, and "Z" re- 
presents the effective excitaiion fraquericies which are shifted from the 
specified frequencies by :!: n, the rotational frequency. 


The user specifies the components of the base acceleration vector R as 
functions of frequency. The program computes the inertial loads -H^R and trans- 
forms them to appropriate frequency- dependent circumferential harmonic components. 


2 . Application of Inter-Segment Compatibility Constraints 


As shown in Section 4.5.1 of Reference 2, equations (4) are used to derive 
the compatibility conditions relating the circumferential haniicnic component 
degrees of freedom on the two sides of a rotational ly cyclic sector: 
side 2 side 1 






(k - o) 


2 


rks _ 


u^‘* cos(ka) sin(ka) 

-u^*' sin(ka) + cos(ka) 


Wk = 1. 2. .... k^) V (12) 


and 


^ 1/2 

^2 




(k = N/2) 


In order to apply these constraint relationships for any given harmonic k, 

an independent set u' consisting of the circumferential harmonic component 

(cosine and sine) degrees of freedom from the interior and side 1 of the cyclic 

1/ 

sector is defined, u is selected from the 'analysis' set degrees of freedom, 
and is defined as 


= G^^(k) and 

^s = G3|,(k) ^ . 


(13) 


- — 1 / c 

ir and u each contain all (and only) the 'analysis' set degrees of 
freedom from the interior and both sides of the cyclic sector. Equations (12) 
are used to define some of the elements of the transformation matrices G^,j, and 
Gsk' For k = 0 and N/2, the matrix G, j, is null. 
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3 . Solution or Freguency-De o e_i i done iiErn'Oii'i c ui., placements 

/ 

For a given harmonic k, the introduction of IT' in the equations of motion, 
(1), results ill the trar.sfonnch equations of .’notion (Reference 3) 


where 



+ B 

!<.H< 

u 

+ 1^/' 

P^ 


# = 


h" 

“ck 

+ 

Cvu m" 

^sk 


‘^cl; 

b" 

^ck 

T 

,-T nH 

“sk 


T 

k'' 

^ck 

f 

pT ./D 

^sk 

^s'- 

II 


p'kc 

■ + ( 

,T 

^sk 

pks 



(14) 


('iKi 

\ t 


As discussed in subsection 1 of Section 1,2.2, P' and P are the transfo-^iGd 
frequency-dependent circumferential harmonic components of the directly applied v.nd 
base acceleration loads. 

At any excitation frequency o/, let 


^ and accordingly. 




;i6) 


=K 

where P and u are complex quantities. Equation (14) can be rewritten as 


r .2 ttK 

l-w M 




The excitation frequency w" is given by 


(17) 


tij' = w for all directly applied and axial t j6 acceleration 
loads, and 

= w±il for lateral base acceleration loads. 


(18) 


Equation (17) is solved for u for all excitation frequencies and ?.11 harmonics 
as specified by the user. Tiie cosine and sine harmonic components of displacements 
are recovered using equations (13). 


4 . Recovery of Frequency-Depende nt Dis placements in Various Segments 

This step is carried out only when the applied loads are specified on the 
various seernents of the complete structure. 
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For loadsjpecifliid as functions of time, equation (7) is used to obtain the 
displacements in various segments with = 0; £c, is, £ = 1 , 2, . . . . z 

nidx 


varied from k . to k 

inin max 


I he user specifies 


The circumferential harmonic !: i; 

^max’ '^min *^max* 

^ For loads specified as^functions of frequency, equation (11) is used to 
ODtain the displacements ^^ 11 “ in various segments wi*!; 'wv representing the 
excitation frequencies. The circumferential harmonic is varied from user speci- 

The recovery of other responses such as stresses, internal forces, etc., is 
identical to that currently existent in NASTRAN. 
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FORCED, VIBRATIOfI ANALYSIS OF ROTATING CYCLIC STRUCTURES 

2.1 Introduction 

Forced vibration analysis of cyclic structures rotating about their axis of 
symmetry can be conducted using this capability. An extensive package of ALTERS 
including new Functional Modules and PARAMeters has been developed and used to 
modify the existing Displacement Approach Rigid Format 8 (Direct Frequency and 
Random Response), Series R, in NASTRAN Level 17.7. Example 2 of the Demonstration 
Manual (Section 4) illustrates the use of this ALTER package. The Theoretical 
and Programmer's Manuals are described in Sections 1 and 3, respectively. 

2.2 NASTRAN Model 

The user models one rotational ly cyclic sector (segment) of the entire 
structure as shown by the 12-bladed disc example in Figure 1. ATI NASTRAN 
coordinate systems (basic, location and displacement) are considered fixed to 
the rotating structure. The only additional requirement is that the X-axis of 
the basic coordinate system be coincident with the axis of rotation. A positive 
value of the rotational speed (PARAM RPS) indicates a clockwise sense of rotation 
when the structure is viewed at in the positive basic X direction. 

Except for the special features discussed in this section, the general 
rules of modelling rotationally cyclic structures in NASTRAN (e.g. CYJOIN, 

NSEGS, etc.) have been maintained. 

The rotating structure can be loaded with steady-state sinusoidal or general 
periodic loads classified as : 

1, directly applied loads moving with the structure, and 

2. inertial loads due to the translational acceleration of the axis of 
rotation ('base' acceleration). 

The sinusoidal loads are specified as functions of frequency using the 
RLOADi.bulk data cards. The general periodic loads are specified as functions 
of time using the TLOADi bulk data cards. 

The following notes apply when using TLOADi bulk data cards: 

1. Time delay t must be set to zero. 

2. In conjunction with the TSTEP bulk data card, TLOADi information is 
used to discretely define P(t) at M time instances as or p'^‘ 

(m - 1, 2, .... M), as discussed in Section 1.2.2 of the Theoretical 



Manual. 

N(1) of TS7EP bulk data card = M-2 
DT(1) of TSTEP bulk data card = (72 - 71 )/M 

. 3. P(t) is defined in the interval [71, 72] with (72 - 71 ) as the period. 

4. Only one physical 7S7EP bulk data card is allowed, i.e. continuation of 

the 7S7EP card is not permitted. 

7he following options are provided to specify the form of excitation; 

Directly applied loads specified as: 

- periodic functions of time on various segments (PARAH CYCIO = +1) 

- periodic functions of time for various circumferential harmonic indices 
(PARAM CYCIO = -1) 

- functions of frequency on various segments (PARAM CYCIO = +1 ) r 

~ functions of frequency for various circumferential harmonic indices 
(PARAM CYCIO - -1) 

Base acceleration specified as: 

- function of frequency for circumferential harmonic indices 0 (axial) 
and 1 (lateral) (PARAM CYCIO = -1 ) 

7he base acceleration refers to the translational acceleration of the axis 
of rotation, and is specified in an inertial coordinate system. 7he user defines 
a rectangular inertial coordinate system v/ith its X-axis parallel to and in the 
direction of the basic X axis, as shown in Figure 1. 7he definition of this 
inertial system, otherwise, is arbitrary. 7he user specifies the X, Y, Z com- 
ponents (magnitude and phase) of the base acceleration vector as functions of 
frequency on 7ABLEDi bulk data cards. 7he use of these tables is activated by 
the PARAMs BX7ID, BXP7ID, BY7I0, BYP7ID, BZ7ID and BZP7ID. 

7he user is provided with two options to include damping by specifying the 
form of the matrices B^^ and in the Functional module GKAD as per 

equations 16 through 21, pages 9.3-7 and 9.3-3, Section 9.3.3 of the NAS7RAN 
Level 17.7 7heoretical Manual. 7he PARAMeters GKAD and LGKAD have been defined 
for this purpose. 

Section 2.4.4 of this manual describes all the PARAMeters applicable with 
this new capability. 
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2.3 Subcase Definitions 

The PARAHeters CYCIO (=-1) and Kf-lAX (>0,$NSEGS/2 for even NSEGS, ^(MSEGS-1 )/2 
for odd NSEGS) determine the number, order and meaning of subcases as follows: 

CYCIQ=+1 

The number of subcases is equal to NSEGS, independent of KMAX. 

SUBCASE 1 (SEGJOT. MO. 1) 

SUBCASE 2 (SEGMENT NO. 2) 


SUBCASE NSEGS (SEGMENT NO. NSEGS) 

CYCI0=-1 

The number of subcases is equal to FKfiAX, where 
FKHAX = 1, if Kf'lAX = 0, 

« 1 + 2 . KH.AX, if 0<KMAX<(NSEGS-l)/2, NSEGS odd, 

= 1 + 2 • KMAX, if 0<KMAX^(NSEG3-2)/2 NSEGS even, and 
= NSEGS, -if YJm = NSEGS/2, NSEGS even. 

SUBCASE 1 Ck' = 0) 

SUBCASE 2 ('k' == Ic) 

SUBCASE 3 ('k' = Is) 

SUBCASE 4 Ck’ = 2c) 

SUBCASE 5 ('k' = 2s) 


SUBCASE FKMAX ('k' = KMAXs) 

In the event that NSEGS is even and KMAX = NSEGS/2, Subcase FKIWX will 
represent 'k' = KHAXc as KMAXs does not exist. 

Directly applied loads on various segments (CYCI0=+1) or their circumferential 
harmonic components (CYCI0=-1) are specified under the appropriate subcases. 

With RLOADi bulk data cards, null loads need not be specified by the user. With 
TLOADi bulk data cards, the user is required to provide information to generate 
null loads where applicable. 
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Base acceleration 
PAR/>J-5eters BXTID etc., 
and assigned to = ( 


Is included only when CYCro=-l. 
the corresponding inertial loads 
. Ic and Is as applicable. 


Based on the activating 
are internally calculated 












^ M O' O »*> 0^ V) f 


0RIGH^± ?AC-'i rS 
OF P003 QUALfTf 

2.4 Rigid Format Dgscription 

2.4.1 Rigid Format Alters to Displecemant SOL 8 

1 ^ i 

i LEvii.NfJlNu Gl- KF ALTER £61 - f<F 8 / SEKitS R iLl7,7) / l-28-o2 / M.G. S 
S i 

i i 

FURPCSE - TO MC-UiFY THE OiRtCT FREQUENCY ANL. RANOOM RESPONSE RiGIU S 

EOUMAl TO ENAtJLE THE (JSEk TC PtHFOKH A FCRCEU ViaRATION S 

KcSPGNSE ANALYSIS OF ROTATING CYCLIC STKUClURtS. i 

& 

EXECLTIVL CECK INPCT - 


1 . SOL e 
2* R»f • ALTLkS 


CASE CONTROL OLcK INPUT - 




1. ALL SPC AND MFC CONSIRAINIS MLST GE ABUVE THE SUbCASE LEVEL. 

2. EITHER FkEUUENCY Ok I STEP MUSI 3E SELEcTEi; ANJ MUST BE AbOVE 
IHL SUBCASE LtVEL. 

j. IF SLLE'CTED. FhEuUtNCY MUST BE USED TO SELECT UNE ANu CNLY 
ONE FrEw. FKfwi CK FKEU2 CARD FRcM THE BULK DATA OEOK ANU 
MUST BE UEFU.tO AbOVE THL SUBCASE LEVEL. 

IF SELECrOU, TSJLP MUST BE USEU TO SELECT THE TIMt-STEPS TO BE 
UStU FOR LOaC DLFTMTILN ANU MUST Be. UEFINEU ABOVE THE SUbCASL 


LEVEL. 

S. UIRCCr IhPUT fiAIKiLES ARE NOT ALLO'nLO. 

O. OFREUUENCY MUgT NOT BE USLU. 

7. A SuFARAlc GROUP OF SUBCASES MUST BL uEFINLC FLK EACH SYMMETRIC 
SECMLNT. 

B. ULOAD Must bL USLu TO DEFINE A FrECUlNCY CR T I ME-UE PENOENT 
LOAUING CuNDITIUN FOR EACH SoBCASE. 

FUR rRLCCLNCY-UbPENDLNT LOADS. jUbCASES rITHCJT LOADS NEED NOT 
rlfer to a uloau card. 

FOR TlME-UEPENUtNI LCAuS, SUbCAScS wITHOUl EOAUS MUST REFER TO 
A clOhu card ihat generates a null load. . 

■J. AN ALTEKNATL LUAUUmu METHOD IS TO UEFINE A SEl'ARATE GROUP OF 
SUBCASES FUR LrOH HARMuNIC INDEX, K. IHE PARAMETER CYCIO IS 
INoLOUED AND THL LUAO COMPONENTS FOR EACH INDEX ARE DEFINED 
OIKEOTLY RlThlN EACH GROUP FUR THe VARIOUS LOaUING CONCIIIonS. 


1 
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ORIGir'JAL >h{ 
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C U N T R C L 


■ t'*. 


DECK 


ECHO 



1. SUFUR] oULK DATA CARuS ARE WOT ALLOWED. 

-i. EPUINT BULK DATA CAROS ARE NUT ALLOWED. 

3. SPOINT BULK UAIA CARDS ARE NOT ALLCwED. 

4, CYjOIN bulk DATA CARDS ARE KEgUIREO. 

b. IF A I STEF CARD IS USED THEN IT MUST NUT BE CUNTINJED SINCE 
ONLY ONE CNIFCRH TIME STEP INTERVAL MUST BE SPECIFIED. 

THE SKIP FACTOR fu« CUTFCT, NG . OK THE ISTEP CARD MUST BE 1. 

0. PARAMETERS OSE-U ARE - 

A, NSEGS - KEuUiREO - IME INTEGER VALUE OF THIS PARAMETER 
IS THE NUMBER CF IDENTICAL SEGMENTS IN THE 
SIRUCTUKAE MUOEL, 

fa. CYC 10 - REQUIRED - THE IWIEGER VALUE UF THIS PARAMETER 

SPECIFIES THE FORM OF THE INPUT AND OUTPUT DATA. 

A VALUE OF Cl IS USED TO SPECIFY PHYSICAL SEGMENT 
REPRESENTATION. A VALUE CF -1 IS USED TU SPECIFY 
CYCLIC TRANSFCkMATIUN REPRESENTATION. THERE IS No 
uEFAULif A VALUE MOST oE INPUT. 

L. CYCSEj - FIXED - THE INTcDER VALUE UF THIS PARAMETER 
SPECIFIES THE PROCEDURE FOR SEQUENCING THE 
EQUATIONS IN THE SOlUMCN SET. THE VALUE OF THIS 
PARAHETlr has been set to -1 TC SPECIFY 
ALTERNATINL. COSINE ANU SINE TERMS. 

D. CTYPE - FIXED - THE aCD VALUE OF THIS PAKAMETEk 

DEFiliES THE TYPE CF CYCLIC SYMMETRY. THE VALUE 
Of IfUS PARAMETER HAS BEEN SET TC -kOT- FOR 
ROTATIONAL SYMMETRY. 

E. KMAX - RECjIRED - THE INTEGFR VALUE OF THIS PARAMETER 

SPEuIFILS The MAXIMUM VALUE OF THE HARMONIC 
INDEX. THERE IS NC UEFAULT FOR THIS PARAMETER. 

THE MAXIMUM VALlL fHAT LAN fat SPECIFIED IS NSEGS/2 

F. RMIN - OPTIONAL - THE INTEGER VALUE OF THIS PARAMETER 

SPElIFIES THE MINIMUM VALUE OF THE HARMONIC 
INDEX TC BE USED If. THE SOLUTIC.N LuuP . KMIN CAN 
EQUAL KMAX. THE uEFACLT VALUE IS 0. 

0. LMAX - OPTIONAL - TfiE INTEGER VALUE OF THIS PARAMETlR 

SPECIFIES THE MAXIMUM TIME HARMONIC INDEX. THE 
DEFAULT VALUE IS NTSIEPS/2. wHERE NTSILPS EQUALS 
N (FROM TSTLP CARUl PLUS 2. 

H. NLUAD - FIXED - THE INTtGuR VALUE CF THIS PARAMEI ER 
IS rUL NUMBER OF LOADING CONDITIONS. THL VALUE 
UF THIS PAKAMEIEK is internally CALCULATED. 

1. RPS - UPTILNAL - THE RLAt VALUE CF THIS PARAMETER 

ULFINLS THE RGTATIUNAL SPEED CF THE STRUCTURE 
IN kEVCLoTICNS per uM I TIMl. I hE DEFAULT VALUE 
IS O.C . 

J. dXTiD - OPTIONAL - JHE POSITIVE INTLDER VALUES OF THcSL 
BYIiD PARA.METEkS DEFINE THE SET lULNT iP ICaI I UN NUMBERS 

BZTID UF THE TABLEUl BULK DAIA CARDS iiHICH DEFINE THE 

BXPTiU COMPCNtNTS CF THt uaSL ACCELERATION VECTOR. THE 
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CF POOR QUALfT/ 

NASTKAN tXECUriVE CONTKQL OECK ECHO 


tlYPT Ib 

llZPTIO 


K. HUKPKT 


L. GKUPN7 


M. HfMASS 


N. CUJHMASS 
Li c uK/VD 


?. LbKAU 


R . ri3 


TAocSS KEFEREU 10 oV BXUu, dyUi) ANU BZT 10 
DtFlWt MAGNiruJt(LT-2) AND ''HE TABLES KEFERLDIO 
BY QXPIIOt bYPTiO AND BZPTIC DEFINE PHASE ( DEGREE I , 
IHL BEFALL r VALUES ARE -I WHICH PEANS IHAI 1 HE 
kESPECTIVE terms are IGivOREU. 
nPTIGNAL - AN INTEGER VALUE CF f.l FOR THIS 
PAkaME TCK HILL CAUSt THE CUrtKENT HARMONIC INDEX e 
XINDE/. r TO BE printed AT THt TOP OF THE HARMONIC 
LOUP, the default VALUE IS Cl. 

OPTIONAL - POSITIVE INTEGER VALUE OF THIS 
PAKAPcTuR WILL CaUSE THE GRID POINT rtEIGHT 
GENERATOR TO BE EXECUTED AND THE RESULT I Nu WEIoHT 

balance infukmaticn to be printed, default is -1. 

OPTIONAL - THE TERMS OF THE STRUCTURAL MASS 
matrix AiiE MCLTIFLILO BY THE HEAL VALUE CF THIS 
PARAMETER wHEN THEY ARE GENERATED IN EMG. THE 
DtFAULT IS l.C . 

FIXED - ONLY LUMPED MASS MAIKICES MUST BE USED. 
OPIlUNAL - THE BCD VALUE CF THIS PAKAHElcR IS 
USED TO TELL THE C-KAD MODULE THE DESIRED FORM OF 
MATRICfcS KDD t BUD AND MUD. THE BCD VALUE CAN BE 
FRCDkESP or TRANRESP. the DEFAULT IS TkANRESP. 

ROTL - REMEMBER TO DEFINE PARAMETERS G. W3 AND W4. 
iLE SECTION 9.3.3 (DIRECT DYNAhIC MATRIX 
ASSEMULY) PAGES 9.3-7 AND 9.3-8 OF THE 
NASTRAN IHEORE TIDAL MANUAL. 

OPTIONAL - THE INTEGER VALUE OF THIS PARAMtIEK 
IS USED IN CCNJUNCTILN WITH PARAMETER GKAD. IF 
GKAu=FP^EQRE SP THEN jLT LGriAD=l. IF Gn AD = T KANkESP 
IHEN AET LGKA0=-1. JHE DEFAULT VALUE IS -1. 
optional - THE REAL VAlUE Cs- THIS PARAMETER IS 
USED AS A UNIFORM STRUCTURAL DAMPING COtFF IC I ENT 
IN THE DIRECT FCRMULATICN OF DYNAMICS PROBLEMS. 
OPTiCNAL - ThL real VALUE LF THIi PARAMETER IS 
UStU Ai A PIVOTAL FREwUENCY FCR UNIFORM STRUCTUAL 
DAMPlfJL, IF PARAMEILK uKAD=T:!ANRESP. IN THIS CASE 
W3 iS KECUIRED IF U.VlFUKSED >TRuLTUAL DAMPING IS 
UEsIRbO, the default VALUE IS 0.0 . 

OPTIONAL - THE REAL VALUE CF THIS PARAMETER IS 
USED AS A PIVOTAL FKLUUENCY FOR ELEMENT STRUCTUAL 
damping IF PAkAMEItk GAAD=TRANRt;GP. IN THIS CASE 
WA iS KtCUlRLD IF STRUCTUAL DAMPING IS DESIRED FOp. 
ANY CF THE STRUCTUAL ELLMEM3. DEFAULT IS O.C . 




KLM.ARnS - 

1 . THL A.'i Al YSIj wILL i_ JOP 
KSNDlA = RMIN TO KMaX, 


TiTRU 


RANGE 


ThL CYCLlu TNDlX i 
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$6$$ $$($$»$$$$£>£$&$$ :>S 

ALTER 3 £ 

PILL : UXVF=APPtND/PD (=ARRENO/PO-APPfcND $ 
i PLKHGKfi INITIAL LKKCR CHECKS uN NStCS AND KMAX. 

CONU EkRtiRClfNShuS i IP USER HAS NOT SPECIPIED NSEGS. 

ULNU EKROkCI iKMAX i IF CSEk HAS NOT SPE Cl P I ED KMAX * 

PARAM //C.N.6U /VfN.CYCIOERR /V,Y,LYCIU=0 /L.N.O $ 

CCNU EHKORCl.CYO iOEKR $ IF USER HAS NOT SPECIFIED CVCIO. 

PARaM //CtiMfDlV /VtN,NS£G2 /V.Y.NSEGS /C.N.^ £ NSEG2 = NSEGS/2 
PARAM //C.NfSUB / V ,i>. » KMAXERR /V.N.NSEG2 /V.Y.KMAX i 
cCNO LRRORCl . kMAaErk > IF KMAX .GT. NSEGS/2 

i StT CEFACLIS FOR PARAMETEkS. 

PARAM //CfNtNOP /V . Y,NCkPkT=C1 /V t Y , LG KAD=-i S 

i calculate omega, 2-OHESA and 0KEGA**2 FROM KPS. SET DEFAULT hPS . 

PAKAMR //C.N.MPY /\,,N, OMEGA /V,Y,RPS=0.0 /C ,N,c. 233183 i 
PAKAMf\ //O.N.MPY /V,N,QMdGA2 /C,N,2.0 /V»N, OMEGA S> 

PARAMR y/C,N,MPY /V .N .OMEGASUR /V,N, OMEGA /VtN.CMcGA S 
4 CENEkATE NCkPS FLAG IF KPS IS ZERO. 

PAKAMk //C,U,E*. //V.YjkPS /C,N,C.O ////V ,N ,NUKPS i 
i MAKE SURE CUuPLEU MASSES HAVE NOT OEEN RECJESTED. 

PARAM //C.N.NUT /V.N.NJLUMP /V, Y,CCUPHASS=-i $ 

CONU LRROKC2 .NOLLMP 4 

ALTEk 21,21 * ADD SET TO OUTPUT FDR TRLG,. 

Gp 3 oEUMi,EwEXlN,uELi.-i2 / SLT.GPTT / V.N.NQGRAV 4 

CFRPNf SLTfGPTT $ 
alter 23 -t 

i SINCE MULTIPLE cONST^A^.^lTS ARE NUT ALLUDED EXECUTE uP4 iNUW SO THAI 
». MURE tRKGR ChEC^S LAN liE MADE BEFORE ELEMENT GENERAIION. 

J) AuO YS NtLDEO FOR PoF isEcCVEKY IN SSG2. 

PAKAM //C»N,MPY /v,N,NSKlP /C,N,0 /C,N,0 S 

CASEcC.GEOMA ,Ec’EX1N,GPUT,BGPGT,CSTM/K&, YS.USET ,rtSET/ V , N . LUSET / 
S,N,MPCF 1/S ,N,MPCF2/S,N, SINGLE/S ,N,OM 17/S, N, REACT/ S, N.NiK IP/ 
S,N,REPEAr/S,N,NDSET/S,N,r.GC/S,N,NDA/C,y ,ASETCUr/S,Y. AUTcSPC 3> 
PURGE GM, GMD/MPcF l/Gc,GCD/CMI r/M-S ,PSF ,CPC/S1NCLE $ 
ciiKPi'iT, oH,Ci'lUfKCfGO,ijCDiKFS,PSF,wPG,USET,YS $ 
t SUPvjRT BULK GAIA IS imuT ALLURED. 

PARaM //C, IV, not /v.n.kealdata /v.n, react $ 

C'JNU E RRGRC 3 f r> E Ac UA TA 4 

4 tXL<.ClE GPU NOR SU EHLCrS CAN Bt MADE. ADD IRE Tj OUTPUT DATA BLOCKS. 
DPG OYNA/UCS,GPL,SIL,USEr / GPLD , SI ED , USETU , T FPOGL ,ULT , PS DL , FRE, , 

TRL,.EwGYN / V , N , LUSt I/S , N , LUSEl J/V , N ,NG1 FL/S , N. NUDE T/ 

S ,N ffvbPSDL / S iN.ivOFRE / / ,N,NGNLF T / j ,N, NUT RC / V , N » I'lGE E p/ c » N , / 
StNf'VUUE 4 

i MUST HAVE EITHER FKLw OR TSTEP BULK DATA. 

PARAM / /c ,iv ,AND/V ,N ,r TEkK /V,N,NuFRL /V.N.NGTKL 4 

CUNO EKRCRcb » F I LRK 4 NO FKEw 0^ TSTEP BULK DATA. 

i ONLY FREcU'cNCY OR TSIEP IS ALLU«ED IN THE CASE CCnTKCL 
PAkmML CAbLCC //Cifv.uTl /C , N r 1 /c,N,l't / / V »N ,FREuiSFT $ 

PARAML CaSLi-C //C,N,uTI /c«N,l /C.N,3d / / V , f,, T IMESE T i 
PAkaM //CtNiMPY /V ,N,rREcTI ME /V,iN,FK£OSET /V,N»TIM£SET S 
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PARAM //CfNiNUT / V . N i f- VtiUl I /VnN.FREQTI ME i 
PAKAM //CfNfLfc /V.U.i^OFKEU /VM'J^FREuSs: T /L,N,0 
PARAM - //C.N.Lfc /V.RfNUnHE /V . w » 7 1 fliiSE T /C,N,0 i 

CCNU EKRORCo.FTERkI $ BQIH FREQ AND TSTEP IN CASE CUNTRUL DECK. 

t EPCINF oCLK CATA NGI ALtCwCO 

PARAM //CcN.NUr /V tN.tXTRAPTS /V.N.NOUE t 

CUND ERKUKC4f eXTRAPTS S 

* GtWtiUlt DATA FOR CYCT2 MODULE, 

GPCYL - GeOM4,EUUYN tUStru /CYCUD /VjN«CTYPt=RuI /ScN.NGGU S 
CCNU ERRUKCl ,NGGC i 

CHRPN7 CYCDC S 
ALTER 32 & 

S Hk£-PURGi: DATA uLuCnS FriAT WILL NOT BE OEKERATEU 
PARaM //CtN, or /VtNvNOdMl /VfN»NOMOG /VtNtNUKPS S 
FURoL blGb.MloG /NOUMl i 
POKGL H2GL,M2bAS£XG /NOMGC $ 
alter Ji> i 

}. CLNLKATL DATA BLOCKS FRLX, dlGGt MlGu, M2GG AND BASLGX. 
i LENEKATE PaKAMETcKS fkmax ang ncbasex. 

LiUMMUUl CASLlL.oGPDI ,CSTM,DI T ,FRL,MGG. t / FRLX , B1 GG. MIGG. 

M2GG,dASLAG,PU2EKG, , / V,N .NOMGG/ Y , Y, C YC 1 0/ V, Yt NSEG5/ 
V,Y,KMAX/StNirKMAX/V,Y,BXTID=-l/V, Y,BXPnO=-l/ 
V,Y,BYTID=-JL/V,Y,CVPTID=-l/V,Y.BZTID=-l/ 
Viy,BZPTlLl=-l/S*N.NUQASEX/VtN.NOFREy/V,N, OMEGA $ 

PARAML FRLX //C.N.PKESlNlE ////V,N,NCFRLX $ 


CUND LliLfKLX ,NUf-Kt.X $ 

LULiV FKLX*Fi<L i 

LABEL LbLFKLX $ 

ChKPrJT F KL I B 1 GG I M 1 GC f M 2uG » BA SE XG i 

ALTLR 42 $ 

PrthAM //L.NtADU /VfN.NuBGG /V,N,N0QM1 /C,N,0 J< RESET NUUGG. 
ALTER S2 $ 

$ t'LCEi INl BGu and KGG, 
cOND LBLllA.NOBMl $ 

PARaMR //C,N .complex // V,N,0MEGA2 / C.N.O.U / Y.N.CMPLXl 5 

PARAMk //L.N.SLb / V .N .MCM 1 .GASQ / C.N.U.O / V.N.OMEGASOR $ 
PAkAMR //C.N .CuMPLEX // V.N.MOMEGASl / L.N.u.O / Y.N.CMPLX2 S 
ALL uGG.BlGu / UoGI / C » N » ( I • 0 . 0, 01 / V.N.CMPLXl S 

Ego 1 Y oGGl.BGG $ 

ADC KCG.MIGCj / KGGi / L.NtIi,Uf0,0} / V.N.CMPLX2 S 

EGLIV kGCI.KGG $ 


LHKPN'T BcG.KGGi 

LaELL LBL I lA 

ALTLR SJ.bS S GP ^ HAS BLEN MOVcD— UP, 

ALILR ua.dfi i DPD HAS BEEN MOVEO-UP. 

ALTER iI4 i PARAM AND E,0IV LDGIl DEPENDING L)N LGKAD FUR FREO DR IRAN. 
PARAM //C ,AND/V,N,KDEKA/V.N,NCU£/V,N.N0K2PP i 

CLND LGKmUI.LGKAD i. dRANLH IN NOT FREtRLSP. 

ALTER IID i SEE aLTEk 114 CUMHENT. 

JUMP LGKAU2 i 
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LAlit'L LUAUl $ 

two i V « iPP , M2UO /NCA /t ?.P P , k 200/MGA/KP.PP r K2 DD/NOA/ KAA, MDO/ MDEMA/ 

KAA.KOD/KGEKA i 

ChKPNT K2PP, H2PP tH2PP pK.iDu, H<UD, 8?.0D tKOD, PDO $ 

LALLL LGKAU2 $ 

ALTER 117,117 t> /<UU PARAMETERS GKAU. l\(3 AND TC GKAO. 

GRAD uSElDtGM ,GU ,NAA AA , MAA , K4AA t K2PP , M2 PP, B2PP/KDD, BOD, MOO, GMD, 

000,K200,H2UO,B200/e ,y,GRAO= TRANKE5P/C,N,D1SP/C,N, DIRECT/ 
C.Y,G=0.0/C ,Y,r.3 = 0,0/C,y,N4 = 0.0/V,N,N0K2PP/V,N,NUM2PP/ 
V,N,NUB2PP/V,N,MPCPl/V,N,SINGLt;/V,N,OMIT/V,N,NOUt/V. N.NUKAGG/ 
V,N,NUbGG/V,N,l'>0f-K2/CrN,-l $ 

ALTER 118 $. SEE ALTER 11*, COMMENT. 

CUNU LCKA02,LGKAD J. BRANCH IF NOT FRELRESP. 

ALTLR 119 £ SEE ALTER ll^ COMMENT. 

JUMP LCKAL)^ F 

LABEL LCXAD3 * 

tuny b200,BDD/NuGPOr/M20O,HOO/NaSlMP/R20U,KU0/K0EK2 4 

LAttL LGKAOA S 

ALTLR 120, 123 4 
^ NEW suLurioN Logic 

4 GLNERATt TIML-DEPENULM LOAOS IF TSTtP mAS REQUESTED IN CASE CONTROL. 
LCNC LELlkLitNOTIME 4 

i LOUP THRU All subcases Fur llME-UEPENOFNT LOADS. 

//C,N,NPY /V,N,R£PEA7T /C,N,i /C,N.-1 4 

//C,N,AuD /V,N,APPFLG /C.N,l /C,N,0 4 INITIALIZE FOR SOKl. 
TRLGLUQP 4 
TkLGLGOP 4 

CASLCC, /CASEYY/C ,N,TRAN/S,N,kEPEA1 T/S ,N .NOLOCPI 4 
CASLYY 4 

//C,N,MPY /y.N.NCOL /C.N.O /C.N.l 4 

LASLYY,USETo,ULr,SLr,aGFUT,SIL,CsrM, IRL.UIT.GM u.GUO, ,ES f ,MGG/ 

, ,PLT IfPD I, ,TOL/ V,N,NCSET/S,N,POLPGU/V.N,NCCL $ 

IKL ,P0T1, ,,,,,,, , / ,POI, /V.N.APPfLG/C ,N, DYNAMICS 4 
TRL.PDI / ,P0 . /V,N,APPFLG/C,N, DYNAMICS 4 

//C,N,AUU /V,N,APPFLl. /V.N.APPFLG /C.N.l 4 APPFLG= APPFLuEI . 
IRLCUONE .RcPEATT » 

TRLGlOCP.IOC 4 
LRRCK3 $ 

TRLGUONL 4 

pcr.po, roL 4 

PD.PJT/PCEPUO 4 
PCT 4 

TOL,.,,,,, / FKCZ.ruLZ.KEChUERl ,RLCRUER2 , / 
y , Y,NSEGS/V, Y,C YC ID/i>. Y,LMA;>=-1/ V .N.FKMAX/ 
S,N,FLMAX/S,N.NrSTEPS/S,N,NCRCl/S,N,N0RC2 4 
FKlZ.FRL // RULZ.FCL $ 

FRL ,FUL,kE0K0tKl,RECRDER2 4 
JUMP CHCERL2 4 

LACLL LuLlRLl i 

4 GtNt.K«TL FRLUUENCY-ULPENOENT LOADS IF FRtUUcNCY WAS SELECTED IN CC. 


PARAM 

PAhAM 

JUMP 

LAL'EL 

CASE 

ChKPNT 

PAkAM 

TKLG 

SDRl 

SUki 

PARAM 

CCNG 

kEP! 

JUMP 

LABEL 

ChKPNT 

LUUU 

LIikI’NT 

DJMiMuU^ 


LOUI V 


3 
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l-RLG CAStAX,USETO,OLT,7Rf.,GM!J,GC!> a)i T t / RPF tPSFt PUFtFOL. PHFftun / 
CtN,t)Ift£L7/V.N,i-!>EUV/C,,N,FR£0 6 

CCNJ LbLHiLXl.NUFkLX S ZERO CUT LOAD COLLMK'S IF FKLX WAS GENERATED. 
MPYAO FPF.PUZEkO. / PPFX /C.N.O S 

LULIV PPFX.PPF i 

LALtL LELFRLXl S 

i lORfJ NEW LOAUS. 

CLNl LbLFRLl .NObASEX S 

MPVAC M2GC.bASLXG, / M2BASEXG /C.N.O S 

AOC PPF,M2bASEXG / PPFl /C , N . ( 1 . 0 .0. 0 ) /C .N « (-1 .0 .0.0^ $ 

EJlilV PPFl.PPF i 

CCNO LBLbASEl. NOSET i 

SSC2 USETO.GMD .VS.KFS.GOO. .PPF / .PODUMl . PSFl ,POFi t 

tWLiV PSFl.PSF // PDF1.PDF $ 

LAEEL LULBaSEI $ 

LABEL LBLFRLi 5> 

tOLiV PPF,?!jF/NDSL I i 

CHKPM PPF.PSF.POF ,FOL i> 

i LuAuS AKt FKEOUENLV-uEPENUENT 

i PLPFCKI-t CYCLIC TRANSFLRMATi GN ON LOADS IF CVCIC=£1. 

PAHAML Pbf- //C.N.TKAILcK yC.N.l /y.N.PDFCOLS $ 
i LALCLLmIE IHE NUMbLR OF LOADS FCR CYCI0=-1. 

PAKAR //C.N.DiV /V.N.NLLAD /Y.N.PDFCOLS /V.N.FxMAX $ NlOAU = NF/FKKAX 
tv,Li V PCF,PXF/CYC IL. i 

COND LaLPDONE.CYLlU t 

<. LAlCULAIL the NOHULK CF LOADS FCR CYCI0=1. 

PARAM //C.N.DiY /V.I^.NlCAD /V.N.POFCCLS /V.Y.NSEGS i NLOAu = NF/NSEGS 
CYU I PDF / PXF.GLYCFl /V.N.CIYPE /C.N.FORE / V , Y , NSEGS=- 1 / 
V.Y,XKAX=-l / V.N.NLOaD /S.N.NeGO 4> 

CGNU LRRORlI.NUGO 5 

CFKPNI PxF i 

JUPP LuLPbONE & 

LAEEL LoLFkLZ i 

I. LOADS AKE TIME-DEPLNDENT 

PAKAM //C.N, Not /V.N.NOTLYCIO /V.Y.CYCIC $ 

^ tisA.'JCH DEPENDING oN VALDL OF C Yu I C 

CCMJ LGLTKLZ.NOToYL iC i 

S CYU0=-1 

LDOlV PDT .PDTKZ 1/NORUl !■ 

CONU LBLnO lA .NGRCl S> 

MFYAJ POr.REOkJERl , / PDTRZi / C.N.D i 

LAEEL LbLkCIA J. 

LYCri PDTKZl / PXTkZ I ,G0YCF2 / V .N . CT YPE/C . N , FORE / V . N.NT S T c PS/ 
V.y.LMAX/V.N.FxMAX/S.W.NCGQ i 
CCND LRRuHC 1 . NOGU t 

ohlsPNl PXTKZl S 

EUUIY PXTK/l.PXFZi/NokOZ # 

CCND LtJLROZA.NOKCX 1 

HPYAD PXTkZ 1 .kLOKULkZ . / PXFZl /C.N.O L , . 

LAtlEL LBLKUZA $ 
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N A S T i< A N 


OR'amAl TASS ig’ 
OF POOR QUALITV 


t M £ C U T I V £ 


CONTROL 


DECK 


ECHO 


i 


EJUIV PXPZl,PXFi $ 

ChKPNT PXFl i 

JUMP LBLTKL3 i ‘ 

LAEEL LULTRL2 i 

$ CYcIG = Cl 

MPYAU PUT, KtOKUERl, / PDTRZ2 / CtN.O S 

CYCTl PDTKZ2 /PXTk 2Z,GCYCF3 / Y. N.CT YPt/C .N ,f DKE/V « N, NTSTEPS/ V, Y t LMAX/ 
V.Y,NS£GS/S,K,NCGC S 
CGNO . ERkCK C I « NGGC $ 

ChKPKT PXTRZ2 i 

LQUIV PXTKZZf PXTR2/NUR02 t 

.’.UNO LBLkOiO.NURUZ i 

MPYAU PXTk22,hLlJl<UbR2, / PX1R2 /C.N.O S 

LALLl LbLk02B 1. 

CYCTl rXTK2 / PX(-Z2,GCYCF4 / V.N,C TYPE/C ,N , FORE/ V » Y . NSEGS/ VtY ,KMAX / 
V,N«FLMAX/S,N,NUGC $ 

CGNU EcvRuRCl^NGGC $ 

EGUIV PXF22.PXF 1 $ 

CHkPNT PXFl t 

LAUtL 1.DLTRL3 i 

> I irtfc-uEPLNUENT LOADS ARE REAL. MAKE LOADS COMPLEX To COKRESPONu 
$ TU FktGUEMCY UtPcNULNT LGAOS. ALSC SUR2 EXPECTS LuACS TO BE COMPLEX 
* IN F«cOhtSP TYPE Problems. 

LuPir PXFl / PXF2 i CCNVEKT REAL PXFl TC LOMPLEX P'^F. 

AUf PXF1.PXE2 / PX» / CtN,(O.S,l.O) / C.N. (0.5 .-1 .0 ) 5 

i LEFIM NLOAD FCk CYoT2. 

//Cf.ifAUl) /V.N.iNLOAD /V.N.FlMAX /C.W.O $ 

LBLPDGNE i 

//otN.AUui ZV.N.KINOEX /V.Y.KMIN=U /C.N.U 


PAHAM 
LABEL 
PA RAM 
i 

i IN n lAL ize 

i 

PARAM 
lONu 


NLCAU = FLMAX 
. INTITIALIZE KINUhX, 


UXVf IF KMIN IS NOT ZERO. 

/V.Y.ivMIN /L.N.-l 


PAKAM 

JUMP 

LABEL 

CYCT2 

oLNu 

4CC 

CYCT2 

oCNu 

PAkAM 

KEPT 

LABEL 

i 

JUMP 

label 


//C.NfAUU /V,N,Ki'lINL 
WUKM. INL .KMINl 5. 

//C.NfAUU /V.N.KMINV /C.N.O /C.N.O i 
KMIiMI.OUP I 
KM I NX OOP $ 

lYCOU, , ,P xE , , /..PKFZ,, / C.N.FXRE/V.Y.NSEGB/ 
V,N,KMINV/V.N,CVCSEU/V,N,NLCAU/S .N.NCGO $ 
ERRGRCl.NOGC i 

PkFZ. / UKVFZ / O.N.(O.C.O.O) S 

CYCLj,, .lKVF^.,. /..UXVF.. / C.N.BACK/V.y.NSEGS/ 
VfN#KMIi\V/Vfi'^fLYCSEU/V*NfNLCAU/S iNtNCCC $ 

EkRORL 1 tl.'JGC i 

KMINV /v.N.KMINV /C.N.l $ 


//C.N.AUU /V,N, 
KM INLOOP , KMINL 
/40KM1NL i 


lOPCYC 
rCPC YC 


L0(’P on K INDEX 
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N A S T R A N 


0?.KHJh quauty 

E X E C U i I M Q CONTROL 


DEC 


ECHO 


I 

a 


CCNU 
PHTHAKM 

label 

ever 2 

CONU 
CHKRNT 

pakam 

i 

CYCT2 

PAkA,-l 
CC N D 
LhkPNT 
A SOLUTION 


nokprt.nokpkt a 

//CfN,0 /CtN.KlNOex 5 
NCKPRT i 

CYCuUiKDUfMOU, , , /KRKf » HKKF i . . /C i N f FCUE/ V tV »NSEGS f 
VfNtMNOEX/V.N,CirCSEy=“I/V»N.NLOAD/S ,N,NOGC % 
tRRURLl.NOGC A . . 

K.KKF,HKKF $ 

//C.N.SYST //C,;..5b /C,N.2 $ METHOD 3T IN CYCT2 PRODUCES 

LNOERFLUWS FOR PXF. USE METHOO 2, 
LrCUU,iiCU,,PXF, , /DKKF f ,PKF , * / C «N .FURE/V t Y » NSEGS/ 

V tN , K INUt X/ V . N .C YC S£ 0/ V , N , NLG AU/ S * N , NOGC A 

//C.NfSYSr //C,N,53 /LtNtO 5 RESET HPYAD HETMUU CCNTRCL. 

EKKL'RC 1 « NQGC A 
BkXF.PKi- $ 


FRRU2 

CHKPNT 

LYCT2 

CCNU 

LhkPM 

P AR AM 

PARAH 

LCND 

i;tPT 

JUMP 

LAbEL 

hwUl V 

CHKPNT 

CCNJ 

cYCTl 


KKKF ,OKKF .MKKF . f PKF .FOL / UKVF /C t N . 0. 0/C . N.O .0/ C. Nt -1 . 0 5 
UKVF $ 

CYCUJ,, ,oKVF, f /s.UXVF,, /C , N »BACK/V . Y . MSECS/ V .N,K INDEX/ 
V.N.CYCSEJ/VfNfNLCAD/SfNfNCGC $ 

LKRORCi,NOGG A 
UXVF A 

//C,NtAL.D /VfN.KINDEX/VtN.KINDEX/C.N.l A KINDEX = KINDtX 


C 1 


//C,N,SLO /VfNiDUNt / V, 
LCYCZ.DONt I IF KINDEX 
TLPLYLtlCC $ • 

ERROR 3 A 
cCYC2 A 

UXVHtUDVF / CYCiu A 
UUVF A 

LCYC3,LYCIU a if lycio 


Y.KMAX / V.N.KlNOEX 
.GI. RHAX THEN EXIT 


.GE. 0 THEi'. TRANSFORM To PHYSICAL. 


/ V.N.CTYPE/C.N.HACK/ V, Y.NSEGS/ V. Y.KMAX/ 


UXVI / UDVF.oLYOal 
VfN.NLOAJ A 
CHKPNI UCVF A 

LABEL LCYL3 A 

oCND LuLTRLA.NUTIKh A 

cGGlV PXFtPUF/ / CYCIC S 

CCND LLYC4,CYlIU 1 iF L Yo 1 0 .GE. 0 THEN TRANSFORM TO PHYSICAL. 

CYCTl PXF / PCF2,ooYCr)2 / V , N ,C TYPE/C . N . BACK/ V , Y , NS EGS/V , Y , KMAX/ 

V.NiNLOAD A 
LABEL LCYC^ A 

A IF loads wLlsE r IMt-uEPENUENT THEN RECOVER PPF AND PSF FROM PXF. 

SDRl USETU.tPDF/,, .gUD.oHD., .. / PPF2,. /L,N,1 /C t N .DYN AM I CS i 

SSC2 JSeTD.GMD,Y£,^l■S,GCU,,PPF2 / . PCDUM . PSF2 , PLDUM i 

EOUIV PPF/.,PPH // PSFZ.PSF A 

uHKPNT HPF.PSF A 

LAtcLLbLIRL4A 

Alter 1/4.124 A USE fOL INSTEAD OF PPF TO GET OUTPUT FREuJENCY LIST. 
VDR CASlXX.twDYN .USETDfUDVF ,FGL , X YCDU . /OUDVCi ,/ C . N , FRE ORES P/C, N, 

D IRECI/S fN.NUSOR J 2/S NCU /S , N .NOP/C , N ,0 A 

alter I4C, 1'»0 a use fUL INSTEAD OF PPF TEi oET OUTPUT FkEwJENCY LIST. 
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E C H 0 






N A S T R A N 


E X c C U T 


C.'J.'GJrJAL PAG^ [g 
a" POOR QUALir/ 

CONTROL 




DECK 




TV 

t 3 


fr 

ii 

il 


buK^i CASEXXiLbTM.MRT idIT,£Qi)YfwS£i.L>c t cEGf'DP.FLLtWFC.UPVCf EST iXYCOQt 
PPF/UPPLl lOGPLl »aURVC.l« OeCCI «Ci:C Ci .PUPVa/C.N.rREORbSP/ 
StN.NGSUkrZ i 

160 li ' ADO LAdfcL Fu« ERkQRSj 
ERRORS $ 

103,166 $ REMOVE ERRCnl ANU EHHOR^t;, 

168 $ FURCEd VlbRAllCN ERRORS 

ERRORCl $ CMLOK NSEGS, RMAX AND OTHER CYCLIC OATA. 

//C,H,-7 /C.iwlYCSTATICS S 

ALLUtvLO. 


'<! 


ALTER 
LAEEL 
ALTER 
ALTtR 
LAEEL 
PrTPAkM 
LABEL 
PklPAKM 
JUMP 
label 

PRTPAmM 
LABEL 
PR I P ARM 
JUPP 
label 

PR TP ARM 
PR TP ARM 
JUPP 
label 

P.< 1 P ARM 
PRTPAkM 
JUMP 

enlalt tk 


$ CGuPLEd HASS NOT 
/C.V.CCiUPHASS S 


ERR0RC2 
//C.N.O 
FINIS t 

CkR0RC3 i SUPOkT oUlR DATA NOT ALLCaEO, 
//C,N,-6 /C ,(/,C YUSTATICS $ 

ERRORCA i> EPuiNT BULK DATA 
//l»N,0 /CvN,NuUE i 
FINIS $ 

EKKURL'j 
//C.iJ.O 
//C,N,0 
FINIS $ 

ERkL’kLtj 
//C,N,0 
//C,N,0 
FINIS I 


i NEITHER FREW 
/C,N,N0FRL $ 
/CtN.NOTRL $ 

i BOTH rKEU ANU 
/C,N,NQFREU i 
/C,i\,NuTIM£ 5 


NCr ALLOWED. 


Ok ISIEP WERE IN bulk DATA DECK. 


TSTEP WERE SELECTED IN CASE CUNTRdL. 


n 
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0;iiG5?;Al 

01- POOH QUALiTif 


FORCED VIBRATION AMALYSIS OF ROTATING CYCLIC STRUCTURES 
2. A. 2 DHAP Sequence for Forced V ibration A ndysis of Rotating Cyclic Structures 


TICNS n EFFECT GO ERR^-2 NOLI ST NOOECK NQREF NOOSCAR 

1 BEGIN N0,8 FOKCEO VIBRATICNS OF ROTATING CYCUC STRUCTURES - SERIES R 

2 PRLCHK ALL S 


3 Ml E 

4 f IL E 

5 CONt) 

6 CU'4 0 

7 P A? AM 
fi co^^n 
9 t» AM 

10 PAR AM 

11 CO'JD 

12 ? A^ AM 

13 PAR AMR 
lA PARAMR 

15 PARAMR 

16 PARAMR 
I r 3 AR AH 

18 COYU 

19 PAR AM 

20 CPI 


KGGX = TAPE/KGG=rAPE /GUU = SA VE/GMD=SAVE/MDa=SAVE/80U=SAV E 6 
UXVF=APPEND/P0T=APPENL)/PD=APP£N0 $ 

ERRLRCLfNSEGS 5 IF USER HAS NCT SPECIFIED NSEGS. 

ERR'JRC IfKMAX $ IF USER HAS NCT SPECIFIED KMAX. 

//C,N,E3 /Vf N»CYCICERP. /V,Y,CVC10=0 /C,H,0 {■ 

ERRORCl.CYCtOERR S IF USER HAS NOT SPECIFIED CYC 10. 
//CrNtUlV /VfN,NSEG2 /VfY, NSEGS /C.N,2 $ NSFG2 = NSEGS/2 
//C,N,SUD /V,N»KMAXERR /V,N,NSEG2 /V«Y»KMAX S 
ERRURCUKMAXERR ^ IP KMAX .GT. NSEGS/2 
//CfNfNUP /V,Y,NOKPRT*a /V, Y , LGKAO=-l i 
//CrNfMPY /YiN, OMEGA / V t Y i PP S=C. 0 /C tN ,6 , 283 18 5 t 
//CtNfMPY /V,N,CIMEGA2 /C,N,2.0 /V,N,CMEGA $ 

//CtNfHPY /V .NfCMEGAScP /V,N, OMEGA /V,N,CM£GA $ 

//C»N,EQ //VfY,RPS /C»N,0.0 ////V,N,NORPS $ 

//C,N,NOT /VfNfNOLUMF /V, Y ,C0UPMASS=-1 $ 

ERRORC 2, NOLUMP i 
//AMPYf/CAKDNO/U/C S> 

GEUMltCECM2, /GPL ,EUEX 1 N ,GPD T » C S T M , BGPDT tS I L/ S • Nt LUSET / S, N, 
NUGPD7 R 


PLTTRAN BGPUTf SIL/BuPDP,SiP/LLScT/StN,LUSEP i 

22 PURGE USET, GM.GOf KAA ,8 A A ,MAA , KAAA tKF S , PSF ,QPC ,E S T t ECT , PLTSETX , PLT P AR , 
GPSETStELSET S/NUGPU I i 




Oi- POOR QUALITY 


LE\/EL 2.0 NASTRAN OMAP COMPILER - SCU.'i'!CE LISTING 


m 


I 


r 

r 

[ 


23 

CONL) 

24 

GP2 

2£ 

PAR AML 

26 

PUR GE 

27 

CO'IO 

28 

PLTSEI 

29 

PPTMSG 

30 

P AR AM 

31 

P AR AM 

32 

CU'JO 

33 

PL3T 

34 

PRTMSG 

35 

LABEL 

36 

GP3 

37 

CH<PN T 

38 

1 A1 

39 

^URGE 

40 

PAR AM 

41 

GP4 


r 

^2 PU^GE 
’ 43 CfKPNT 


LBL5,N0GP0T f> 

GEQM 2,EQEXIN/ECf 6 
PCOB//«PR£S<=////NOPCOB S 
PLTS£rx,PLTPAR,bPSfcTSeELSers/NQPCOB S 
PUNUPCOE $ 

PCOH,EqEXIN,ECT/PLTS£rX.PLTPAR,GPS£TS,ELSETS/S .N.NSIL/ S,N, 
JUHPPLOT*--! S 

PLTSETX// S . 

//»MPY«=/PLTFLG/l/l $ 

//«MPY*/PPIL6/0/0 $ 

PljJUMPPLOT S 

PLTPAR,GPSeTS,GLSETS,CASECC,BGPDTtEQEXlN,SIL,,ECT ,,/PLfJTKi/ 
NSIL/LUSET/S,NtJUMPPLOT/S,N,PLTFLG/S>NcPFlLE i 

PL'JTXl//$ 

?1 i 

GEUM3,EQEXIN ,GEfJM2 / SLT.GPIT / V,N,NOGXAV $ 

SLTfGPTT $ 

ECTf bPT .EGP0 7 .SIL.OPTT ,CSTM/EST,GEI ,GPEC T , ,/ LUSET/S , N. NCS I MP = 
-1/ l/S,N,N{.IGENL=-l/£tNtGENEL V 

X4GG,GPSTf(iGPSr,MGG ,BGG ,K4NN ,K4FF ,K4AA ,.1NA tMFF f MAA, BNN, HFFi UAA 
KGGX/NCJSIMP/CGPST /GcNL-L i 

//CfNfMPY /V,N,NSKIP /CfN,0 /C,N,0 i 

CASECC ,GE0M4 ,EUEXI NfGPOT.GGPIJT ,CSTP/RG,VS i OS ET , AS tT/ V • N, LU S ET/ 
Sf N, MPCF 1/S , X,MPCF2/S,N, S INGLE /StN ,DMI F/ S ,N » RE AC T/S t N, NSKI P/ 
S.Nf KEPEAT/S ,N,NUSF I/S ,N f NOL/ 5 » N , NOA/C , Y , ASETOUT /S f Y t AUT US PC 4 


GMf GMD/MPCM/GO,GGU/OMir/KFS,PSF ,QPC/SINGLE S 
GM^oMOiRGtGU ,GC»tKF StPSF ,CPC -USE T,yS i- 
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ORKlif^IAL 13 

OF POOR QliAinY 

% 


LEVEL 2.0 NASTRAN DMA F COMP J Lci< - SOURCE LISTING 


•'*4 

P AH AH 

//CfNsNOT /V fN,F£ACnAT.A 

/V^HjP-EACT S . 

45 

CGNO 

ERHURC3,REAC0ATA S 


*♦6 

UPO 

OYS'AMICSiGPLfSIl.tUSET / 

GP LD, SI LD fUSETD»TP POOL, QLT, PS DLt HRl. 


TRL»fEQDYN / V ,N f LUSET/S,N t LOSE TD/V ,N. KDTFL/S , Nt NODLT/ 
S,N,NGPSOL/S»N,NOFRL/V,NrNONLET/S,N,NOT(a/V»N,NOEED/C.N»/ 


4 7 

P A-l AM 

SiiNtNOUE $ 

//C,N,ANG/V,N,FrERR /V,NtNCPRL /V.NfNUTRL i 


4« 

CUNO 

ERR0RC5fFTERR $ NC FREQ CR TSTEP BULK DATA. 


49 

PAH AML 

CASECC //CfN,UT! /C,N,i /C,N,14 //V »N tERECSET $ 


‘>'J 

P ARAML 

CASeCC //CjNfDTl /CtNa /CtNfiB /.' V , TI MESET & 


51 

PAR AM 

//C,N,MPV /V,N, FREQTIME /V,N,FREQSET /V»N,TIMESET 

& 

52 

P AH AM 

//C,N,NOT /VfN.FTERRl /V,N,F RtGTI ME S 


53 

P A.H AM 

//CjNtLE /y.NjNOFPEQ /V.N.FREQSET /C,Nf<3 % 


54 

PAH AM 

//CiNiLL /VfN.NGTIPE /ViNfTIMESET /CtNtO $ 


55 

C'JMD 

ERRURC6,FTERR1 i BOTH FKEG AND TSTEP IN CASF CU.NTROL DECK, 

56 

P AH AH 

//Cf.NtNOT /V,Nff;XTRAPTS /V,N,NCUE i 


6 7 

CUND 

ERKURC4,EXTH AP TS S 


58 

GPCYC 

G£JM4,£QCYNf USETD /CYCDO / V ,N ,C TYPE-ROT /SfN.NCGG 

$ 

59 

CUN D 

ERRGRCl.NOGU i 


bO 

CIKPM r 

CYCOD i 


61 

C'NO 

LBLl.NOSIMP $ 


62 

P AP. AM 

//*AOO*/NUKGGX/i/C i 


63 

PAR AM 

//=»AUO*/NOKGG/ I/O $ 


64 

P AH AM 

//♦AOCi*/NUBGG=-l/i/C S 


6 5 

P AH AM 

//=»A0D<=/N0K4GG/l/0. $ 


66 

EMG 

EST, CSTM ,MP r ,DI T ,GE(JM2 r/KLLM ,KDICT ,HELM, MD ICT , DbLM 

,BOICT/ 


NtN'JKGGX/S tN sriOr(GG/S»N,NQtlGu/S ,N,NCK4GC//C,y tCOUPMASS/C» Yt 
CPBAR/C, Y,CPP.UD/C ,V,CPUUADi /C , Y ,CPGUAD2/ C , Y r CPTR lAl/C, Y, 
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uitiGiNAL r:i’ 
OF POOR QUALITY 


LEVEL 2.0 NASTRAN OHAF COf-iP» LgK — SOURCE LISTING 


^ 67 CO'JD 
^ 68 6SA 
69 LABEL 
P 70 PAR AM 
« Jl PURGE 
^ 72 PURGE 
P 73 CONO 

Q . 

74 EMA 
I 75 LABEL 
_ 76 DUMMOOl 


77 

P AR AM L 

i 

CUNO 

79 

0 

EQU IV 

U ISO 

LABEL 

1 

CH<PN T 

y 

fl2 

CONU 

I ” 

EMA 

84 

LAB EL 

I 85 

CUNIU 


LMA 

8 7 

LAB EL 

I Ofi 

PUR GE 


CPTR IA 2 /C,VjCP 7 U 3 E/C»V,CfHlOPLT/C»VrCPTRPLT/C,YsCPTRBSC $ 
lblkggx.nokggk $ 

CPECTfKOICT,KELM/KGGX,G?ST $ 

LBLKGGX $ 

//C»^',QP /VtNfNOBMl /VsN»NOHGG /VjNjNOKPS S 
OlGGf/UGG /NGBMl $ 

H2GG,H23ASEXG /NONGG S 
LBLMGGjNCMGG $ 

GPECT,MUICT»MELH/MGG,/-l/CrY,6TMASS=l.O J 

LBLMGG $ 

.•?2u,G,BASEXG,PUZERU, , /V , A ,KUMGG/ V , V ,C YCI 0/ V , Y .RR EGS/ 

V,Y,r.MAX/S,N,FKMAX/V,V, 8 XTI 0 ^-l/V y|HXPriD=-l/ 
V,Y.BYTID=-l/V,Y,BYPTIO=-i/V,Y,BzflO=-l/ 

Y t Yf BZP I ID=- I/SiN »NG8A SEX / V»N f NCFREG/V »N f CPEGA S 

FRLX //C. iN fPRESENCE ////VtN»NOFRLX $ 

LBLFRLX.NOFRLX i 
FRLX.FRL J 
LBLFRLX $ 

FRL.a lGGiMiGGfM2GG,f?AS£XG $ 

LfiLBGGfNOBGG i 
GPECT,B01CTtBELM/BGG, i 
LBLOGG 1. 

LBLK 4GG,NaK4GG i 

GPEC T,KDICT f KELM/K4GG t/NUKAGG $ 

LBLK4GG t . 

MNN,MFE,MAA/(\U.VGG i 
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°®GW,U 


LEVEL 2.0 NASTRAN OtiAP CU?^?IL£R - SOur.CE LI5VING 


7U 

i 


PAH AM //C»N,ACiD /V;M*MCoGG /V»N#NOCKl /CjMfU S RESfeT NaOGG, 

PURGE BNNf OFErBAA/liUeCG 5 

tOMO LBLlfGPOPNT £ 


CONU EKROR<>,NOMGG £ 

GPWG BGPDP,CSTH,EQEXINfMGG/CGPhG/VtY»GP.DPNT*^-l/CfY»WT,iASS $ 


□FP aGPWG,,t,,//S»N,CARUMO i 

LAO EL LOLl £ 

tyuiv KGCXfKGG/NQGENL S 


CO^JU 


LBLllfNOGENL $ 


SMA3 GEI»KGGX/KGG/LUSE7/ACGENL/N0SI HP £ 

LAO EL LBLli £ 

CD'JG LBL11A,N03MI S 

PAH AMP //C,N, COMPLEX // V,N,QMEGA2 / C,N,0.0 / V,M,CMPLX1 £ 

PARAMR //C»M,SUB / VtN,MOHEGASQ / CtNfO.O / V.N.CMEGASUR £ 
PARAMR //C,N, COMPLEX // V,NtMCMcGASC / C,N,0.0 / V,N,CMPLX? £ 

ADO OGG.BIGG / OGGI / C,N,(i. 0,0.0) / V.N.CMPLXl £ 

EQUIV BGGl.BGG £ 

AUIJ KGG.MICG / KGGl / C ,N, (1. C,0.0) / V,N,CMPLX2 £ 

EOU IV KGGl.KGC S 

CHXPNT BGGtKGG 5 


I 


LAD EL LHLllA 


COMO LBLA.GENEL £ 


CONI) 


LBL^jNOSIMP i 


PAH AM 


//<'E0<=/GPSPr LG/AUTOSPC/0 $ 


COMO 


LUL4,GPSPFLG $ 





LE\/EL 2.0 NASTRAN OMAP C{jr-!Pn.C-f^ ~ SOURCE LISTING 


114 

GPSP 

GPLiGPST,US£7fS[L/0GPSl/SvNr!ICGPST S 

115 

COND 

LBL4sNrGPST $ 

116 

3FP 

OGPST,,,,,//S,N,CAi<ONC S 

117 

LA8EL 

LBL4 $ 

118 

EQJ IV 

KGG»KNN/«PCf- i/MGG,MNN/riPCFU BGG ,BNN/MPCF 1/K4GG, K4NN/MPCFI $ 

119 

COND 

LBL2fMPCFl $ 

120 

MCEl 

U£ET,RG/GM i 

121 

'1CE2 

USET,GM,KGGf.MGG.BGGTK4GG/KNN,MNN,B.NNtK4NN $ 

122 

LABEL 

L BL2 $ 

123 

EQJ IV 

KNN,.KFF/SINGLE/HNN,.«FF/SINGLE/BNNtBFF/SI NGLE/K4 NNtK4FF/S INGLE 

124 

LONG 

LBL3. SINGLE $ 

125 

SCE 1 

USET.KNN.MNN ,BNN,K4NN/KFF ,KFS , (HFF ,BFFvK'»FF S 

126 

LABEL 

LBL3 i 

12 7 

EQUIV 

KFF,KAA/OMIT t 

128 

EQU IV 

HFF,MAA/GMIT S 

129 

EQJ IV 

BFF,BAA/0H1T i 

130 

EQJ IV . 

K4FF ,K4AA/0MI T $ 

U1 

CO'JO 

LBLSfOMIT $ 

132 

SM3 1 

USET,KFF ,, ,/GC',KAA,KCC,LCO,« ,, , $ 

133 

CCiNO 

LBLMfNOMGG $ 

134 

iMP 2 

USET,GO,PFF/MAA $ 

13 5 

LABEL 

LBLM $ 

13 6 

CUN J 

LBLB rNQBGG S 

137 

SMP 2 

USET ,G0 ,BFF/OAA S. 

138 

label 

LBLB t 
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1 

CQND 

lAO 

E? 

SM’2 

i 1^1 

LABEL 

1? 

L-QJ IV 

P»3 

PAS AM 

1 l-VA 

PAP AM 


BMG 

I 

l'*6 

PAR AM 

1 WT 

PUS GE 

1^9 

r 

CO'J 0 

1 l',9 

MTSXIN 

r ISO 

LABEL 

1 

LSI 

JUMP 

[ 152 

LABEL 

1S3 

f. 

PURGE 

f 154 

CASE 

I- 

155 

MTRX IN 

f 


' 156 

PAS AM 

' 157 

P AP AM 

1 

150 

EQJ IV 

159 

A 00 5 

160 

CUNO 

161 

Tr\ SP 


0? POOR QUAin;^ 


LEVEL 2.0 NASTRAN OMAP COMPILER - SOURCE LISTING 


LeL5,NQK4GG S 
USET ,GO,KAFP /K^AA S 
LBL5 f. 

GG,C.OO/NGUE/GH,GHD/NCLc 5 
//OAUO* /NEVER /I/O $ 

//«=MPy*/REPEATF/-l/l $ 

MATPar.L,BGPn7,fcOEXIN,CSTP/3DFCCL/S,N,NCKHFL/S,N, NOAbFL/ S,Nt 
MFACT i 

//«'AND*/N0FL/NGA8FL/M0KaFL 6 
•<BFL/NOKBFL/ ABFL/NCAflFL i 
LBLFL3.NCFL $ 

,BUP00L,EQDYN, t/ABFL,KBFL,/LUSETO/S,N,NOABFL/S ,N,N0KBFL/0 i 
LBLFL3 i ; 

LBL13 i ' 

LBL 13 S 

OUOVC l.UL0VC2,XYPLTFA,UPPCl ,CGPCl ,CUPVC1 ,CESC1 ,OEFCl , UPPC"' • 
QQPC2,0UPVC2,0ESC2,G£FC2, XyPLTrtPSOF ,AUTC , X YPLTR t K2PP,M2PP, 

Q2PP fK2J0,M200fU20U /NEVER S 

CAStCCtPSDL/CASCXX/^FREO^/S.N.REPEAlF/S.N, NOLOOP i . 

CASEXX.MATPOOL.EyOYN, ,TFPCCL/K20PP,P2DPP.e2PP/LUSETD/5 »N, 
N0K20PP/S,N, NGM20PP/S,NiNC82PP i 

//*AN0>»/N0M2I'P/NC'ABFL/NCM20PP $ 

//♦AND-i'/NOKZPP/NaFL /NCK2DPP t 

M2UPP, M2PP/NUABFL S 

ABFL ,i<0FL,K2nPP, , /K2PP/(-l .0,0.0) J. 

L!iLFL2,NCAUFl. i 

ABFL/AOFLT i 
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QRlGtMAL FAGS :u 
OF POOR QUALITY 


LEVEL 2.0 NASTRAN DfiAP CGKtP* LEU - SCURCE LISTING 



i 

A OD 

If. 3 

LABEL 

1 16^ 

PAR AH 

B 

PAR AH 

lf>(> 

PAP AM 

I 

PURGE 

If. 8 

PAR AM 

B 169 

CO'f D 

3 

EQJ IV 

171 

JUMP 

y 172 

LAB EL 

i 

EQJ IV 

PI 

CHKPNT 

175 

LABEL 

1 

CONO 

177 

GKAO 

^ 178 

LABEL 

1 

COHO 

130 

EQU IV 

• 

cc 

. 

JUMP . 

■ 132 

LABEL 

« 103 

EOU IV 



ABFLTf M2Df»P/.*U'?p/MFACI 5 
LQLFL2 $ 

//«'AND<‘/BDEBA/N0UE/hCB2PP $ 

//^ANOO/KDEK^/NOGENL/NOSIMP S 
//<=AND<‘/MDEMA/NaUe/NCM2RP i 
K20D/N0K2PP/H2DD/N0N2PP/B20D/N0B2PP 6 
//C»N*AND/V,NfK0EKA/VtN,N0UE/V»N*NCK2PP t 
LGKADl.LGKAO $ BRANCH IN NCT FPEQRESP. 

M2P? ,M2D0/NQA/Q2PP»B2DD/NLiA/K2PPrK2D0/N0A/MAA,MDD/MUEKA/BAA. 
800/BDEBA S 

LGKAD2 S 

LGKAOl $ 

r'(2PP*M2[)O/NUA/t52PP,02DO/NOA/K2PP,K2DD/NQA/MAAf MDD/MOEMA/ 

KAAt KDU/KDE kA $ 

K2PP fM2PP.32PP,K2D0»M2DDf B2DD,Kt;U,MDD t 
LGKA02 > 

LBLlR.NOGPOr $ 

U SE TU, GM, GO, KA A. B A A, MA A, AA, K2PP, M2 PPf 02 PP/ KDD, 8DD, HOD, GMO, 
GOO, K2DD ,M20D,B20n/C , Y ,GKAO= TR ANRES P/C ,N,OISP/C, N, DIRECT/ 
C,y,G=0 .0/C , Y,W3 = 0.C/C,Y, h'^ = 0.0/V,N ,N0K2PP/V ,N,N0M2PP/ 
V,N,N002PP/V,N,MPCFl/V,N, Si NGLE / V , N ,UM I T / V ,N, NOJ E/V , N , N0K4GG/ 
V,N,HGHGG/V,N,KDEK2/C,N,-l i 

LBLIO $ 

LGKAD3,LGKAJ $ BRANCH IF NOT FREQRESP. 

82D0 ,B0D/N03GG/ H2DD , MOO/NOSI MP/ K2 00 , KDD/ K0EK2 i 
LGKA04 S 
LGKA0 3 ;» 

B2DD,BDU/NOGPDT/M2DO,HDO/NGSIMP/K2DO,KCO/KDEK2 $ 
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LEVEL 2.0 NASTRAN DHAP COMPILER - SOURCE LISTING 


a 

LABEL 

fn 

CUMD 

y ii!6 

3 A.M 

18 7 

P AR AM 

y 

IF a 

JUMP 

p IS 9 

LAHEL 

190 

CAS E 

i 191 

CH<PMT 

n I'^Z 

PAR AM 

19 3 

TRLG 

1 19 4 

SDR 1 

3 

SDRl 

19fa 

P AR AM 

I 

COMD 


R£PT 

a 19 9 

JUMP 

i 

LABEL 

201 

CHKPNT 

1 202 

LQJ IV 

203 

CM<.PNT 

□ 204 

UUMH002 

20 5 

n 

LOU IV 

y 

CH<PNT 

20 . 

J UM P 


LGKADA $ 

L 8 LTP.Ll 9 .MOTi HE S 

//CfNjKPV /VfNjREPEATT /C»Na /C,Nr~l $ 

//C.MfAJD /V,N,APP!-LG /C,N,L /CfN,0 5 INITIALIZE FOR SDRl. 
TRLGLCQP $ 

7RLGLGUP $> 

CAS£CC»/CASEVY/C.NtTRAN/S»N»Rl:PEATT/S#NfNCLtJOPl $ 

CASEYY 4 

//C,rj,MPY /y,N,NCOL /C,NfO /C,Ntl $ 

CASL YYtUSE TO ,DL T , SL T»BGPDT ,S IL f CSTM ,TRL , 01 T ,GHDt GODt , EST.MGG/ 
, ,POTl,POi, , TCL/ VtN,NCSET/S»N,PnEPDC/V,N,NCOL i 

TRLfPDTltt,, ,c.f » / »PDT, / V ,N , APPFLG/C , N ,UY.NAMI CS S 

TRL»PD1 / »PD , /V,N,APPFLG/C,N,DYRAHICS $ 

1 

//C,N,AOD /V.N.APPFLG /VtN.APPFLG /C,N,1 $ APPF LG= APPFLGC 1 . 
TRLGDONEfREPEATT $ 

TRLGLOGP.IOO $ 

ERRORS 5 
TRLGOO.ME % 

POTfPO, TOL S 
POfPUT/PDEPOG i 
PDT $ 

TOL,t»»».f / FRLZ.FCLZ,RECRlJERlfRECRDE.R2,,,, / 

V,y, NSLGS/ V, Y ,C YC IC/S,Y,LMAX---1/ V,NtFKMAX/ 
S,N,FLMAX/SfN,NTSTL'P5/SiNfNCRQl/S ,N,N0R02 4 

FRLZ ,FRL // FOLZtl CL ?■ 

FRL,FDL,REORUERlf RECRGERZ $ 

LRLFRL2 $ 
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LEVEL 2.0 NASTPAN DMAP CaHP|i.ER - SOURCE LISTING 


|2« 

LAB El 

LBLTRLl $ 

209 

i 

ERL G 

CASEXX,USETU ,OLT,rRLvGnatGGD,OIT» / PPF,PSF,PDF,FOL»PHFOUn / 
C,N,DIRECT/VtN,FREQY/C,N,FREQ $ 

210 

CO'JO 

LBLFRLXlfNUFRLX S ZERO OCT LOAD COLUMNS IF FRLX HAS GENERATED. 

^211 

MPy AO 

PPFjPOZERUt / PPFX /C,N,0 £ 


tQJ IV 

PPFXf PPF $ 

H 

la 213 

LABEL 

LOLFRLXl i 

1 

CUNO 

L3LFRLl»N08ASEX S 

215 

MPVAl) 

M2GG,BASEXGt / M2BASEXG /C,NtO £ 

I 

ADD 

PPF,M2BASEXG / PPFi /C ,N , 11 . 0 ,0. 0) /C t N , ( - 1 -0 ,0 .0 ) 5 

f 

EQJ IV 

PPri.PPF 5 

»2 218 

CU'IO 

LBLBASElrNOSET $ 

f? 219 

SSG2 

USETO,GMD,YS ,KFS.GCO,,PPF / , PODUMl , PSF t , POFl $ 

y 

22 0 

tQJ IV 

PSFl.PSF // POEl.PDF 4 

1 

LABEL 

LBLBASEl $ 

222 

LAPEL 

LBLFRLl £ 

li 223 

tUJ IV 

PPF, POF/NOSE T £ 

i 

CH<PNT 

PPF,PSF,POP ,FOL £ 

^ 22 5 

PAR AML 

PDF //C.N, TRAILFR /C,N,1 /V,N,POFCGLS t. 

i 

? A^ AM 

//CfN.OIV /V.NfNLOAD / V , N , PLif C CL S /V.N.FKMAX £ NLOAD = NFZFKMAX 

22 7 

EQUIV 

PUF,PXF/CYCI C 4 

y 22 0 

CUN 0 

LBL’UUNE.C YC 10 £ 

n """ 

P ARAM 

//C,N,OIV /V,N,NLUAI) /V,N,PUFCCLS /V.Y.WSEGS £ NLDAU = NF/NSEGS 

y 

230 

CVCTl 

PDF / PXF.CCVCFl /V.N.CTYPF /C.N,FCKF. /V , Y ,NSl-GS = -1 / 
V,Y,KMAX=-1 / V,N,NLCAU /S,N,NOOO £ 

F's 

y 231 

>1 ■ 

CON U 

ERRURC i ,NOGO £ 

2.25 
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LEVEL 2.0 NASTRAN DMAP COMPUKR - SOURCE LISTING 


222 

CHKPNT 

PXF J. 

233 

JUMP 

LBLPDQNE S 

23A 

LABEL 

LBLFRL2 S ' 

235 

PAP AM 

//CtNtNUT /V,N,NGTCYCiO /V.Y.CYCIO $ 

236 

CllNO 

LBLTRL2tNCTC YCIO S 

237 

EQU IV 

POT,PDTRZl/NCRai t- 

233 

CCNO 

LBLRUlAfNORUl J 

23S 

Mpy AO 

POT.kEL'RDERl , / PD1R21 / CtN,0 $ 

240 

LABEL 

LbLRGlA i 

241 

CYCTl 

POTRZl / PXrR2l,GCVCF2 /V t N,C T YPE/C »N ,FORE/V ,N,NTSTCPS/ 
Vf Y,LMAX/VfN ,Ft<MAX/S«N,NUGO $ 

242 

ct:"< C 

ERRURC ItNOGO $ 

24 3 

CHKPNT 

PXTR21 $ 

244 

EQ'J IV 

PXTRZl.PXFZl /NCK02 1 

24 5 

CO'JO 

LBLPU2A ,NORL 2 $ 

246 

MPYAO 

PXTKZl,RE0R0eR2. / PXFZl /CfNrO S 

24 7 

LABEL 

LBLR02A $ 

248 

EUU IV 

PXFZ If PXFl 1. 

249 

CHKPNT 

PXFl $ 

250 

JUMP 

LBLTRL3 $ 

251 

LABEL 

LBLTRL2 $ 

252 

MPY AO 

PUT, REORDER 1 t / P0TPZ2 / C,N,0 $ 

253 

CYCTl 

P0TRZ2 /PXTP 22,GCYCF3 /V, N ,CT YPE/C ,N ,FCRE/V, N, NT ST EPS/ V , Y , LM AX / 
VfYfNSLGS/SfNfNCGCl » 

254 

CUV D 

ERROR ClfNuGO $ 

255 

CriAPNT 

PXTRZ2 $ 


2 . 


fO 




LEt'EL 2aO NASTRAN DMA? CCMPilEM - SOURCE LISTI.'JG 

256 EQUIV PKTRZ2,P«TR2/NCR02 £ 

257 CCMD LGLR028.N0R02 $ 

258 MPyAO PXTRZ2,REORDcR2j / PXIR2 /C,NtO £ 

259 LABEL L3LR023 3 ' 

260 CYCTl PXTR2 / PXFZ2,GCYCF9 / V, N»C TYPE/C fN ,FOF<E/V,Y, NSEGS/Vt Y.KHAX/ 

VfNTFLMAX/S,N,NCGO S 

2ol CONU ERRORClfNOGO S 

262 EUUIV PXFZ2,PXF1 i 

263 CH<PMT PXFl S 

26A LABEL LBLTRL3 $ 

265 COPY PXFl / PXP2 $ CONVERT REAL PXFl TC CCMPLEX PXf, 

266 ADD PXFi,PXF2 / PXF / C,N,(0.5,U0J / C t N , (0 . 5 ,-l . 0 ) i 

267 3A<Art //CtNfAOD /ViN,NLOAD /VfN,FLMAX /CtNiO $ NLOAD = FLMAX 

268 LABEL LBLPDONE 4 ' . ^ 

269 PARAM //CfNtADO /ViN.KINUEX /V,Y,KHIN=0 /C.NtO $ INWTIALIZE KINOcX. 

%,'• 

270 PARAM //CfNfAOO /VtN,KMINL / V t Y , KP I N VC , N f-1 1 

271 CGNO NOKM INL iKMIWL $ 

272 PARAM //CfN.ADO /V,N,KMINV /C,N,0 /C ,NV0 4 

273 JUMP KMINLOUP 5 ’ 

27A LABEL KMirJLOUP 4 

275 CYCT2 CYCUO , , ,PXF , , /,,PKFZ,, / C . N ,FORE/ V ,Y ,NS EGS/ 

V,N,KMINV/VtN,CYCSEQ/V,N,NLGAD/S,N»NCGC 4 

276 CUM C ERRUMClfNaCU 4 

277 ADO PKFZ, / CKVFZ / C ,U , (0.0 ,0.0! $ 

278 CYCT2 CYCD 0 , , , UKVF Z , , /,,UXVF,, / C , N , 8 AC K/ V , Y , NSEG5/ 

V,N,KMINV/V,N ,CYCSE>;/V,N,f>.Li:AD/S ,N,NCGC 4 

279 CUMD ERRJRClfNOGO 4 
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r 

LEVEL 2 

.0 NASTRAN DHAP COMPILER - SOURCE LISTING 

j =^^0 

PAR AM 

//C,N,AOO /V tNiKMINy /VtNjKMINV /CtN,l $ 

281 

T 

SE3T 

KHINLOOP.KKI NL ?. 

<M 

OC9 

LABEL 

NOKMINL 6 

T 283 

JUMP 

TUPCYC S 


LABEL 

TOPCYC $ LOOP ON KINDEX 

I 

COViD 

NOKPRTfNOKPR T $ 

286 

T 

P RT P ARM 

//CtN.O /CtNfKINOEX $ 

i 28 7 

LABEL 

NOKPRT $ 

I 

CyC72 

CYCOD »KOOf MOOf t » /KKRF , MKKF f , , /C ,N f PURL- / V » Y tNSE US / 
VtNf KINUEX/VfN,C VCSEa=-i/V,NfNLaAD/S,N,NCGC $ 

289 

T 

CO'IO 

EPRuRClfNOGU $ 

<1 Z^c 

CH<PNT 

KKKF,MKKF $ 

J 

P A? AM 

//C,N,SYST //C,N,53 /C,N,2 i METHCD 3T IN CYCT2 PRODUCES 

2^2 

If 

CYCT2 

CYCOO» BUD, f P XF , , /BKKFtfPKFf, / C »N rFORE/ V f Y t NSEGS/ 

VfN,KlNDEX/V,N,CYCSEQ/V,N,NLCAD/S,NfNOGa $ 

“ 2‘0 

PAR AM 

//CtUjSYST //CfN»58 /CtNfO i RESET MPYAO METHOu CONTROL. 

1 

tONU 

ERRORC l.NOGD i 

295 

CHKPM 

BKKFtPKF $ 

1 296 

f RR02 

KKKF ,BKKF ,MKKF , ,PKF ,FCL / LKVF /C t N ,0 . 0/ C t N f 0 . 0/ C . N, -1 . 0 

r» 29 7 

CHKPN T 

UKVF $ 

S 298 

CYCT2 

M.' ’ /C ,N ,BACK/V ,Y ,NS EGS/V . N, K INDEX / 

VfNf CYCSFO/V »NfNLCAC/b»N,NCGC $ 

1 299 

C.0VU 

ERRCRCl.NOGU S 

1 

CH<PNT 

UXVF i 

* 33 1 

P AR AM 

//CiNfADO / V f N jK I NDE X/ V f N f KI NOE X/C ,N , 1 J KINUEX = KINOEX 

j 3u2 

PAR AM 

//CfNjSUB /VfNjDONE / V»Y»!<FAX / VfNfKlNOEX 5 

3)3 

f 

CUMP 

LCYC2fUUNE S IF KINOEX .GT. KMAX THEN EXIT 



w.->v »*• * 

OF POOR QUAL1T\' 

lEyEL 2.0 NASTRAN DMA P CUiiPI LCR - SOURCE LISTING 


30 

REPT 

TOPCYCtiCO $ 

30 5 

JUMP 

EBP0R3 $ 

.306 

LABEL 

LCVC 2 $ 

30 7 

EQU IV 

UXVF.UUVF / CYCIO 5 

308 

CfKPMT 

UOVF $ 

309 

CUND 

LCYC3»CYCin $ IF CYCIO .GE. 0 THEN TRANSFORM TO PHYSICAL. 

310 

CYCTl 

* V,N,CTyPE/C,N,BACfVV,Y,NSEGS/V.Y,KMAX/ 

V»NfNLOAC $ 

311 

CHKPNT 

UOVF $ 

312 

LABEL 

LCYC3 S 

■»13 

CUVO 

LBLTRLA.NOTl PE $ 

31^ 

6UJ IV 

PXF,PDF2 / CYCin % 

315 

CUV D 

LCYC4fCYCI0 $ IF CYCIQ .GE. 0 THEN TRANSFORM TO PHYSICAL. 

316 

CYCTl 

PXF / P0F2,GCYC»2 / V ,N ,C T YPE /C ,N , 0ACK/V , Y ,NSEGS/V ,Y , KMAX/ 
VfNfNLUAD S 


31 7 

LAB EL 

LCYC4 $ 

318 

s Oi^ 1 

USETOt tPDF 2f j fOnO iGMOt t » t / PPFZ,t /CfN,l / C »N f DY NAM I CS $ 

319 

SSG2 

USETO fGHDt YS f KF S »GOu -tPPFZ / » PCDUP , PSF Z t PLOUM $ 

320 

fcOJ IV 

PPFZ ,PPF // PSFZ.PSI- 4 

32 1 

UKPNT 

PPFjPSF S 

322 

LABEL 

LBLTRL4 $ 

32 3 

V OR 

CASEXX.EQD YN fUSL'TD »CUVF »F CL » X YCOB ,/ QUDVC 1 i/C»N,FREUPLS P/Ci N. 
DlPGCT/SfN.N0S0PT2/S,N, NOD/S, N,NCP/C,N,0 4 

324 

LUND 

LBL15,N00 $ 

32 5 

COND 

LBL 1 5A ,NUSGK T2 i 

326 

SDR 3 

0UUVCl,,,,,/CUUVC2,,,,, $ 

32 7 

0 FP 

0U0VC2, , ,, ,//S,NrCARONC i 
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, ■QUALITY 

t-EVEL 2.0 NASTkAN DMAP COMPtLEk -* SOURCE LISTING 




I 


328 

KYTPAN 

XYC0BfGU0VC2f ,, 5 /Xy(>LTFA/<‘FREQ*/«=OSET«/S »N »PF I LE/S t N» C ARDNG S 

329 

X VP LOT 

XYPL TFA// $ 

330 

JUMP 

LBL15 $ 

331 

LABEL 

LBL15A 5 

332 

QFP 

OUDVCl, , ,r »//SrN,CARl)NC S 

333 

lABEL 

LBL15 A 

334 

COND 

LBL20,NUP 5 

333 

LOU IV 

UDVF.UPVC/NOA $ 

336 

COND 

LBL19,NUA $ 

33 7 

S OR 1 

USETD» fCDVFt *fGllDtCMO,PSr,KFSf ,/UPVCttQPC/l/«‘DYNAMICS* $ 

338 

LABEL 

LBLL9 i 

339 

.SD-12 

CASEXX tC STM, MPTf 0 1 T ,EQDYN , SI LO » t f BGPOP »F UL »OPC ,tJ PVC, EST , XV COB, 
PPF/OPPC l.QuPCUOUPVCl ,CESC1 ,CEFC 1 ,PUPVCi/C, N, FRcURFSP/ 
S,N,N0SURT2 $ 

3AO 

CL NO 

LBL17,NUSURT2 S ' 

3'tl 

SDR 3 

OPPC 1 , UOPC 1 , CUPVC i .UESCl , OEFC 1 , /0PPC2 ,0QPC2 , CJ PVC2 , OFSCP, 

□ EFC 2, 4 ■ ' , 

342 

OF? 

0PPC2,UQPC2,CUPVC2,CEFC2,CFSC2 ,//S,N,CARnNC f 

343 

XYTk AN 

XYCl)B,UPPC2,CQPC2,tiCPVC2,GESC2,CEFC2/<VPL7F/<=FREJ*/*PSEl V S. 

N,PFILE/S,N,CAFnNQ i .jui / 5 , 

344 

XYPLOT 

XYPL TF// $ - 

34 5 

CUNO 

LBL16,NQPS0L S 

346 

FAN OL.M 

XYCOBfO lT,PS0L,tJUPVC2fCPPC2,CCPC2, LESC2 ,0EFC2 ,CAS6XX/PS0F, AUTU/ 
S,N,NCRD i 

34 7 

CUNO 

LBL16,N0R0 1 

34 8 

XYTRAN 

XYCUB ,PSLf-,AUTU,, ,/XYPL7R/*P.AND<=/*P5FT VS ,N,PFILE/ S.N, 

CAR UNO i 

349 

<YPLOr 

XYPLTR// 5 
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OF FO03 QDALm^ 


LEVEL 2.0 NASTRAh! OMAP COMPILER - SCURCE LISTINC- 


350 


i 


351 

352 
J53 
35A 


35 5 


356 
35 7 
35 3 


fcl 


359 


350 


1 

362 





36 A 
3o5 
366 
36 7 
366 
3o9 
370 


I 

372 
I 373 
37 ^ 


* JU'IP L8L16 S 

LAO EL LBL17 $ 

'JF.» UUPVClfOPPClvOQPCl,CEFCl»CE5Cl ,//S,KtCAr<!)AC $ 

LABEL LBL16 * 

COviD LBL20, JUMPPLQT i 

PLOT PLTPAR ,GPS£TS,ELSETS,CASEXX,BOPDT,EaEXIW,SIP,rPJPVClf 
OESC 1/PL0TX2/NSI L/LUStP/JLMPPLGT/PLTFLG/ StNfPFILE $ 

PRTMSG PL0TX2// t 

LABEL LBL20 t 

CUVU FINI S-.REPtATF I. 

REPT LBL13, 100 & 

LABEL ERROR 3 $ 

PRTPARM //-3/*0IRFRKD<‘ t 

JUMP FINIS $ 

LABEL ERRGRA $ 

PRTPARM //-A/«DIP.FRRO« li 

LABEL ERRURCl $ CHECK NSEGS, KMAX AND OTHER CYCLIC DATA. 

PRTPARM //C,N,-7 /CiK.CYCSTATICS i 

LABEL EKRJRC2 $ COUPLED MASS NOT ALLChED. 

PRTPARM //CfNfO /C • Y tCClJP PA SS S 
JUMP FINIS i 

label ERRURC3 i SUPCPxT BULK DATA NCT ALLOWED, 

PRTPARM //CfN,-6 /C r N ,C YC ST A TI C S 1- 

LABEL ERRORCA S EPCjINT BULK DATA NuT ALLCWED. 

PRTPARM //CfN,0 /C,N,NOUE i ' 


GPECT* 


JUMP 


F INI S $ 



JLuf<;4 



PAGH 53 

POOR QUALsT/ 


LcVEL 2.0 fMSTRAN DMAP COMPILER - SCcRCE LISTING 


! 375 

LABEL 

ERR0RC5 

5 neither FREQ OR TSTEP HERE IN BULK DATA DECK. 

376 

PRTP ARM 

//C.NrO 

/C,NtNCFRL S ■ 

; 377 

PRTPARM 

//C,N,0 

/CtN.NCTRL $ 

• 378 

JUMP 

FINIS $ 


379 

LABEL 

ERRDRC6 

S BOTH FREC AND TSTEP WERE SELECTED IN CASE CONTROL. 

38 0 

PRTPARM 

//C , N t 0 

/CtUtNCFREQ $ 

381 

’RTPARM 

//C,N»0 

/C»N,NUTIME 5 

. 38 2 

JUMP 

FINI S i 


333 

LABEL 

FINIS $ 


38A 

END 
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FORCED VIBRfT;;^^: or n:-TATING CYCLIC STRIICTURES 

2.4.3 Description of DHAF Ocgrationj for Forced Vibration Analysis of Rotating 
Cyclic Structures 

5. Go to DMAP No. 365 if user has not specified parameter MSEGS. 

6. Go to DMAP No. 365 if user has not specified parameter KfiAX. 

8. Go to DMAP No. 365 if user has not specified parameter CYCIO. 

11. Go to DMAP No. 365 if KMAX > NSEGS/2. ' . 

18. Go to DMAP No. 367 if user has requested consistent mass. . 

20. GPl generates coordinate system transformation matrices, tables of grid point 
locations, and tables to relate internal to external grid point numbers. 

23. Go to DMAP No. 141 if only Direct Matrix Input. 

24. GP2 generates Element Connection Table with internal indices. 

27. Go to DMAP No. 35 if no plot output is requested. 

28. PLTSET transforms user input into a form used to drive structure plotter. 

29. PRTHSG prints error messages associated with structure plotter. 

32. Go to DMAP No. 35 if no undeformed structure plots are requested. 

33. PLi3T generates all requested undeformed structure plots. 








II 


i 


34. PRTMSG prints plotter data and engineering data for each undeformed plot 
generated. 

36. 

GP3 generates Grid Point Temperature Table. 

38. TAl generates element tables for use in matrix assembly and stress recovery, 

41. GP4 generates flags defining members of various displacement sets (USET) and 

forms multipoint constraint equations [Rg] {u^^} = 0. 

45. Go to DMAP No. 370 and print error message if free-body supports are present. 

46. DPD generates flags defining members of various displacement sets used in dynamic 
analysis (USETD), tables relating internal and external grid point numbers, 
including extra points introduced for dynamic analysis, and prepares Transfer 
Function Pool, Dynamics Load Table, Poiver Spectral Density List and Frequency 
Response List. 

48. Go to DMAP No. 375 and print parameters NOFRL and NOTRL if there was no FREQ 
or TSTEP bulk data. 

55. • Go to DMAP No. 379 and print parameters NOFREQ and NQTIME if both FREQUENCY and 
TSTEP were requested in the Case Control deck. 
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57. Go to DMAP No. 372 and print paranieter NOUE if extra points are present. 

58. GPCYC prepares segment bcimdary table (CYCDD). 

59. Go to DiiAP No. 365 and print error messoge if CYJOIN data is inconsistent.’ 

61. Go to DiMAP No. 95 if there are no structural elements. 

66. EMQ generates structural element stiffness, mass, and damping matrix tables 
and dictionaries for later assembly. 

67. Go to DWP No. 69 if no stiffness matrix is to be assembled. 

68. EMA assembles stiffness matrix [•‘^gg3 Grid Point Singularity Table. 

73. Go to DliAP No. 75 if no mass matrix is to be assembled. 

74. EHA assembles mass matrix 

•- ggJ 

76. DUMMODl generates modified Frequency Response List, FRLX, Coriolis acceleration 

coefficient matrix [BIGG ], centripetal coefficient matrix [HIGG Base 

yU 99 rr 

acceleration coefficient matrix [K2GGggJ, Base acceleration matrix [BASEXG^^] and 

load modification matrix, [PDZEROg], for base acceleration problems. 

79. Equivalence FRLX to FRL if FRLX was generated by DUKMODl . 

82. Go to DMAP No. 84 if no viscous damping matrix is to be assembled. 

83. EMA assembles viscous damping matrix [BggJ. 

85. Go to DMAP No. 87 if no structural damping matrix is to be assembled. 

86. EMA assembles structural damping matrix C^'ggJ* 

91. Go to DMAP No. 95 if no weight and balance is requested. 

92. Go to DMAP No. 363 and print error message if no mass matrix exists. 

93. GPWG generates weight and balance infonnation. 

94. 0FP formats weight and balance information prepared by GPWG and places it on the 
system output file for printing. 

96. Equivalence [K J to [K J if no general elements. 

97. Go to DMAP No. 99 if no general elements. 

98. SMA3 adds general elements to [Kg^] to obtain stiffness matrix [K^g]. 

100. Go to DMAP No, 109 if parameter RPS = 0.0 orif no mass matrix is present. 

104. ADD assembles the Coriolis acceleration matrix into the viscous damping matrix 

= [Bgg] ♦(■!-. RPS) [BlGGgg] 
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105. 

106. 

107. 

no. 

111 . 

114. 

115. 

116. 

118. 

119. 

120 . 

121 . 


123. 


I' 


' ' ■ GRIGiKAL hAGi 

- OF POOR QUALfT!' 

Equivalence [BGGl „] to 

99 9a 

ADD assembles the centripetal acceleration matrix into the stiffnesv matrix. 
[KGGlgg] = [KggO ~ ( 2 7T . R R ^ ^ [MlQGgg] 

Equivalence [KGGlgg] to 

Go to OMAP No. 117 if general elements present. 

Go to D.W No. 117 if no structural elements. 

GPSP determines if possible grid point singularities remain. 

Go to DHAP No. 117 if no grid point singularities exist. 

0FP formats the tableof possible grid point singularities prepared by GPSP and 
places it on the system output file for printing. 

Equivalence [Kg^] to [M^g] to [Bgg] to [B^J and [K^gl to [kJj 

if no multipoint constraints. 

Go to DMAP No. 122 if MCEl and NCE2 have already been executed for current set 
of multipoint constraints. 

fiCEl partitions multipoint constraint equations [RgJ = j R„3 and solves for 
multipoint constraint transformation matrix [G^] = 

MCE2 partitions stiffness, mass and damping matrices 



'^nn 

* ^'nml 


r 

nn 

'■U 

V = 

— _ 
^mn 1 *^mrn 


7 

Jmn 



B 

1 B 



K^l 


nn ^ nm 

8 1 B 

_mn 1 mm. 

• ' 

nn 

_mn 

nm 

V" 

mm_ 


and performs matrix reductions 


[K„„] - 

[H„g ■ 

[B 1 = [B 3 + [G^][B 3 + [B^ 3[6 3 + [G^][B 3[G ], 

^ nn-' nrr '• m-’'- mn-^ '■ mn-”- m-* m-^*- mm-"- m-*’ 

[K^ 3 - [K^ 3 + [rJl[K^ ] + 3 ’’’[g 1 + [g'''3[k^ 3[g ]. 

^ nn-* '■ nn-* m'*- mn-* mn-* *- m' *■ m-''- mm-”- m-* 


Equivalence [Kj^^3 to [K^^3 . to to [B^^3 and [K^^3 to [kJ^3 

if no singlepoint constraints. 


2.35 



124. 

125. 


127. 

128. 

129. 

130. 

131. 

132. 


133. 

134. 


136. 


Go to DHAP No. 126 if no single-point constraints. 
SCEl partitions out single-point constraints 


POOR QUALITY 





”ff 1 "fs 

nn-* 

- + - 
b‘ 1 =ii 

1 ^fs 

a' 

1 4s 

.4 T- 4 
^'sf 1 hs 


y|“i 

and [k; 1 = 
'• nn** 


Equivalence to if no omitted coordinates. 
Equivalence to if no omitted coordinates. 
Equivalence to if "o omitted coordinates. 
Equivalence to if no omitted coordinates. 
Go to DMAP No. 141 if no omitted coordinates. 

SMPl partitions constrained stiffness matrix 


[Kf,] = 

solves for transformation matrix [G ] = 
and performs matrix reduction , 

[>^13 ■ 

Go to DMAP No. 135 if n' mass matrix. 

SMP2 partitions constrained mass matrix 


^aa I f^ao 
— .[_ — 

l^oa 


-C'=o<.-T'P<oa,] 


H 

1 

H 

aa 


ao 

— 


— 

M 


M 

oa 

1 

00 


and performs matrix reduction 


Go to DMAP No. 138 if no viscous damping matrix. 
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137. 


SMP2 partitions ccnSts'oin'fid viscoits damping matrix 

■f* 


GRiCSrAL 'i- 

OF POOft Q'JAl-iTY 


[B,,J 


^aa < »ao 


oa 


B. 


00 


139. 

140 


and performs reduction 

[bL] ■ * Cc,„][G„j . [b^„g„]T . [gJ][b„„][g„] 

Go to DiiAP No. 141 if no structural damping matrix. 

SHP2 partitions constrained structural damping matrix 


c4f] = 


4. I 

aa . ao 


_ J_ 

111 I '^ooj 


l.4’ 


142. 

145. 

148. 

14; 

151. 

152. 

154. 

155. 

158. 

159. 


and performs matrix reduction 

.. [Gj]rK^„KG„] 

Equivalence [Gq] to [G^] and [G^^jJ to [G^ if no extra points introduced for 
dynamic analysis. 

BHG generates BMIG card images describing the interconnection of the fluid 
and the structure. 

Go to DMAP No. 150 if no fluid structure interface is defined. 

MTRXIH generates fluid boundary matrices [Aj^ and [K^ if a fluid 
structure interface is defined. The matrix [Kj^ is generated only for a 
nonzero gravity in the fluid. 

Go to next DMAP instruction if cold start or modified restart. LBL13 will 
be altered by the Executive System to the proper location inside the loop for 
unmodified starts within the loop. 

Beginning of loop for additional sets of direct input matrices. 

CASE extracts user requests from CASECC for current loop. 

HTRXIN selects the direct input matrices for the current loop, 

. ro2 T pp-J* >■ ppJ 

and [Bppl. 

Equivalence to [H-p] if no [Aj^^^j.]. 

ADDS adds end [Kpp,] and subtracts [A^ from, them to form 
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160. Go to DflAP No. 163 if r - 


OivlQJK.AL 

OF POOR QUALITY 


161. Transpose ^^1 to obtain 

162. ADD assembles input matrix D’pp] ~ HFACT [Aj^ + CMppj. 


169. Go to DMAP No. 172 if transient type GKAD matrices are to be generated. 

170. Equivalence [MppJ to [BppJ to [8^^] and [KppJ to if no constraints 

applied, to it no direct input mass matrices and no extra points 

and [Bgg.] to if no direct input damping matrices and no extra points. 

172. Go to DMAP No. 175. 

173. Equivalence [Mpp] to [H^^l, [BppJ to [B^^l and [Kpp] to if no constraints 

applied, [M^gJ to it no direct input mass matrices and no extra points, 

and [Kgg.] to iF no direct input stiffness matrices and no extra points. 

176. Go to DMAP No. 178 if only extra points are defined. 

177. GKAD assembles stiffness, mass, and damping matrices for use in Direct Frequency 
Response, if parameter GKAD = FREQRESP. 

[K„j] - (1 t . 

CBj,] = pJdi + t“dd^ • ^ ; / 

Direct input matrices may bo complex. 

[13 

GKAD asserribles stiffness, mass, and damping matrices for use in Direct Transient 
Response if parameter GKAD = TRANRESP. 

= [kJj] * , 

tMdd] = t [B^jl . 

and [B,,j . [bJj] * [6^,] * A [kJ,] + i , 


where 


*^aa j 

0 

o 

0 

U. 1 

•J 

r" 1 

"1 


0 

aa 1 


0 ' 0 


I 


^ n4dJ • 


2.38 



and I 0 


j n ^ • j 
‘•‘'dir 




All matrices are real. 

Go to DMAP No. 1Q2 if transient type GKAD matrices were generated. 

Equivalence to if all stiffness is Direct Matrix Input, [M^j] to 

mass is Direct Matrix Input and [B^^] to if all damping 

is Direct Matrix Input. 

Go to DMAP No. 134. 

Equivalence [B^^] to [B^jjl if all damping is Direct Matrix Input, [M^^J to 
[Mdd] ii" all i^ass is Direct Matrix Input and to is all stiffness 

is Direct Matrix Input. 

Go to DMAP No. 203 if loading is frequency-dependent. 

Beginning of loop for additional subcases for time-dependent loads. 

CASE extracts user requests from CASECC for th current loop. 

TRIG generates matrices of loads versus time. {Pl^l is generated with one 
column per output time step. tPl^} is generated with one column per solution 
time step, and the Transient Output List (TOL) is a list of output time steps. 

SDRl appends {Pljl to {P^}. 

SDRl appends (Plj) to (Pj). 

Go to DMAP No. 200 if no additional time-dependent loads need to be processed. 

Go to DMAP No. 189 if additional time-dependent loads need to be processed. 

Go to DMAP No. 360 and print message if more then 100 loops. 

Equivalence {P^} to {P^} if the output times are the same as the solution times. 

DUMM0D2 generates a Frequency Response List (FRLZ), and a Frequency Output List 
(FOL), from the Transient Output List (TOL). Load reordering matrices REORDERl 
and RE0RDER2 are generated based on parameter values. This module, in effect, 
generates data blocks necessary to convert time-dependent loads into frequency 
dependant loads. 


Equivalence FRLZ to FRL and FOLZ to FOL. 
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207. Go to DMAP No. 234. v, . 

209. FRLG forms the dynamic load vectors (Pp) , (P^) and Frequency Output List 

(FOL) for frequency-dependent ‘loads. 

210. Go to DMAP No. 213 if FRLX was not generated by DUMMODl . 

211. HPYAD uses PDZERO from DUMMODl to zero out selected columns of {P‘} in base 

P 

acceleration problems. 

212. Equivalence {PXp} to {Pp}. 

214. Go to DMAP No. 222 if not a base acceleration problem. 

215. MPYAD forms the complete base acceleration matrix, {MZBASEXG^} = [M2GG ] • 

f y 99 

{BASEXGM. 
y 
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216. 


217. 

218. 
219. 
228. 

230. 

231. 
233. 

236. 

237. 
233. 
239. 

241. 

242. 
244. 
246. 
248. 
250 

252. 

253. 

254. 

256. 

257. 

258. 
260. 
261. 
262. 
265. 


ADD assembles the freq'^.,,^ 




aUus and the loads due to base acceleration 


{Pip} = {Pp> {M23A3F.X8J} ' , 

Note that the p-set and g-set are the .sa.ma because no extra points are allowed. 
Equivalence {Pip} to {Pp). ... t 


Go to DMAP No. 221 if there are no .SPC’s, MPC's or OMITS. 

SSG2 applies constraints to {Pp}. 

Go to DMAP No., 268 if parameter CYC 10 = -1. 

CYCTl transforms loads on analyses points to symmetric components. 

Go to DMAP No. 365 and print error message if a CYCTl error was found. 


Go to DHAP No. 268. 

Go to DM.AP No. 251 if parameter CYCIO = +1. 

Equivalence {P^} and {PDTR7.1 } if REORDERl was not generated by DUMM0D2. 

Go to DMAP No. 240 if REORDERl was not generated. 

HPYAD reorders columns of {P^}. 

CYCTl transforms loads on analysis points to symmetric components, in time. 
Go to DMAP No. 365 and print error message if a CYCTl error was found. 
Equivalence {PXTRZl} and {PXFZl} if RE0RDER2 was not generated by DUMM0D2. 
MPYAD reorders columns of {PXTRZl}. 

Equivalence {PXFZl} to {PXFl}. 

Go to D.MAP No. 264. 

MPYAD reorders columns of {P^l. 

CYCTl transforms loads on analysis points to summetric components, in time. 
Go to DMAP No. 365 and print error message if a CYCTl error was found. 
Equivalence {PXTRZ2} to {PXTR2} if RE0RDER2 was not generated. 

Go to DMAP No. 259 if RE0RDER2 was not generated. 

MPYAD reorders columns of [PXTRZ2}. 

CYCTl transform symmetric components, in time, to sumnetric components. 

Go to DMAP No. 365 and print error message if a CYCTl error was found. 
Equivalence (PXFZ2} to {PXFl}. 

COPY makes a physical copy of iPXFl; called {PFX2}. 
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265. ADD makes loads complex,, since SDR2 fci/.p'icts complex loads in a frequency response 
problem. Time-dependent loads are real > 

{PFX} == (0.5, i.Oj) * O'FXl) (0.5, -1.0^) • {PFX2} 

271. Go to DHAP No. 282 if 0. , ^ 

27'4. Beginning of loop to create It'lIN null columns of {UV^} for KIMDEX = 0 to 

(KfiIN-1). Tiiese leading null columns are necessary because CYCTl expects 
columns for KINDEX = 0‘tc K.M.Aa. 

275. CYCT2 transforms loads from symmetric components to solution set for rota- 
tional symmetry. This operation is necessary to get a correct size matrix for 
generating null {UVf} columns. 

276. Go to DMAP No. 365 and print error massage if a CYCT2 error was found. 

277. ADD generates a null vector {UVZ^} = {PZ^}* 0.0. 

278. CYCT2 finds symnetric components of displacements from solution set data and 
appends it to {UVj[}. 

279. Go to DMAP No. 365 and print error message if a CYCT2 error was found. 

280. PARAM increments the value of KMINV=KMINV+1 . 

281. Go to DMAP No. 274 if more null vectors are to be generated for {UV^}. If 
the initial {UV^} for KINDEX values 0 to (KMIN-1) has been completed then 
go to DMAP No. 282. 

284. Beginning of loop for cyclic index value (KINDEX), for values KINDEX = KHIM 
to KMAX. 

283. CYCT2 transforms stiffness and mass matrix from symmetric components to 
solution set for rotational symmetry by the equation: 

291. CYCT2 transforms damping and loads from symmetric components to solution 
set for rotational symTietry by the equations: 

[PfcJ = [gJ] tPy t [gJ] (P^) 

294. Go to DMAP No. 365 and print error message if CYCT2 error was found. 

296. FRRD2 solves for the displacements using the following equation: 

+ '®dd" * ' 'V- 

298. CYCT2 finds symmetric components of displacement from solution set data and 

appends to output for each KINDEX. . . 

299, Go to DMAP No.. 365 and print, error message if CYCT2 error was found. 



301. PARAM increnents the value of KlfJDEX = KINDEX + 1. 

303. Go to DMAP No. 306 if all cyclic index values are complete. 

304. Go to DMAP No. 284 if cdditional index values are needed. 

305. Go to DMAP No. 360 and print error message if more than 100 loops on KINDEX. 

307. Equivalence {U^} to {11^} if parameter CYCIO = -1 . 

309. Go to DMAP No. 312 if parameter CYCIO = -1. 

310. CYCTl transforms displacements from symmocrical components to physical components 

313. Go to DMAP No. 322 if loads wart- frequency-dependent. 

314. Equivalence {P^} to {P2^} if parameter CYCIO = -1 . 

315. Go to DMAP No. 317 if parameter CYCIO = -1. 

316. CYCTl transforms loads from symm.et 'ical components to physical car.ponents if 
loads v/ere time-dependent. 

f 

318. SDRl recovers dependent loads {PZp}. 


319. SSG2 applies constraints to {PZ ^1 to form {PZ^}. 

P ■ S ' 

320. Equivalence {PZ p} to {Pp} and {PzJ} to {Pg}. 

323. VDR prepares displacements, sorted by frequency, for output using only the 
independent degrees of freedom. 

324. Go to DMAP No. 333 if no output request for the independent degrees of freedom. 

325. Go to DMAP No. 331 if no output request for independent displacements sorted 
by point number. 

326. SDR3 sorts the independent displacements by point number. 

327. 0FP formats the requested independent displacements, sorted by point number, 
prepared by SDR3 and places them on the system output file for printing. 

328. XYTRAM prepares the input -^or X-Y plotting of the independent displacements 
vs. frequency. 

329. XYPL0T prepares the requested X-Y plots of the independent displacements 
vs. frequency. 

332. 0FP fonnats the requested independent displacements , sorted by. frequency , 
prepared by VDR nd places them on the system output file for printing. 
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■334. Go to DHAP No.. 357 if no dependent degrees of freedom 

for forces ind stresses. 

335. Equivalence {u^} to {Up} if no constraints applied. 

336. Go to DMAP No. 338 if no constraints applied. 

337. SDRl recovers independent cor;po.o£rits of displacements 



and recovers single-point forces of constraing {q^} = -(Pg) + [K^ 2 ]{uj:} . 

339. SDR2 calculates element forces (0EFC1) and stresses (0ESC1) and prepares load 
vectors (0PPC1), displacement vectors (0UPVC1), and single-point forces of 
constraint (0QPC1) for output sorted by frequency. 

340. Go to DMAP No. 351 if no output requests sorted by point number of element 
number. 

341. SDR3 prepares req'iested output sorted by point number or element number. 

342. 0FP formats tables prepared by SDR3, sorted by point number or element number, 
and places them on the system output file for printing. 

343. XYTRAN prepares the input for requested X-Y plots. 

344. XYPL0T prepares the requested X-Y plots of displacements, forces, stresses, 

loads or single-point forces of constraint vs. frequency. . 

345. Go to DMAP No. 353 if no Power Spectral Density List 

346. RANDOM calculates power spectral density functions (PSDF) and autocorrelation 
functions (AUT0) using the previously calculated frequency response. 

347. Go to DMAP No. 353 if no RAND0M calculations requested. 

343. XYTRAN prepares the input for requested X-Y plots of the RAND0M output. 

349. XYPL0T prepares the requested X-Y plots of autocorrelation functions and 

power spectra! density functions. 

350. Go to DMAP No. 353 if no frequency response output requests sorted by frequency 
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352. 0FP formats frequer-y /osp'aso requests prepared by. SDR2, sorted fay. 

frequency, and places them on the syst-H output file for printing. 

354. Go to DMAP No. 357- if no deformed structure plotj are requested. 

355. plot generates all requested d&fomvad plots. 

356. PRTMSG prints plotter data and engineering data for each deformed plot generated. 

358. Go to DMAP No. 383 if no additional sets of direct input matrices need to be - 
processed. 

359. Go to DMAP No. 152 if additional sets of direct input matrices need to be 
processed. 

361. DIRECT FREQUENCY AND RAMD0H RESP0NSE ERR0R MESSAGE NO. 3 - ATTEMPT T0 EXECUTE 
M0RE THAN 100 L00PS. 

364. DIRECT FREQUENCY AND RANDi3M RESPONSE ERR0R MESSAGE N0. 4 - MASS MATRIX REQUIRED 
F0R WEIGHT AND BALANCE CALCULATIONS. 

366. STATICS WITH CYCLIC SYMMETRY ERR0R MESSAGE N0. 7 - CYCLIC SYMMETRY DATA ERROR. 

368. Coupled mass is not allowed - Print parameter COUPMASS. 

371. STATICS WITH CYCLIC SYMMETRY ERR0R MESSAGE N0. 6 - FREE-B0DY SUPPORTS NOT 
ALLOWED. 

373. EPoint bulk data not allowed - Print parameter NOUE. 

376. Neither FREQ or TSTEP were in bulk data - Print parameters NOFRL and NOTRL. 

300. Both FREQ and TSTEP were selected in case control - Print parameters NOFREQ 
and NOTIME. 

384. END of DMAP sequence. 

1 

I 
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2 . 4.4 CASE CONTROL DECK Arllj ANALYSIS OF 

ROTATING CYCLIC STRUCTURES 


The following items relate to subcase definition and data selection for 

Forced Vibration and Random Response of Rotating Cyclic Structures; 

1.. The SPC and MPC request must appear above the subcase level and may hot be 
changed. 

2. Either FREQUENCY or TSTEP must be selected and must be above the subcase level. 

3. If selected, FREQUENCY must be used to select one and only one FREQ, FREQl or 
FREQ2 card from the Sulk Data deck. 

4. If selected, TSTEP must be used to select the time-steps to be used for load 
definition via a TSTEP Sulk Data card and must be defined above the subcase 
level. 

5. Direct input matrices are not allcwed. 

6. OFREQ must not be used. 

7. A separate group of subcases must be defined for each symmetric segment. 

8. DLOAD must be used to define a frequency or time-dependent loading condition 
for each subcase. For frequency-dependent loads, subcases without loads need 
not refer to a DLOAD card. For time-dependent loads, subcases without loads 
must refer to a DLOAD card that explicitly generates a null load. 

9. An alternate loading method is to define a separate group of subcases for 

each harmonic index, k. The parameter CYCI0 is included and the load components 
for each index are defined directly within each group for the various loading 
conditions. 

10. If Random Response calculations are desired, RANDOM must be used to select 
RANDPS and PANDTi cards from the Bulk Data Deck. 

The following printed output, sorted by freqeuncy (S0RT1) or by point number 
or element number. (S0RT2) , is available, either as real and imaginary parts or 
magnitude and phase angle (0° - 360° lead), for the list of fequencies specified: 

1. Displacements, velocities, and accelerations for a list of PHYSICAL points 
(grid points and extra scalar points introducca for dynamic analysis) or 
SOLUTION points (points used in formulation of the general K system). 

2. Nonzero components of the applied load vector and single-point forces of 
constraint for a list of PHYSICAL points. 

3. Stresses and forces in selected elements (ALL available only for SOPJl). 


The following plotter output is available for Frequency Response 
calculations: 
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1 . 


lindeformed plot of ths strUctuval Kiodel. , 

2. X-Y plot of any component of displacenvent. velocity, or acceleration of a 
PHYSICAL point or SCLUTIgN point. 

3. X-Y plot of any component of the applied load vector or single-point force 
of constraint. 

4. X-Y plot of any stress or fores component for an element. 


. The following plotter output is available for Random Response 
calculations: 

1. X-Y plot of the power spectral density versus frequency for the response of 
selected components for points or elements. 

2. X-Y plot of the autocorrelation versus time lag for the response of selected 
components for points or elements. 

The data used for preparing X-Y plots may be punched or printed in tabular form 
(see Section 4.3). This is the only form of printed output that is available 
for Random Response. Also, a printed summary is prepared for eacli X-Y plot which 
includes the maximum and rnininium values . of the plotted function. 


The following items relate to Bulk Data restrictions: 

1. SUPORT cards are not allowed. 

2. EPOINT cards are not allowed. i i 

f 

3. SPOINT cards are not allowed. 

t , 

4. CYJOIN cards are required. 

5. If a TSTEP card is used then it must not be continued since only one uniform 
time step interval must be specified. The skip factor for output, NO, on 
the TSTEP card must be 1. 


The following parameters are used in Forced Vibration and Random Response 

of Rotating Cyclic Structures: 

1. GRDPNT - optional - A positive integer value of this parameter will cause the 
Grid Point Weight Generator to be executed and the resulting weight and 
balance information to be printed. All fluid related masses are ignored. 

2. WTMASS - optional - The terms of the structural mass matrix are multiplied 
by the real value of this parameter when they are generated in EMA. Not 
recommended for use in hydroelastic problems. 

3. COUPMASS - fixed - Only lumped mass matrices must be used. 
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4. GKAD - optional - The SCD valne of tfivs paramster is used to tell the GKAD 
module the desired form of matrices KDD.BDD and MOD. The BCD value can be 
FREQRESP or TRANRESP. The default is TRANRESP. 

NOTE: Remember to define parameters G, W3 and W4. See Section 9.3.3 (DIRECT 

DYNAMIC MATRIX ASSEMBLY) .Pages 5.3-7 and 9.3-8 of the NASTRAN theoretical 
manual for further details . . 

5. LGKAD - optional - The integer value of this parameter is used in conjunction 
with parameter GICAD. If GKAD = FREQRESP then set LGKAD = 1, if GKAD = TRANRESP 
set LGKAD = -1. The default value is -1. 

6. G - optional - The real value of this parameter is used as a uniform structural 
damping coefficient in the direct formulation of dynamics problems. Not recom- 
mended foruse in hydroslastic problems (use GE on MATl). 

7. W - optional - The real value of this parameter is used as a pivotal frequency 

for unifonn structural damping if parameter GKAD = TRANRESP. In this case 

W3 is required if uniform structural damping is desired. The default value 
is 0.0. 

8. W4 - optional - The real value of this parameter is used as a pivotal frequency 

for element structural damping if parameter GKAD = TRANRESP. In this case W4 

is required if structural damping is desired for any of the structural 
elements. The default value is 0.0. 

9. NSEGS - required - The integer value of this parameter is the number of 
identical segments in the structural model. 

10. CYCI0 - required - The integer value of this parameter specifies the fonn of 
the input and output data. A value of +1 is used to specify physical segment 
representation, and a value of -1 for cyclic transform representation. There 
is no default. 

11. CYCSEQ - fixed - The integer value of this parameter specifies the procedure 
for sequencing the equations in the solution set. A value of +1 specifies 
that all cosine terms should be sequenced before all sine terms, and a value 
of -1 for alternating the cosine and sine terms. The value of CYCSEQ has 
been set to -1 . 

t 

12. CTYPE - fixed - The BCD value of this parameter defines the type of cyclic 
summetry as follows: 

(1) R0T - rotational symmetry 

13. Kf'lAX - required - The integer value of this parameter specifies the maximum 
value of the hatinonic index. There is no default for this parameter. The 
maximum value that can be specified is MSEGS/2. 

14. Kf'ilN - optional - The integer value of this parameter specifies the minimum 
value of the harmonic index to be used in the solution loop. KMIN can 
equal KflAX. The default is 0. 

15. NL0AD - fixed - The integer value of this parameter is the number of static 
loading conditions. The value of NLOAD is internally computed. 
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16. NOKPR T - optional - An integer value of -)1 for this parameter will cause the 
current harmonic index, KINDEX, to be printed at the top of the harmonic loop. 
The default is +1 . 

17. LMAX - optional - The integer value of this parameter specifies the maximum 
harmonic in the fourier decomposition of periodic, time-dependent loads. 

The default value is NTSTEPS/2, where NTSTEPS = N+2 where N is from the 
TSTEP bulk data card. 

18. RPS - optional - The real value of this parameter defines the rotational 
speed of the structure in revolutions per unit time. The default is 0.0. 

19. BXTID, BYTID, BZTID, BXPTID, BYPTID, BZPTID - optional - The positive integer 
values of these parameters define the set identification numbers of the TABLED 
bulk data cards which define the components of the base acceleration vector. 
The tables referred to by BXTID, BYTID and BZTID define magnitude (LT-2) and 
tables referred to by BXPTID, BYPTID and BZPTID define phase (degrees). The 
default values are -1, which means that the respective terms will be ignored. 
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3.1 ISATA BLOCK AND TV^flLE DESCRIPTION 

3.1.1 Data Blocks Outp uS fr om L ' Ciiule O'JMMODl 

3. 1.1.1 FRLX (TABLE) 

Description 

Frequency Response List 

The FRLX contains one log.,... record for each different set defined in the 
bulk data. Each record contains a sorted list of frequencies defined in the set 


Table Format 

Record 

Word 

Tvee 

Item 

0 

1, 2 

BCD 

Data block name 


3 

I 

Set ID^ 


2+n 

I 

SetID„ 

1 

1-m 

R 

Radian frequencies belonging to set ID, 
(w = 27 iF) ' 

* 



n 

1-k 

R 

Radian frequencies belonging to set ID„ 

nfl 



(w = 2 nF) " 

End-of-file 


Table Trailer 


Word 1 = number of rrequency sets 
Word 2-6 = zero 

3. 1.1. 2 BIGG (MATRIX) 

Description 

[BTGGgg] = Coriolis acceleration coefficient matrix - g set. 

3. 1.1. 3 MIGG (MATRIX) 

Description 

[MlGGgg] - Centripetal acceleration coefficient matrix - g set. 

3. 1.1. 4 M2GG (MATRIX) 


Description 


[M2GGgg] 


- Base acceleration coefficient matrix - g set. 



3. 1.1. 5 BASEXG (MATRIX) 
Description 


[BASEXGg] - Base acceleration matrix - g 'et. 

3. 1.1. 6 PDZERO (MATRIX) ’ 

Description 

F 

[PDZEROgJ - Load modification matrix in base acceleration problems - g set. 


3.1,2 Data Blocks Output from Module DUMMQD2 
3. 1.2.1 FRL (TABLE) 

Description 

Frequency Response List 

The FRL output by DUMflODZ contains one logical record. This logical record 
contains a sorted list of frequencies. 


Table Format 

Record Word Type 

0 1.2 BCD 

3 I 

1 1-m R 

2 

Table Trailer 

Word 1 = 1 
Word 2-6 = zero 

3. 1.2. 2 FOL (TABLE) 

Descripticn 

Frequency Response Output List 
Table Format 

Record Word Type 

0 1-2 BCD 

3-NFREQ+2 R 

1 


Item 

Data Block Name 
Set ID = 1 

Radian frequencies (w = 2ttF) 
End-of-file: 


Item 

Table Name 

Frequencies F (u = 2^iF) 
End-of-file 
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Word 1 = Number of frequcnci^Ss (.Ni-flEQ) 

Word 2 = Frequency set record lUimber (=1) 

Word 3-6 - zero 

3. 1.2. 3 REORDER! (MATRIX) 

Description 

[REORDER! 1 - Load reordering matrix in time-dependent frequency response problems 
for cyclic structures. 

Matrix Trailer 

Number of columns = y!TSTEPS*FKMAX, if CYCIO = -! 

Number of rows = J NTSTEPS*NSEGS, if CYCIO = +! 

Form = square 

Typi = real singer precision 

3. 1.2. 4 RE0RDER2 (MATRIX) 

Description 

[RE0RDER2] - Load reordering matrix in time-dependent frequency response problems 
for cyclic structures. 

Matrix Trailer 

Number of columns = ■) FLMAX*FKMAX, if CYCIO = -1 

Number of rows = J FLM,AX*NSE6S, if CYCIO = +1 

Form = square 

Type = real, single precision 



3.Z ■ FuUCYIONAL MGODLcS 

3.2.1 Functional Module DUMilODl 

3. 2. 1.1 Entry Point : DUMODl 

3.2.1 .2 Purpose 

To generate the Coriolis, centripetal and base acceleration coefficient , 
matrices and the base acceleration matrix for a forced vibration response 
analysis of rotating structures. 

3. 2. 1.3 DHAP Calling Sequence 

DUMMODl CASECC, BGPDT, CS7M, DIT, FRL, HGG,,/FRLX, BIGG, MIGG, M2GG, 

BASEXG, PDZER0,,/V, N. N0MG6/V, Y, CYCIO/V, Y, NSEGS/V, Y, 

KMAX/V, N, FKMAX/V, Y, BXTID=-1/V, Y, BXPTID=-1/V, Y, BYTIO=-l/ 

V, Y, BYTID=-1/V, Y, BYPTID=-1/V, Y, BZTID=-1/V, Y, BZPTID=-1/ 

V. N. HOBASEX/y, N, h'OFREQ/V, fl, OMEGA S 

3. 2. 1.4 Input Data Blocks 
CASECC - Case Control . 

BGPDT - Basic Grid Point Definition Table., 

CSTM - Coordinate System Transformation Matrices. 

, DIT - Direct Input Tables 
FRL - Frequency Response List (radians). 

HGG - Partition of mass matrix (g-set). 

Notes: 1. All input data blocks can be purged if only parameters FKMAX and 

NOBASEX are to be computed. 

2. CASECC, DIT and FRL can be purged if output data blocks FRLX and 
BASEXG are purged. 

3. 2. 1.5 Output Data Blocks 

FRLX - Frequency Response List (modified) 

BIGG - Coriolis acceleration coefficient matrix (g-set) 

HIGG - Centripetal acceleration coefficient matrix (g-set) 

M2GG - Base Acceleration coefficient matrix (g-set) 

BASEXG - Base acceleration matrix (g-set .x f) 

PDZERO - Load modification matrix in base acceleration problems (g-set x f) 
Notes." 1. All output data blocks can be purged if parameter N0MGG=-1. 

2. BIGG and MIGG can be purged if N0K6G=-1 or if 0HEGA=0.0. 
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3. FRLX and Pn?rRC..c?'-. b? if 0ME6A=0.0. \ 

4. FRLX, PDZERO, !?2GG and 3ASEXG can be purged if N0MGG=-1 or if 
N0FREQ=-1 or if CYCiO=i-l or if all three parameters BXTID=SVTID= 

‘ BZTID=-1. 

3. 2. 1.6 Parameters 

NOHGG - Input-integer-no default. M6G v.'as not generated if NOHGG®-!. 

CYCIO - Input-integer-no default. This parameter specifies the form of the 
input and output data from cyclic structures. A value of +1 is used 
to specify physical segment representation and a value of -1 for 
cyclic transformation representation. 

NSEGS - Input- integer-no default. The number of identical segments in the 
structural model. 

KI^AX - Input-integer-no default. KflAX specifies the maximum value of the 
harmonic index. The maximum value that can be specified for KI-lAX 
is NSEGS/2. 

FKfiAX - Output- integer-no default. FKMAX is a function of KJiAX. 

NOBASEX - Output-integer-no default. H0BASEX=-1 if data block BASEX6 is not 
generated. 

NOFREQ - Input-integer-no default. NOFREQ—1 ;if FREQUENCY was not selected 
in the Case Control deck. 

OMEGA - Input-real-no default. Rotational speed of the structure in radians. 
OMEGA = 27T.RPS. 

BXTID - Input-integer-defaults. The values of these parameters define the 
BYTID set identification numbers of the TABLEDi Bui k Data cards which define 
BZTID the components of the base acceleration vector. The tables referred 

BXPTID to by BXTID, BYTID and BZTID define magnitude (LT-2) and the tables 

BYPTID referred to by BXPTID, BYPTID and BZPTID define phase (degrees). The 

BZPTID default values are -1 v/hich means that the respective terms are 

ignored. 

3. 2. 1.7 Method 

Parameters NOBASEX and FKf'lAX are computed depending on the values of various 
input parameters. Parameter NOBASEX is set equal to -1 if parameters N0MGG=-1 
or CYCI0=t-l or N0FREQ=-1 or if parameters BXTID=BYTID=BZTID=-1 , othenvise NOBASEX 
is set equal to +1 indicating that base acceleration data blocks are to be 
generated. 
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If parameter CYCI0=-1 . then parametei^ FhW\X is computed as follows. If NSEGS 
is odd then FKHAX=2*KHAX+i i if NStGS is even and KHAX=NSEGS/2, then FKJ-1AX=NSEGS, 
otherwise FKfiAX=2*Kf-1AX+l . 

If parameter M0MGG=-1 then no data blocks are generated and an exit is made 
from module DUH’iODl , othenvise computations proceed in throe phases. In the first 
phase BIGG and M166 are generated unless parameter 0MEGA=0.0, M2GG is generated 
if parameter N08ASEX=+1 . The second and third phases generate data blocks asso- 
ciated with base acceleration problems and are only executed if f'iOBASEX=+l . In the 
second phase FRLX and PDZERO are generated unless parameter 0MEGA=0.0. Data 
block BASEXG is generated and output in phase three. 

3. 2. 1.7.1 Phase 1 - Generation of BIGG, HIGG and M2GG 

Phase one begins with a request for open core and buffer allocation. If 
0MEGA=0.0 then BIGG, anf M166 are not output and their buffers, IBUF3 and IBUF4, 
are not allocated and IBUF5 is set equal to IDUF3. If coordinate system trans- 
formations exist then the C3TH data block is open and the coordinate system in- 
formation is placed in core and readied for use by subroutine PRETRO. 

The primary loop in phase one is controlled by the number of grid points in 
the Basic Grid Point Definition Table (BGPDT), scalar points are not allowed by 
DUMKOOl . Each grid point in the BGPDT is considered in order and the corresponding 
columns of the mass matrix, MGG, are processed to form BIGG, MIGG and M2GG. When 
all grid points have been processed the necessary trailers are written. For the 
ith grid point in the BGPDT the corresponding translational terms of MGG are un- 
packed and the diagonal terms are isolated into a 3 x 3 matrix [m|]. If the 
grid point is not in the basic system then subroutine TRANSD calculates the 3x3 
transformation matrix [T^.] from global coordinates to basic coordinates for the 
grid point and [m|] is transformed to the basic system to form [M-]. The average 
of the three diagonal terms of [M-3 is then used to form [BTj], [MT.j and [M2.J. 

I III 

These three submatrices are then transformed back to the global coordinate system, 
if necessary. The 3x3 matrices [BT^], [mT^] and [M2j] are then packed into the 
BIGG, MIGG and H2GG matrices. 


(a) 


[MGGJ 

gxg 




0 

CM„3 


where n = the total number of grid 
points. 
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where 


and 


[H..] = 
. ^6x6 


[m]j = 

'3x3 


i 


ior v' ^' 


r“ 4 


't - 

! r»2 
i 




^T1 

m I ‘ 0 0 

0 rnT^ . 0 


0 


r3 


Or T<jQTi 


^b) Transform [M^j from global to basic coordinate system 


(c) Compute average of [M^.] 


■3 . -K . 

" _ r ^ wnere the mass (in the basic coordinate system) 

’ k-1 ^ point of the total of 'n' nodes 

3>o the direction. 


in 


(d) Form BIGG 


[BIGG] = 
gxg 


[Bl,] 


[Bio] 


0 


[Bln] 


where 


[BIJ 

^6x6 


[Blj] ! 

h 1 - - -I 

; [0] 


and 


[Blj] = [T,]^ 

’3x3 ' 


0 0 0 

0 0 -m. 

0 iii. 0 


[T,] 


(e) Form HIGQ 


[HIGG] 

gxg 


[Ml,] 


[Mlpl 
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where 


: C'"i3 


6x6 


t M-J » 1 r 
i ■ 


i t“ii 


OrtrOfiVAL FAG-' fw 
POOR 


and 




lT \r. 


0 0 


n 


[HI j = [TJ‘ O n. 0 

3x3 ' ’ « 

0 m. 


a^i 


(f) Form M2GG 


[M2GG] 

gxg 


v/here 


[M2J 


'6x6 


[H2J 


[M2J 


iM2|3| I 


[ 0 ] 


CH2n] 


and 


[M2h = [T.f 
’3x3 ’ 


m\ 0 0 

0 m. fn. 

0 




3. 2. 1.7. 2 Phase 2 - Generation of FRLX and PDZERO 

In this phase the FREQUENCY set selected in the Case Control deck is located 
in data block FRL and stored in core. If parameter 0MEGA=0.0 of if parameters 
BYTID=-1 and BZTID=-1 then phase tv/o is complete, otherwise phase two processing 
continues. The Frequency Response List must be modified to include an expanded 
set of frequencies. Read and copy from FRL to FRLX record 0 and all logical 

records up to the selected frequency set. The only set that will be modified in 

FRL is the selected frequency set. Once the set of selected frequencies have 

been found and stored in core a vector for FRLX and PDZERO are generated using 

the FRL frequencies stored in core and parameter OMEGA. Lot for i = 1, NFREQ 
be the frequencies (in radians) from FRL. 

If w^./O.O, create 3 entrios, 0.0, 1.0 and 0.0 for PDZERO and create 3 entries, 
1 u.- 0MEGA|, oj., and loi.fOMEGAl for FRLX. 
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If = 0.0, create 2 entries I -C and G.G for PDZERO and create 2 entries, 0.0 
and Iomega! for trlx. 

After the expanded list of freciucncier. is generated call routine DUMOIE to sort 
it in ascending order. DUMOIE also returns a sorting index so other vectors may 
be sorted the same as FRLX. Sort PDZERO using this sorting index. Output this 
FRLX vector and continue copying tfio romstning records of FRL to FRLX. Output 
data block PDZERO by writing out the PDZERO vector FKMAX times, thus creating 
FKflAX columns. The original unexpanded frequencies from FRL and the sorting 
index stored in core are retained for phase 3 processing. 

3.2. 7.3 Phase 3 - Generation of BASEXG. 

If NCBASEX=-1 then this phase is skipped, othenvise processing continues. 

A unique list of table IDs using parameters BXTID, BYTID, 3ZTID, BXPTID, BYPTID 
and BZPTID is generated and a call to PRET.AB is made so that tables TABLEDl , 

TABLED2, TABLED3 and TABLED4 can be interpolated by calls to TAB. Routines DUMOIA, 
DUMOIB, DUMOIC and DUMOID are used to generate data block BASEXG. Routine DUMOIA calls 
the routines to generate the BASE table and outputs the BASEXG matrix. The BASE 
table is used to generate up to three groups of flFREQX columns, where NFREQX is 
the number of expanded frequencies from phase two, in the BASEXG matrix. Routine 
DUMOIB is called to generate the BASE table if the original FRL frequency list 
was not expanded, see phase two, otherwise routine DUMOIC is called. Routine 
DUMOID sorts the columns of the BASE table so that they are arranged in the same 
order as tlie modified frequency set if FRLX was generated in phase two. The 
following is a mathematical description of matrix DASEXG. 

(a) Let Xg(f^.), 0j^(f^). Y^(f.), 0y(f.), Z^j(f^.), 02(f,-) be Input via frequency 
dependent tables TABLEDi where the table IDs are defined by parameters 
BXTID, BXPTID, BYTID, BYPTID, BZTID and BZPTID respectively. X^, Y^ and 

Z are magnitudes in L7-2 units while 0 , 0 and 0 are phase angles in degrees. 

O A jr Z . 

(b) Define control flag M0DFRL. 

If parameter 0MEGA=0.0 or parameters BYTID=-1 and BZTID=-1 then set MODFP.L 
to false, otherwise tIODFRL is true. 

(c) Let FRL be a vector of NF frequencies (in radians). 

rP.L - , ^ 2 ^ “^3’ ''' ^NF^ 

(d) If MODFRL is false then generate complex base table BASE of order 3 x NF. 

[BASE] = [{3ase(f,)} • • • {BASE(f,,J} J 

3;<NF .3x1 3x1 


■7 

9* 
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where f| = w./2tt for V- 1, 2, .... NF 
and 


OF pOOt? QUAt. J 


{BASE(f5)} 

1 o 


3x1 






3x1 


(e) If HODFRL is true then generate complex base table BASE of order 3 x NFX 
where NFX is an expanded number of frequencies as defined below. 

[BASE] = r[BASE(f.)3 [BASE(fp)j .... [BASE(f,.p)]l 
3xNFX L ’ Hr j 

where f^ = u^/2tt for i -• 1 , 2, . . . , NF 

and each {BASE(f.)} is either 3 x 2 if w- = 0.0 or 3 x 3 if w. 0.0. 

(e.l) If oj. = 0.0, then [BASE(f.)] is defined as follows: 

’ ^ 3x2 


CBAS'^(f.)3 
* ^ 


3x2 


where 


A 0 
0 B 
0 C 


SGN = 1.0 if parameter OMEGA ^0.0, otherwise SGN=-1.0 and 
A = XQ(f,.)-e'^^^’' 

B = Y^(f.) •cos(0y(f.)) - i.SGN.Zjf^) •cos(0^(f.) 

C = Zg(f^)vcos(0^(f^)) + i-SGN-Yp(f^)-cos(0y(f.)) 

(e.2) If 0.0, then [BASE(f.)3 is defined as follows: 

’ 3x3 


[BASE(f.)3 
* ^ 


3x3 


0 A 0 

8 0 C 

DOE 


where 


and 


SGNA = 1.0 if (w^ - OMEGA) > 0.0, otherwise SGNA = -1.0 
SGNB = 1.0 if (w^. + OMEGA) i.0.0, otherwise SGIB = -1.0 
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ORiQiWA- PACiS 13 

' OF ^'OO^^QUAUTy 

C-. 0.5 

0 - 0.5 . [SGOA . V„(fi) . ,i-S0SA(S,(f,)-S00A-t) , •SGm.S,(fi )j 

E = 0.5 . [-SGHB.Y„(f,.) , ,f-SG»(0,(f,)-SGHi!. f) , . ^LSONB-OzCf ( )j 

(f) Define the cciiiplex base acceleration matrix SASEXG of order G x {NF*FKfiAX) 
as follows: 

Let NF be the number of frequencies in the BASE matrix, i.e., let MF = HF 
if MODFRL was false or NF = NFX if MODFRL was true. 

[EASEXG] =r[BASEXGh [BASEXG^] [EASEXG^] . . . [BASEXg'"'^'’^'’^] j 

qx(NF-FKHAk) qxNF gxNF qxNF gxNFj 


where 


[BASEXG^J 

gxNF 


[BASEX^l 

6xNF 

[BASEX’3 

6xNF 


[BASEX’] 

5xNF 


for i = 1 , 2 and 3 


and 


[BASEXG^] = [0] for i = 4,5,6,..., FKflAX 
gxNF 

NOTE: [BASEX\l is repeated N times where N ='g/5 and g is the g-set size. Scalar 

points are not allowed so each node has 6 degrees of freedom. 


[BASEX^] 

6xNF 


BASE(1,1) BASE(1,2) .... BASE(1,NF)‘ 


0 

0 

0 

0 

0 


0 

0 

0 

0 

0 


0 

0 

0 

0 

0 
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[BASEX^J 

6xNF 


.0 : . .. 0 ... 0 ■ 
BASE(2^V) BASh(2,;:!} BASE(2,NF) 

BASEOVl) BAS[-:(3,2) BASE(3,NF) 

0 - O ' ■ 0 

0 0 0 

0 0 0 


[BASEX^J 

6xNF 


■ 0 0 
BASE(3,1} BASE{3,2) 

-BASE(2,1) -BASE(2,2) 

0 
0 
0 


0 

0 

0 


0 

BASE(3,NF) 

-BASE(2,NF) 

0 

0 

0 


3. 2. 1.8 Subroutines 

Utility subroutines GKMATO, PRETRD, TRANSD, PRETAB and TAB are used. 
See subroutine descriptions. Section 3 of NASTRAN Programmer's flanual. 


3. 2. 1.8.1 Subroutine Name; DUMOIA 

1. Entry Point: DUMODIA 

2. Purpose: To define and output the complex single precision base 

acceleration matrix BASEXG. 


3. Calling Sequence: Call DUMOIA (BASE, BASEl , Z, W, BUF, INDEX, MODFRL, 

BASEXG, NROW, NF, NFX, FKMAX, OMEGA) 

BASE - Storage for BASE matrix - complex S.P. - input. 

BASEl - Storage for sorted BASE matrix - complex S.P. - input. . 

Z - Storage for one column of matrix BASEXG - complex S.P. - input. 

W - Frequencies (radians) from data block FRL - real - input. 

BUF - GINO buffer for BASEXG - real - input 
INDEX - Sorting index - integer - input 

MODFRL - Flag to indicate if frequency list was expanded - logical - input 
BASEXG - GINO file number of BASEXG -- integer - input. 

MROW - G-set size - integer - input. 

UF - Number of frequencies in FPL data block - integer - input. 

NFX - Expanded number of frequencies - integer - input. 

FKMAX - Function of parameter KMAX - integer - input. 

OMEGA - Rotational speed of structure in radians - real - input. 
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3. 2. 1.8. 2 Subroutine Nams: DUMOIB 

1. Entry Point: D'JMOID 

2. Purpose; To define the complex single precision BASE matrix used in 
generating the cwnpSete bass acceleration matrix BASEXG. This routine 
is only called if MODFRL is false. 

3. Calling Sequence: CALL DUMOIB (BASE, W, fJF) 

BASE - BASE matrix - complex S.P. - output 
W - Frequencies from data block FRL - real (radians) - input. 

NF - Number of frequencies in 1.' - integer - input. 

COMMCN/CONDAS/PI, TF'OPI, RADEG, DEGP.A, S4PISQ 
C0HM0N/BLAKK/DUM(5), 3XTID, BXPTID, BYTID, BYPTID, BZTID, BZPTID 

3. 2. 1.8. 3 Subroutine Name: DUMOIC 

1. Entry Point: DUMOIC 

2. Purpose: To define the complex single precision BASE matrix used in 
generating the complete base acceleration matrix BASEXG. This routine 
is only called if HODFRL is true. 

3. Calling Sequence: CALL DUMOIC (BASE, W, OMEGA, NF) 

BASE - BASE matrix - complex S.P. - output. 

W - Frequencies from data block FRL - real (radians) - input. 

OMEGA - Rotational spend of the structure in radians - real - input 
NF - Number of frequencies in 11 - integer - input. 

COMMON/CONDAS/PI, TWOPI, RADEG, DEGRA, S4PISQ 
C0MM0N/BLANK/DCM(5), BXTID, BXPTID, BYTID, BYPTID, BZTID, BZPTID 

3. 2. 1.8. 4 Subroutine Name; DUMOID 

1. Entry Point: DUMOID 

2. Purpose: To sort the columns of matrix BASE in the same order as the 
expanded frequencies in data block FRLX. 

3. Calling Sequence; CALL DUMOID (BASE, BASE!, INDEX, NFX) 

BASE - BASE matrix - complex S.P. - input/output 
BASEl - Temporary storage used for sorting matrix BASE - complex 
S.P. - input. 

INDEX - Sorting key - integer - input 
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NFX - Numtei; of columns of- matr-x 3ASE and leiigth of INDEX - 
integer “ fhpiit^ . , 

3. 2. 1.8. 5 Subroutine Name: DlifiOlE 

. 1. Entry Point: DUKOiE 

2. Purpose: To sort the list of expanded frequencies of data block FRLX 
and to supply an index key so these vectors can be sorted the same way 

3. Calling Sequence: CALL DUM01E(A,K,N) 

A - Vector to be sorted - real - input/output. 

K - Sort index key - integer - output 
N - Length of A and K 

3 . 2 . 1 . 9 Design Requirements 

a) Open core is defined at /DUMIXX/ 

b) No scratch files are used 

c) DUMODl resides in LINKNS07 

d) Open core for five GINO buffers is needed. 

e) The loyout for open core is as follows: 

Phase I 

. COHM0N/DUMIXX/ Z ’ , '.i 


Z(1) 

Column of "iGG 

NTVPE*G-set 


FREE 


Z(ICSTM) 

CSTM DATA 

LCSTM 

Z(IBUF5) 

M2GG 

GINO BUFFER 

Z(IBUF4) 

HIGG 

GINO BUFFER 

Z(IBUF3) 

BIGG 

GINO BUFFER 

Z(IBUF2) 

BGPDT 

GINO BUFFER 

Z(IBUFl) 

CSTM/MGG 

GINO BUFFER+1 
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f*7* 


f 


f 


I 

I 
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1 

i 
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Phase II 


COMMON/OUMIXX/2 


Z(IFRL) 
Z( INDEX) 
Z(irRLX) 
Z(IPDZ) 


Z(IBUF3) 

Z(IBUF2) 

Z(IBUFl) 


FRL DATA 


SORT INDEX KEY 


FRLX DATA 


PD2ERO BATA 


FREE 


PDZERO 


CASEI.L/FRLX 


FRL 


Phase III 


COMMON/ DUM I XX/Z 


Z(IFRL) 

Z( INDEX) 

Z(ITAB) 

Z(N1) 

Z(N2) 

Z(N3) 


FRL DATA 


SORT INDEX 


PRETAB TABLE DATA 


BASE MATRIX 


BAS El MATRIX 


COLUMN OF BASEXG 


FREE 


Z(IBUFl) i DIT/BASEXG 


CRlQiivA- 

OF po-oi? QUALFTf. 


NF 

3kNF 

3)frflF 

3$NF 


GINO BUFi^ER 
GINO BUFFER 
GINO BUFFER+1 


NFS 

3*NFSX 
NTABL 

(3*NFSX)*2 ') 

C 

(G- 


[3*NFSX)*2 y 

;G-set.)*2 J 


GINO BUFFER+1 


Complex 


3.2.1.10 Diaqnos \ c Messages 

The following fatal error messages may occur: 
3001, 3002, 3003, 3008 and 3031. 


53 
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3.2.2 Functional Module DUMMQD2 

3. 2. 2.1 Entry Point : DUM0D2 

3. 2. 2. 2 Purpose 

To generate tables FRL and FOL and matrices REORDERl and RE0RDER2 to be 
used in a forced vibration response analysis of rotating cyclic structures. 
Parameters LflAS, fITSTEPS, FLMAX, NOROl and N0R02 are also computed. 

3. 2. 2. 3 DMAP Calling Sequence 

DUMMOOE TOL,,,,,,,/ FRL, FGL, REORDERl, RE0RDER2, , ,,/V,Y,NSEGS/V,Y, 
CYCIO/V,Y,l.MAX=-l/V,N,FmAX/V,N,FLMAX/V,N,NTSTEPS/V,N, NOROl/ 

V,N,NCR02 S 

3. 2. 2. 4 Input Data Blocks 

TOL - Time output list. . 

. NOTES: 1. TOL must be present 

3. 2. 2. 5 Output Data Blocks 

FRL - Frequency Response List 

FOL - Frequency Output List 

REORDERl - Load reordering matrix for time-dependent frequency response problems 
RE0RDER2 - Load reordering, matrix for time-dependent frequency response problems 

NOTES: 1. FRL and FOL cannot be purged. 

I , 

3. 2. 2. 6 Parameters 

NSEGS - Input-integer-no default. NSEGS is the number of identical 
segments in the structural model. 

CYCIO - Input-integer-no default. The value of this parameter specifies 
the form of the input and output data for cyclic structures. A 
value of +1 is used to specify physical segment representation 
and a value of -1 for cyclic transformation representation. 

LHAX - Input/output-integer-default. LMAX specifies the maximum time- 
harmonic index for cyclic structures. The default value is 
NTSTEPS/2, where NTSTEPS is defined below. 

FKf'lAX - Irout-integer-no default. FKMAX is a function of parameter KMAX. 
FLHAX - Output-integer-no default. FLMAX is a function of parameter LMAX. 
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NTSTEPS 

- Output- irivener-rto default. 

The number of time steps from data 

i : 


block TOL. 

• • 1 

■ ; . . ■ ! 

i • 

NOROl 

- Output-intGUGu-no default. 

T10R01=-1 if matrix REORDER! is not 

i ^ 


generated, !-l other.-.'ise. 


I 

M0R02 

- Output- integcr-no defualt. 
generated, +1 otherwise. 

'■i0R02-“l if matrix RE0RDER2 is not 


J 

I 


T 

i 

I 

j 

4- 

T 

I 


1 
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3. 2.2. 7- Method 

Computations proceed In tiiree phases. Parameters NTSTEPS, LMAX and FLMAX 
are computed 1n Piiase I. Data blocKS FRL and FOL are generated and output in 
Phase II -and matrix data blocks P.F0RDER1 and RE0RDER2 and their respective 
parameters NOROl and N0R02 are generated, and output in Phase III. 

3. 2. 2. 7.1 Computation of Parameters NTSTEPS, 1.MAX and FLHAX 

Data block TOL is open and the list of output times is read from the header 
record and stored for use by Phase II. Let. NTIMES be the number of times read. 

a) Parameter NTSTEPS 

If CYCI0=-1, then tiTSTEPS=(MTIMES*FKnAX)/FKMAX 
If CYCI0=-M, then NTSTEPS-(HTIMES*NSEGS)/NSEGS 

b) Parameter Lf-lAX 

If LMAX<0, then the default value of LMAX is set equal to NTSTEPS/2. 

c) Parameter FLMAX 

If NTSTEPS is even and LMAX=NTSTEPS/2, then FLHAX=NTSTEPS, otherwise 
FLMAX=2*LMAX+1 . 

3. 2. 2. 7.2 Generation of tables FOL and FRL 

The list of times read from TOL are now converted to the frequency domain. 
The number of frequencies, NFREQ, is set equal to FLMAX. 

Let PERIOD = TIME(2) + TIME(NTSTEPS) 

then, FOL(l) = 0.0 

FOL(i) = (i-V2)*(^ •’^/PERIOD) for i = 2,4,6, ... ,NFREQ 
and rOL(j) = F0L{j-l) for j = 3,5,7,. ...NFREQ-l 

Data block FOL is then output and data block FRL is then generated om 
FOL by converting the FOL frequencies in hertz to FRL frequencies in ra .ns, 
FRL(i) = F0L(i)*2ii for i = i , NFREQ. 

3. 2. 2.7. 3 Computation of parameters NOROl and NOR02 and matrices REORDERl and 
RE0RDER2. 

REORDER 1 and RE0RDER2 are used for reordering columns of a matrix by post- 
multiplying the matrix whose columns are to be reordered. Routine DUM02A is 
called twice, once to generate and output REORDERl and once to generate and 
output RE0RDER2. See the subroutine description of DUM02A for details. 
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3.2. 2.8 Subroutines • f ^ ■ v 

DUM0D2 uses standard NASTRAIi GliiO Toutlnes and utility routines. 

3. 2. 2.0.'! Subroutine Name; DUH02A 

1. Entry Point: DUM02A 

2. Purpose: To generate and output column reordering matrices REORDERl 

and RE0RDER2 and to compute parameters NOROl and N0R02. 

3. Calling Sequence: CALL DUM02A(FILE, KKl , KK2, NORO, BUFFER) 

FILE - GINO file number of REORDERl or RE0RDER2 - integer - input. 

KKl - Reordering row index - integer - input. 

KK2 - Reordering column index - integer - input. 

NORO - N0R0=+1 if reordering matrix was generated, -1 otherwise - 
integer - output. 

BUFFER -GINO buffer - real - input 

4. Method; If KKl =1 or KK2 = 1 then set parameter N0R0=-1 , otherwise 
set parameter NORO=+lto indicate that the reordering matrix was 
generated. If N0R0=-1, then return otherwise continue processing., . 
Generate a real single precision reordering matrix of order KK1*KK2 by 
KK1*KK2. This matrix can be used to reorder columns of another matrix 
by post-multiplying the matrix vfnose columns are to be reordered. 
Column i of the reordering matrix contains a 1.0 in row j if column j 
is to become column i of the reordered matrix. For example, if 
column 5 is to become column 1 of the new matrix then the reordering 
matrix contains a 1.0 in row 5 of column 1. 

3. 2. 2. 9 Design Requirements 

a) Open core is defined at /DUM2XX/ 

b) DUM0D2 resides in LINKNS07 

c) No scratch files are used 

d) Open core for one BUFFER+1 is required. . . 
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The layout of open core 
COMMON/DUMZXX/ 


Z(ITOL) 

Z(IFOL) 


Z(IBUFI) 



3.2.2.10 Diagnostic Messages 


The following fatal error messages may. occur: 3001 , 3002, 3008, 3037 
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FORCED VIBPAflOrl ANALYSIS 


OP ROTATI 


!5G CYCLIC STRUCTURES 


Frequency Response of a'lS-Rladed Disc 
(Examples 1~5) by the Direct Method 


A. General Description 

Five inter-related examples are presented to illustrate various features 
of this nev; capability to conduct forced vibration analysis of rotating cyclic 
structures. A 12-bladed disc is used for demonstration. 


The capability includes the effects of Coriolis and centripetal accelera- 
tions on the rotating structure which can be loaded with: 

1. directly applied loads moving v;ith tiie structure, and 

2. inertial loads due to the translational acceleration of the axis of 
of rotation {'base* acceleration). 


Example 1 is conducted on a finite element model of the complete structure 
(Figure 1). Examples 2 through 5 use a finite element model of one rotationally 
cyclic sector (Figure 2). Results of example 1 are used to verify soitie of the 
results obtained in the remaining examples. Table 1 summarizes the principal 
features demonstrated by these examples. 


Steady-state frequency-dependent (sinusoidal) or time-dependent (periodic) 
loads are applied to selected grid point degrees of freed' m. The specified 
loads can represent either the physical loads on various segments or their cir- 
cumferential harmonic components. For illustration purposes only, the frequency 
band of excitation, 1700-1920 Hz, due to directly applied loads and base accel- 
eration is selected to include the second bending mode of the disc for a circum- 
ferential harmonic index k = 2. The 'blade-to-blade' distribution of the 
directly appl ied loads also corresponds to k = 2. Table 2 lists the first few 
natural frequencies of the bladed disc for k = 0, 1 and 2. Modes for k = ?. are 
shown in Figure 3. 

B. General Input 

1. Parameters: 

Diameter at blade tip = 19.4 in. 

Diameter at blade root = 14.2 in. 

Shaft diameter = 4.0 in. 
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Disc thickness - rt.2& 'in. : 


mads thickness 
Young's modulus 
Poisson's ratio 
Material density 


0.1 25 in. 

30.0x10^’ Ibf/in^ 

0. 3 

7.4 >t 10"'^ Ibs-sec^/in'l 


Uniform structural damping (g) = 0.02 


2. Constraints: 

All constraints are applied in body- fixed global coordinate system(s). 
All grid points on the shaft diameter are completely fixed. Rotational 
degrees of freedom 0^ at remaining grid points are constrained to zero. 
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A. Description 

This example useSitha direct rv'equoncy response capability in NASTRAN, 
RF8, and forms the basis to verify some of , the results of examples 2 through 5 

B. Input 

1. Parameters: 


Same as general input parameters. 

2. Constraints: 

Same as general input constraints. 

3. Loads: 

P(f;n) = ,A(f) cos (n'^.©-|| ) , , 

where n is the segment number, 

(Drepresent:. k = 2, 

O represents the total number of segments in the blatled 
disc. 

? is specified using RLOADi bulk data cards. 

C. Results - 


Sample plots of grid point displacement and element stress response are 
shown in Figures 4 through 6. The expected behavior about a k = 2 natural 
frequency of the bladed disc can be seen in all these figures. 
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n. Driver Packs and Bulk Data 
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FORCED VIBRATION ANALYSIS OF flOTCi^^p. rvr»,7r 
BLAOEO DISC EXAMPLE i SFULL HaiaLij-f.Evi: 5 

UNDEX 2C TYPE LOADS 


■'iCjORES 


OF POOR QUAUV: 


.iii I 


CARO 

COUNT 

I 


0 N T R 0 L DECK ECHO 


TITLE « FORCED I'l DR AH Oi4 ANALYSIS 0? ROTATING CYCLIC STRUCTURES 
SUBTITLE = BLAOED DISC EXAMPLE I 4FULL HQDEL.FREO fcOADSJ 

LABEL = XSNOEX 2C TYPE tOADS 
$ ■ . 

SPC = 30 
FREQ = I 
DLGAD =1 
OUTPUT 


10 

SET 1 = a, 22, 36, 50,1 

11 

16,30,44,53 

12 

1 S • 32 ,46 ,60 

13 

OLOAO = 1 

14 

DISP(SQRT2,PHASE) 

15 

STRESS! S0RT2. PHASE! 


DISP(SQRT2,PHASEI = ALL 
STRESS! S0RT2* PHASE! = ALL 
OUTPUrCXYPLOTl 

PLOTTER NASTPLT. .MODEL 0,0 
XPAPER = 8.0 
YPAPER =10.5 
XAXIS = YES 
. YAXIS = YES 

XGRID LINES = YES 

YGRID LINES = YES , 

CURVELINE SYMBOL = I , : 

YLOG = YES ! 

XTITLE = FREQUENCY (HERTZ) 

YTITLE = GRID POINT DISPLACEMENTS £ MAGNITUDE, INCH ) 

TCURVE = 14( T3RM) ,10 (T3RM) ,95{T3RMJ 
XYPLOT.XYPRINT OISP RESPONSE / 14 £ T 3RM » . 18 ( T 3RH J ,95( T3RM ) 

YTITLE = ELEMENT STRESSES ( HAGNIIUOE.PS I ) 

TCURVE = 11(31,11(5) ,IU7) .11110) , IK 12)411(14) 
XYPLOT.XYPRINT STRESS RESPONSE / 1 1 ( 3 ) .11 ( 5 ) . 1 1 ( 7 ) , 

lU 10), IK 12).11( 14£ 
TCURVE = 109(3) .109(5) ,109(7) ,109(10) ,1091 12), 109( 14) 
XYPLOT.XYPRINT STRESS RESPONSE / 1091 3 ), 109( 5) , 109i 7) , 

1091 10 » 6 1091 12).109£;14) 

BEGIN BULK 


SER INFORMATION MESSAGE 207, BULK DATA NOT SORTED, XSORT WILL RE-ORDER OECK^ 
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EXATiPLE 



This example uses the forced vibration capability with cyclic synmetry. 

The user input/output data for loads, displacements, stresses, etc,, pertain 
to the physical representation of the various segments of the bladed disc. 

The frequency-dependent applied loads correspond to k = 2, and hence the solution 
loops on the circumferential harmonic index k arc restricted to k -2 only via 
parameters KfilN and KJ1AX. 


B. Input 

1. Parameters: 

In addition to general input parameters, 

CYCIO = +1 physical cyclic input/output data 
KHIN = 2 minimum circumferential harmonic index 
KiMAX = 2 maximum circumferential harmonic index 
NSEGS = 12 number of rotationally cyclic segments 
RPS = 0.0 rotational speed 

GKAD = FREQRESP\ Specify the form in which the damping parameters 
LGKAD = +1 / are used. 

2. Constraints: 


Same as general input constraints. 

3. Loads: 

P"(f) = A(f) cos , 

where n is the segment number, 

(Z) represents k = 2, 

represents the total number of segments in the bladed 
disc. 

P is specified using RLOADi bulk data cards. 

C. Results 

Displacement and stress output results for selected grid points and elements 
are presented in Figures 7 through 10. Agreement between results of Figures 7-8 
and Figure 4, Figure 9 and Figure 5, and Figure 10 and Figure 6 is excellent. 
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d. DLOAC MUST BE USED TO DEFINE A FRCUUENCY CR T IME-DEPENDtNT 
LOADING UlNDiTICN FUR EACH SUBCASE. 

FUR FREObENCY-DEPENDENT LOADS, SUBCASES wITHCUT LOADS NEED NOT 
REFER TO A DLUAU CARD. 

FOR TIME-DEPENDENT LCaDS, SUBCASES WITHOUT LOADS MUST REFER TU 
A CLUAD CARD THAT GENERATES A NULL LOAD. 

V. AN ALTERNATE LOADING METHOD IS TU DEFINE A SEPARATE GRCUP OF 
SUBCASES FUR EACH HAHMUNiU INDL;;, K. THE PAKAMETER CYCIO IS 
INCLUDED AND TmE LOAD CCMPCNcNTS FCR EACH INDEX ARE DEFINED 
DIRECTLY WITHIN EACH GROUP FOR THu VARIOUS LOADING CUNDITIUNS. 
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i 
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i 

$ 

i 
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i 

i. 

$ 
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1. 

Z, 

3« 

4. 

!>« 


6. 


SUFURf 

EPQIN7 

SPG2NT 

CYJUIN 


tJULK 

8ULK 

BULK 

BULK 


DAIA CAROS 
UATA CAROS 
DATA CAROS 
DATA CAROS 


ARE !yjT ALLOWED. 
ARE iSlOF ALLOWED. 
ARE NOT ALLONEO, 
ARE REOUIKEO. 


IF A TSrtP CARO IS USED THEN IT MUST NOT 
UNLV UNE UNIFORM TIME STEP INTERVAL MUST 


UE COMINJEO SINCE 
BE SPEClrl cO. - 


THE SKIP FAClUR 
PARAMETERS USED 


FUR OUTPul» 
ARE - 


NCe Ut'l THE TSIEP CARO MUST BE 1. 


A. NSEGS 


B. CYCIO 


C. CYCSEQ 

D. CTYPE 

E. KMA,< 

F. KM IN 
C. UMAX 

H. NLOAO 

I. RPS 


J . BXTIO 
BYHO 
B2TIO 
BXPTID 


REQUIRED - THE INTEGER VALUE OF THIS PARAMETER 
iS THE .NUMBER OF lOENTiCAL SEGMENTS IN THE 
STRUCTURAL MODEL. 

REwuIREO ~ THE INTEGER VALUE OF THIS PARAMETER 
SPECIFIES THE FORM OF THE INPUT ANU OUTPUT OATA. 

A VALUE OF £2 IS USED TO SPECIFY PHYSICAL SEGMENT 
REPRESENIATICN. A VALUE OF -I Ib USED IQ SPECIFY 
CYCLIC TRANSFCRMATIUN REPRESENTAT ION. THERE IS NO 
OEFACLT, A VALUE MUST BE INPUT. 

FIXED - THE INTEGER VALUE OF THIS PARAMETER 
SPECIFIES THE PftUCEUURE FOR SE-UENCING THE 
EQUATIONS IN THE SOLUTION SET. THE VALUE OF THIS 
PARAMETER HAS UEEN SET TO -1 TC SPECIFY 
ALTERNATING COSINE ANU SINE TERMS. 

Fl/EO - THE BCD VALUE OF THIS PARAMETER 

DEFINES THE TYPE CF CYCLIC SYMMETRY. THE VALUE 
OF IHIS PARAMETER HAS BEEN SET TC -ROT- FOR 
ROTATIONAL SYMMETRY. 

REQdlREU - THE INTEGER VALUE UF THIS PARAMETER 
SPECIFIES THE MAXIMUM VALUE CF THE HARMONIC 
INDEX. IciERE IS NC DEFAULT FOR THIS PARAMETER. 

the maximum valce that can be specified is NSEGS/2 

OPTICiMAL - THE INTEGER VALUE OF THIS PARAMETER 
specifies THE MINIMUM VALUE OF THE HARMCMC 
INDEX TC BE USED IN THE SOLUTION LOOP. KMIN CAN 
EOLAL KMaX. the default VALUE IS 0. 

OPTICNAL - IHE INTEGER VALUE CF IhlS PARAMETER 
SPECIFIES THE MAXIMUM TIME HARMDML INDEX. THE 

IS NTSTEPS/2, WHERE NTSTEPS EQUALS 
uARDl PLUS 2. 

INTEGER VALUE OF IHIS PARAMETER 
OF LOADING CCNDITICNS. THE VALUE 
OF This parameter IS INTERNALLY CALLJLAIED. 
uPliCNAL - THE REAL VALUE CF THIS PARAMETER 
DEFINtS THE ROTATIONAL SPEED CF THE STRUCTURE 
IN REVOLUTiCNS PER UNIT TIME. THE. DEFAULT VALUE 
IS C.O . 

CPTIUNAL - the POSITIVE INTEGER VALUES OF THESE 
PARAMETERS DEFINE THE StT I DENT I F IL AT I UN NUMBERS 
UF THE TABLED! BULK DATA CARDS WHICH DEFINE THE 
COMPONENTS CF ThE BASE ACCELERAI ION V ECTUR. THE 


default, value 

N (FROM TSTEP 
FIXED - THE 
IS THE NUMBER 


} 1 
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h 

$ 
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$ 

S. 

$ 

i 

i 

$ 

$ 

$ 

t 

6 

I. 

$ 

$ 

$ 

'j> 

$ 

$ 

$ 
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i 

i 

i 

i 

i 
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I 
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5. 

« 


tlYPTIO 

3ZPTI0 


i\ • NUisPRT 


L. CKJPNT 


M. W7HASS 


N. CUUPHASS 

O. 


P • L Gi'i AD 


3. 


R. 


5. 


TAaCES REf-EREU 10 aV DXTID, BYTiD AfiD bZT 10 
DEFINE HAGWi TUOr:{iT“;i} .* AND IHE TABLES RLFcftED TO 
BY BXPTiUf EYPnO AND BZPTID DEFINE PHASEl DEGREE ) . 
THE DEFAULT VALUES ARE -1 WHICH MEANS THAT THE 
RE SPEC Ti VC TERMS ARE iGNOHED. 

OPTIONAL -- AN INTEGER VALUE UF Cl FOR THIS 
PARAMETER WILL CAUSE THE CURRENT HARMONIC INDEX. 
KINDEX, TO BE PRINTED AT THE TOP OF THE HARMONIC 
LOUP. THE DEFAULT VALUE IS Cl. 

UPTIONAL - A .'’CSITIVE INTEGER VALUE OF THIS 
PARAMETER WILL CAUSE THE GRID POINT WEIGHT 
generator to be executed and the resulting WEIGHT 
BALANCE INPOKHATICN TO BE PRINTED. DEFAULT IS -1. 
UPTIONAL - THE TERMS OF THE STRUCTURAL HASS 
MATRIX ARE HULTIPLItD BY THE REAL VALUE OF THIS 
PARAMETER WHEN THEY ARE GENERATED IN EHG. THE 
DEFAULT IS l.O . 

FIXED - ONLY LUMPED MASS MATRICES MUST BE USED. 
OPTIONAL - THE BCD VALUE OF THIS PARAMETER iS 
used to TELL THE GKAD NODULE THE DESIRED FORM UF 
MATRICES KDD. bOD AND MUD. IHE BCD VALUE CAN BE 
FREJRESP OR TRANRESP. THE DEFAULT IS TRANKeSF. 

NDTt P‘-MEM3Eii TO DEFINE PARAMETERS G, W3 AND WA. 
SL-. SECTION 9.3.3 (DIRECT DYNAMIC MATRIX 
ASSEMBLY) PAGES 9.3-7 AND 9.3-B OF THE 
NASIRAN TKEORETIlAL MANUAL. 

OPTiCNAL - THE INTEGER VALUE OF THIS PARAMETER 
IS USED IN CCNJUNCHCN WITH PARAMETER GKAD. IF 
GKAU=FRECRtSP THEN SEl LGKAD=l. IF GKAD=T RANRESP 
then set LGKAU=-l. THE DEFAULT VALUE IS -1. 
OPTIONAL - THE REAL VALUE OF IHIS PARAMETER IS 
USlD as A UNIFORM STRUCTURAL DAMPING COEFFICIENT 
it. THE DIRECT FORMULAT 1 ON UF DYNAMICS PROBLEMS. 
OPTIONAL - THE REAL VALUE OF THIS PARAMETER IS 
GSCJ AS A PIVOTAL f KEwUENCY FDR UNIFORM STKuCTUAL 

damping if parahcTer gkau-iranresp. in This case . 

WJ IS RECUIRED IF UNIFORMED STROCTUAL DAMPING IS 
DESIRED. THE DEFAULT VALUE IS 0.0 . 

OPTIUNAL - THE REAL VALUE Lr THIS PARAMETER IS 
USED AS A PIVOTAL FREwUENCV FOK ELEMENT STRUCTUAL 
DAPPING IF PARAMETlR GKAD= 1 RANRi.SP. IN THIS CASE 
IS RE I.UI RED IF STRUCTUAL DAMPING IS DESIRED FOP. 
ANY CF THE STRUCTUAL ELEMENTS. DEFADlI IS 0.0 . 


1 


» REMARKS - 
4 i. THE AHALYilS wil.l 
$ K INDEX = KMIN TO 

i 


LOOP 

KMAX. 


THRU A RANGE OF THE CYCLIC INDr.X, 
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I I 
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ALTER 3 1. 

EILE JX\/^=»APPEND/RUT=A^'PE^4iJ/PD=A^’PE^jD $ 

i PERFLfxM I.miAL ERROR CHECKS ON NSEOS ANO KMAX. 


CCND 
CCNU 
PARAM 
CCNU 
PARAM 
PARAM 
CONU 
$ SET 
PARAM //C.N.NOP 
♦ CALCULATE OMEGA, 
RARAMR //C,N,KPY 
PARAMk //C,N,HPV 
PARAMR //C,N,HPY 


ERR(jRCl,NSECA i IF USER HAS NCT SPECIFIED NSEGS. 

ERKUKCl.KHAX t IF USER HAS NOT SPECIFIED KMAX. 

//C,N,EO /V',N,CYCIO£RR /V,V,CYCI0=0 /C,N,0 $ 

EKRURCl.CYv. lULRK 5 IF USER HAS NOT SPECIFIED CVCIfJ. 

//C,ii,DIV /V,N,NSEG2 /V.Y.NSEGS /C.N.2 $ NSEG2= NSEGS/2 
//C,N,SUB /V,N,XHAXERft /V,N,NSEG2 /V,V,KHAX 5 
EKKORCl.RMAXtkK i IF KMAX ,GT. NSEGS/2 
DEFAULTS FOR PARAMETERS. 

/V,Y,NOKPKT=Cl /V»Y,LGKAD=-l S 

2#LMECA AND OHEGA’*-S‘2 FROM RPS. SET UEFAJLT KPS. 
yV.N. OMEGA /V,Y,RPS=0.0 /C.N.o, 283185 $ 

/V,N,OMEGA2 /C,N,2.0 /V.N, OMEGA S 
/V.N.OMEGASOR /V, N, OMEGA /V,N, OMEGA S 
i ULiNERATE NUKPS FLAG IF KPS IS ZERO. 

PARAMK //C,N,EU //V,Y,RPS /C,N.0.0 ////V,N,NORPS b 
$ MAKL SORL OOUPLED MASSES HAVE NOT DEEN REOUESTEDo 
PARAM //C,N, NUT /V.N.NOLUHP /V, Y,COUPMASS=-l b 

CCNu fcRRGRC2,NULLMP $ 

alter 21,21 i AOu SLT TO OUTPUT FOR TRLG. 

GEUK3,EUEXIR,GEuM2 / SLT,GPTT / V.N.NUGRAV 5 
SLT.GPTI 5 
b 

i SINLL MULTIPLE CONSTRAINTS ARE NUT ALLCriEU EXECUTE GP't NO»d SU THAT 
i more error UhC-CKS CAN oE MADE 8EFOKE ELEMENT GENERATION. 

$ ACC YS NEEUEU FUR PSF RtLOVERY IN SSG2, 

//C,N,MPY /V,^,NS^IP /C,N,G /C,N,0 i 

CASECC,oLOM^,EwEXlN,GPDr.BGPUT.CSTM/RG, YS ,USET .ASET/ Y, N. LUSET/ 
S,N,MPCfl/S,N,MPCF2/S,N,SINGLE/S,N.OMIT/S.N,REALl/S.N,NS\lP/ 

S, N, REPEA T/S.N, NO SET/S,N,NCL/S,N,NOA/C,Y,ASErUUT/S.Y,AUTOSPC F 
GM, CMD/MPCF l/UO,GCU/CMI T/KFS.PSf .CPC/SINGLE i 
GM.GMD.kG.GU.uUD.KFS.PSF .CPC.USE f , YS S 
UULK DATA IS NUT ALLOKED. 

//C»N,Nl,T /VtN.KEACDAIA /V.N.REACI $ 

ERKORL3,RLACuATA $ 

UPO NOW SO LhECKS CAN BE MACE. ADU TRL TO OUTPUT DAI A GLfiCKS. 


CP3 

LhKPNT 
ALTER 23 


PARAM 

GPA 


PURGE 
UrlKPNl 
J. SUPURT 
PARAM 
CCNC 

■i EXECgTE 


CPU 


i MUSI 

PARAM 

CCND 

3. ONLY 

PARAML 

PAR AML 

FAkAM 


UYNAM ItS.GPL.SlL.USET / GPLU , SI LU , USETU , TFPQUL .OLT . PSDL . FRL, 
TRL,,EuDYN / V , N , LUSE T /S , N, LUSET U/ V, N ,NGT FL/S , N, NOULT/ . 
S,N,NOPSDL/S.N,NUFRL/V,N,NCNLFT/S.N.NOTKL/V ,N,NO£6.U/C,N./ 
S.N.NOUE $ 

HAVE tITHLr. FREU OR TSTEP 8ULK DATA. 

//C.N.AND/V.N.FTEKR /Y.N.NOFRL /V.N.NOTRL S 
EKRQRCS.FTEkk b NO FREU CR TSTCP BULK DATA. 

(REGUEiNCY UR ISTEP IS ALLURED IN THE CA^E CCNTRUL 
CASFICC //G.N.UTI /C.N.l /C.N.14 //V .fUFREuSET i 
CASEXC //C.N.CTl /C , N , 1 /C,N,38 //V, N. T IMlSET i 
//C.N.MPY /V.N.FRLUnME /V.N.FREwSET /V.N.TIMESET S 
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f^ARAM //C.N.NOT /V.N.FTEKRI /\i vtufHhUU Hi & 

HAkAM //C,N ,LE /V,.M,NUFK£U /Vsfiji’REaSET /C .N»0 & 

PAHAM //C.N.Lt /V,N, NOTIME /V.N.TIMESET /C.N.O S> 

CQNb EHRURC6,rTERRl $ duTH fREO ANU TSTEP IN CASE CONTROL DECK- 
$ tPniNl cillLR data not ALLOWED 
f'AKAM //C.N.NCT /V»N, extracts /V»N,NQUE 5 
LONG EkRORCAfLXTRAPIS $ 

«. GENERATE DATA FOR LYCT2 MODULE* 

CPCYC oEdHAjEUDYN.USETU /CYCOD /Ve NtCTVPE=RUT /S»N,NDGu 5' 

CCNO ERRURCl.NUGC S 

LHkPNT LYCCO S 
ALTEh J2 5 

1. PRI-PURUL DATA aLOCkS THAT WILL NOT D£ DEfiERATEO 
PARaM //C,N,Ok /V.N.NOaMl /V.N.NCMGG /V.N.NORPS & 

PURGE olGG,MlGu /NCBMl i 
PURGE M2GGtM2bASLXG /NUHGG $ 

AL I Ek 3t* & 

i GENERATE DATA ELOLkS FKLX, 5IGG, MlGGi M2GG AND BASFGX, 
s. LtNtkATE PAKAMETEnS FkMAX AND NCDASEX. 

U'JMMUUl LAScCCtoGPOT.CirH.DI I,FRL.MGG»» / FRLX»81GGtMIGG, 

M2GG,bASEX&.PDZERC. , / V , N .NOHGG/ V, Y, L YC I C/ Y t Y # NSEGS/ 

V, Y,KMAX/StN.FKMAX/V,Y,8XTIU=-i/V,Y,0XPTlD=-l/ 
VfY,SYTID=-l/V.Y,aYPnD=-l/V,Y,BZTIU=-l/ 

Vf YtaZPriD=-l/StN,NUdASEX/V.N,NObR£y/VtNfbMEGA S 
PARAML FkLX //C .N t Pht SENuE ////y.NtNOFP.iX $ 

CuNG LbLFRLX.i'jOFkt.X S 

EyLIV FKLXtFKL i 
LABEL LBLFRLX i 

CHKPnT F KL r B iGG t M i GG 1 M2GG f liA S£ XG S 

Alter ^i2 $ 

PARAM //L.N.AOD /V.N.NOaGG /V.N.NDdMi /CtN.U S RESET NOQGG. 

ALTER fs2 $ 

i REGLFINE bGij AND KGG. 

CUNG LBLliA.NQBHl i 

PARAMk //C.N.COMPLtX // V,N,CPEGA2 / C.N.O.O / kf.N.LMPLXl i 
PARAMK //CtNfSL/B / V.N.HCMEGASG / CtNfU.O / V.N.OMEGASOk S 
PARAMR //C,NtCOHPLEX // V»NtMUK£GASG / C.N.O.O / Y.N.CMPLXZ 5 
AGC BGG.aiGG / Bbul / C ,N . ( i . 0, 0. 0) / V,N,CHPLX1 J. 

EGUIY bGGI«QGC i 

AGC kGGfMIGG / KGGl / C t N ♦ t i « 0 t 0* 0) / VtNtCMPLX2 i 

CULIV kC-GltkGG 

ChKPNT oGGtKGG :* 

LABEL LBLliA 

ALTER ij3,bs J> GP4 HAS BEEN MOVEU-UP. 

ALTER dB,60 i lPU HAS BEEN MOVED-UP. 

alter ilA 1. PARAM AND EGuIV LOGIC DEPENDING ON LGKAD. FOR FREU OR TRAN. 
PARAM //C,N,ANi)/V,N,KULkA/V,N,NuUE/V.NiNGk2PP 4 
LONG LGkAUlfLGKAU 4 .JRANCH IN NGT FRcCRESP. 

ALTER 115 S SEE ALTER il-t COMMENT. 

JUMP LGKAU2 S 
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LAHtL LGKAOl i | 

I EUUIV .M 2 PP,M 2 DD/NCA/a 2 PP,B 20 D/NUA/K 2 PP»Kaun/N 0 A/HAAtMUD/HD£MA/ 1 

I KAA,KxDD/KO£KA $ j 

CHKPNr K2PP»M2PPfD2PP,R2D0,M<il)De8200eK0D,MOi) $ i 

LAUEL LGKAD2 S \ 

ALTER 117,117 $ AOU PARAMETERS bRAU, w3 ANU ir.'i TO GKAO. ] 

CRAU USFTOtGM ,GQ ,KAA , oAA , MAA , KAAA ,K2PP ,M2 PP , B2PP/KDD, BilO. HDU.GH'J, I 

GOO, R2DO, M200, a2UD/C .Y,GKAO=TRANRESP/C,N,OISP/C,N, DIRECT/ ■ 

C,Y,G=0,0/C,Y,k3 = C'.C/CoV,w4 = 0.0/V,N,NOK2PP/V,N»NUH2PP/ ] 

V,N,iMC!B2PP/V,N,MRCFl/VsNrSIf;oLb/V,INi.bMIT/V,l\i,NCUE/V,r<!,KL]!<4G&/ | 

V,N ,NGabb/V,N,.SUEK2/C ,N,-1 $ | 

! Alter iia $ see alter lit comment. » 

j CUND LGKAD3,LGKAD J BRANCH IF NGT FREGRESP. . j 

AUER 119 S SEE ALTER 114 CCMMENT, 1 

I JUMP LGKAU4 i ■! 

j label LGKAU3 i f 

EwUlV b2UD,L.DU/NuGPHT/M20D,MUO/NOSIMP/K2UO,KOD/i<DER2 S t 

LAEEL LGKAJ4 S 

! ALTER 120, 123 $ \ 

i $ NCw 3ULGTILN LUGIL -I 

$ GENcKATE T IME-OEPcNUcN T LOADS IF TSTEP ,<AS REQUESTED IN CASE CONTROL. j 

' CuNC LBLTRLI .NOTIME 5 

i $ LUuP TH>nU all SUBCASES FOR TI ME-UEPENDENl LOADS. ; 

PARAM //C,N,HPY /V,N,REPEATf /C,N,1 /C,N,-l i •; 

, PAKAM //C.N.ADD /V.N.APPFLG /C.N.l /C,N,u $ I NI T I AL I2E FD K SDR 1 . ;• 

1 JUMP TRLCLOCP t . i 

* LABEL IRlGlUOP i i 

CASE CASECC,/CASEYY/C ,N, TKAN/S,N,REPEATT/S,N,NOLOOPl t i 

ChKPNT CASEYY $ 

PARAM //C.N.MPY /V.N.NCOL /C,N.O /C.N.l i '■ 

TRLG CASLYY.USETUvlJLT.SLT.bG PuT ,31L,CSTM,IRL,DIT ,GKD,GUn, ,EST, JIGG/ i 

, ,PDT l.PDl, .TCL/ V.N.NUSEI/S.N.PDEPDO/V.N.NCOL S I 

SDRl TKL tPDTl, I / ,PUT, /V .N, APPFLG/C . N . DYAAMI lS 6 I 

SDRl TRL.PUl / ,PD , / V ,N , APPFLG/C , N , DYNAMl CS $ i 

PAKAM //C.N.ADD /V.N.APPFLG /V,N,APPFLG /C.N.l t APPFLG=APPFLoC i. 

CCND TRLGJlNE .RcPEATT l •’ 

kept 7RLGL0UP . IOC 1 ; 

JUMP LKRCR3 $ j. 

LABEL TRLGDUNE J 

CHKPNT PUT.PO, IuL $ . * 

EQCIV P 0, PUT /PDEPDO i ‘ I. 

I lHnPNT PCT 1 

DJMNUDii TCL,,,.,,, / FRLZ ,fULZ .RECRDERl ,<\ECR0ER2 . . . , / 

' /.Y.NSEGi/V.Y.C YC IC/S,Y,LMAvt=-l/V,N,FKMAX/ 

S , N . F LH A X / S . N . iM r S TE P S/ S . N » NLiRCl / 5 »N . NCR C2 ■■ 

EUUIV. FRLZ.FRL // FOLZ.FCL i> | 

lHKPNT FkL . FCL , h EC^GcK 1 » RE CKDEK2 $ 

JUMP LBLFRL2 $ • 

j LABEL LBLIRLI i I 

i i GtNuiwWL FRECULNCY-DtPENOtNT LOADS IF FREQUENCY vJaS SELECTED IN CC. 
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. 

r£- 

!r 


PKLG 

CQNIJ- 

MPYAD 

EUOlV 

LABEL 


$ FORM NtW LCAUS. 


““ 5 S'RLlA:?AJRE«;cfN:F“ 8 T^ ■' / 

Rp“pDZEk“'''/‘'ppPS I ■'= f"*-* “S CENENAIIU. 

PPFXfPPF i 
LBLFKLXI * 


Co NO 

mpyal 

AUL 

fwUIV 

LCND 

itSGZ 

EJUIV 

LABEL 

LABEL 

LJUIV 

CHRP NT 


LBLFRLlfNObAStX S 
M2GG»aASLXG, / M2UAS5XG /C.AUO S 

PPFUPpff'^ ' '‘"“'"• 0 , 0.01 /C,., . (- 1 . 0 , 0.01 i 


/ vPUOUMioPSFl.POFl £ 


NLUAO = NF/FKMAX 


LBLbASEifNUSET 
USETO tGMO tVS,KF SfGOU » ,PPF 
PSFlfPSF // PDF I, PDF i 
LBLBAStl S 

LBLFRLI $ , 

PPF,PDF/NUSET i 

PPF,PSF,PUF ,FCL 4. 

* LOADS ARE rREQUcNC Y-OEPENDENI 

* PL'RFURM CYCLIC TRANbFuKMATI ON ON LOADS IF CYCIC=&1, 

t^LArri"! /C,N,l /V,N,PUFCGLS $ 

* CALLCLaTE the number of loads for CYCIC=-i. 

LuC /V.N.FKKAX 

CCND LBLPDONE ,CYC 10 i> i 

i CALCULATE THE NUMBcK OF LOADS FOR CYCIU=1. 

/VpN.NLCAU /V.N.PDFCOLS /V.Y.NSEGS £ NLOAU = NF/NScrs 

t v '''•“•'-"PE .'C.N.PUKE /V,r.NiECS.-l / 

COO / VtW.OLOAD /S.N.NCGG S 
ERRORC 1 1 NGU U £ 

PXF £ 

LBLPDuNE i 
LBLFRL2 i 
$ LOADS ARE TIME-DEPENDENT 

PARAH //CtiTfNOT / V # N cNU TC Yu 1 0 /V»Y,CYCIC 1 
£ BRANCH DEPENDING UN VALUE Cf CYCIC 
CCND LBlTRL 2,N0TCYLIU i 

i LYCICi=-l 

Eguiv PDT.PDTkZ l/NUROl i 

CCND L15LRU lA iNUROl $ 

MPYaD PCT.KEOkDERI , / PDTR2I / CtN,0 £ 
label LBLRGIA £ 

CYCTI •;\ff<^l*‘iLYCF2 /V.N.CTYPE/C.N.FCKS/V.NiNTSTFPS/ 

V» Y»lMA// V f N ,rKMAX/S .Nf NCGG £ 

Cr.ND ERRURCl.NOGL J. 

CHKPNT PXTK/1 i 

tUUiV PXTRZl.PXFZ i/NURu2 £ 

cGND LbLK02A,N0kU2 * ■ 

ilPYAD PXTRZl.KEORDEUZt / PXFZl /C,N,0 £ 

LABEL L5LRU2A £ 


CCND 

CFKPNT 

JUMP 

LABEL 
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1 n • 



EQUIv/ PXFZ1,PXP1 $ 

ChKPNT PXFl $ 

JUMP LbLTRL3 t- 

LAbEL LGLlRLa S. 

$ CYUlU = Cl ' 

MPYAU PDT,R£ORUERlf / P0TRZ2 / CtN.O $ 

CYCri POTR/..i /PXTrZZ?GCYCF 3 /VfNfCTYPfc/C.iNiFURE/V.NtNTSTEPS/V.VeLHAX/ 

V» YfNSEUS/SsNtNCUn h 
CCMO EkRORCI.I'iOuU f- 

ChKPNT PXTRZ2 S 

tUUn/ PXrRZ2»PXrR2/NURU2 $ 

CCjNO LBLR02B.NCRU2 i 

MPVAU PXTkZ2fRtOkUhR2, / PXTR2 /C.N.O S 

cAetL LBLk02ii $ 

CYCTL PXTR2 / PXfZ2rGCYCF4 / V.N.C TYPE/C ,N , FORE/ V . Y , NSEGS/ V . Y .KMAX / 
V,NfFLMAX/S»N,NGGO 5 
CONI) ERKUKLl.NUGG J. 

EQUIV PXFZ2.PXF I i. 

LHKPiNT PXFl i 

LABEL LfiLTRL3 $ 

i TIMt-UEPtNOENT LOADS ARh REAL. HAKE LOADS COMPLEX TU CORRESPOND 
i Tu FkbUUENCY DEPtNDbNT LOADS. ALSC SDR2 EXPECTS LOADS TO Bt COMPLEX 
•t IN FKEQRESP TYPE PRObLEMS; ... 

COPY PXFl / PXF2 i CONVERT REAL PXFl TO COMPLEX PXF. 

ADD PXF1.PXF2 / PXF / C , N . I D. 5 . I . 0) / C.N. ( 0.5 ♦-! .0 ) 4 

5CcrINtNL0AUFDKCYC12. 

PAKAM //CtNtADD /VaNtNuJAD /V»N«FLMAX /CaNtO NLOAD * FLMAX. 

LAbEL LbLPDCNE $ 

PARAM //CtNtADD /VtNtKINDEX /VtY,KMIN=0 /CtNtO 4 INTIT I AL IZE KINUEX . 

$ 

i initialise jxvf if nHIN is nut zero. 

$ 

PARAM //CtNtADD /VtNtKMlNL /V.YtKMIN /C,Nt-l S 

LCND NOKMINLtKMINL 4 

PARAM //CtNtADu /V.NtKMINV /CtNtO /CtNtO $ 

JUMP KMINlOUP t 

LAbEL KMINLOUP $ 

v,YCT2 CYCDOt t tPXF t t /ttPKFZtt / CtNtFORE/VtYtNSEGS/ 
VtNtKMINV/VtN,OVCSF.U/V,NtNLCAD/S»NtNL;GO i 
CCND ERKOkCltNOGC 4 

ADL PKFst / dKVFZ / CtNtlO.C.O.OJ $ 

CYCT2 CYCDDtt tlKVf Zt t /ttUXVFtt / C t N t BACK/ V t Y t NS EGS/ 
VtNtKMINV/VtNtCYLSEO/VtNtNLaAD/StNtNGGC $ 

LCND LR.RUKC 1 1 NOOC 4 

PaRAM //CtNtADu /V.NtRMlNV /VtNtKMINV /CtNtl 4 

RbPT KMINLGOPtKMlNL $ 

LAbEL NUKMINL * 

i 

JUMP iUPCVC n 

I ABEL TUPCYC 4 LOOP ON KINUbX 
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N A S T R A N 


■ 0R1G?MAL PACi‘5 sS’ 
OF POOR QUAUTVV 


G U T i V £ 


CONTROL 


DECK 


vCNu NOKPKT tNCKPRT i ' 

PHTPAkM //CcNeO /CtN.KlNDfcX S 

LAbtL NDKPKT i 

CYCV2 CVCOOsKUUtHOU, , , /KKKF ,MKKP , t . /C #N,f-ORE/V . Y. NSEGS / 
ViN,KINL)E.'</V.N,CYCS£U=-l/V»NtNLGAD/S»NtNOGC 6 
CCNU fcRRORCl,NfJGC 3> 

GHKPNT ^^UF,fl^KF t 

PARAH //C»NfSY3T //C.N.Sfi /C.N»2 $ KETHCD 31 IN CYCT2 PRODUCES 
$ UNUERFLUWS FOR PXF. USE METHOD 2. 

t-YCT2 CYCOi),HUD,,PXF», /BKKF,tPKF,» / C .N ,FLiRL/V t Y , NSEGS / 

V.N.MNUEX/V ,N,CYCSEQ/V.N.NLOAD/S,N,NCGC i> 

PARAM //C.NfSYST //t.,N.5a /C.N.O $ RESET HPYAU METHOD CONTRCL. 

CQNU EkKORCltNOGC 5 

CFKPUT BKKF.PkF * 

S SULJTIUN 


FKKD2 KRKF.bKKF .MKKF* ,PKF ,FGL / UKVF /C .N t 0. U/C » N tO .0/ G • N» -1.0 $ 
CHKPNl UKVF i 

CYU2 CYCUOif .UKVF, e /,,UXVF,, /C , N .BACK/ V ,Y , NSEGS/ V ,N,K U4DEX/ 
V.U.CYCSLQ/V.N.NLOAD/S.N.NOGC $ 

CUNO ERKOKG 1 .NOUU $ 

CHKPNT UXVF $ 

PAkAM //t.N.ADu /V,N,KlM)tX/V,N,KlNDEX/C,N,l 4 KINOEX - KINOEX G i 

PARAM //C.N.SUb /V.N.OGNE / V.V.KMAX /V.N, KINOEX 5 

CUNO LCYC2,UDNfc $ IF KINOEX .GT. KHAX THEN EXIT 

RtPT TGPcYC.lGC 4 

JUMP ERRORS i 

LAEEL LCYC2 4 

EwUIV UXVF.UDVr / CYCIU 4 

GhKPNT UDVF i 

GUNl» LCYC3.CYGIU * IF CYciO .GE. 0 THEN TRANSFuKM TO PHYSICAL. 

CYCTI UXVF / UDVF.GCYGol / V ,N ,cT YPE/C . N , BACK/ V, Y .NSEGS/ V . Y , KMAX/ 
V.N.NLOAU 4 
CHKPNT UOVF $ 

LABEL LCYC3 t 

CUNU LBLTRLA.NGTIML S 

EUUIV PXF.PUF2 / CYCIU 4 

CCNO LCYCA, CYCIU 4 IF CYCIU .GE. 0 THEN IKANSFURM TO PHYSICAL. 

CYCTI PXF / PDF2,GCYCa2 / V ,N .C TYPE/C , N , BACK/ V ,Y , NSEGS/V , Y . KHAX/ 
V.N.NLOAU 4 
LABEL LCYC^ 4 

4 IF LUAUS WERE T IME-uEPENOENT THEN RECOVER PPF AND PSf FROM PXF. 

SORI USLTDfrKL'r2f..uOD.GR0.t.» / PFFc*. /C.N.I /C.N.DYNAMI CS 4 

SSG2 USETD.GMD,YS,KFS,oUD.,PPFZ / .PGDUM.PSFZ .PLDUH 4 

LOU IV PPFZ.PPF // PSFZ.PSF S 

CHKPNT PPF.PSF 4 

LAtEL LBLTRLA 4 

Alter IZ^, IZA 4 use TOL instead of PPF to get output FREOJLNCY LiST. 
VOK CASLXX.EwUYNrUSETU.UOVF.FCL.XYCrtb./CJUDVCl./C.N.FREOKESP/U.N, 

DlRcCT/S,N,NGLOXT2/S.N.NCO/j,N.NUP/C,N.O 4 
ALTER 140,140 4 USE FUL INSTEAD OF PPF TU GET OUTPUT FHEOUENCY LIST. 





GR:ffiNAL 13 
Cr POOR OUALffY 


N A S T R A N EXECUTIVE CONTROL DECK ECHO 


SDK2 CASLXX«CSTM,MPT,UlT,EuDVN.SlLDe . .BGRDP.ELLeUHC.UPyCf EST.XVCD3, 

PPF/Of*PCl.GCPCl,UUPVCl,OfcSCi»C£f-Cl,PUPVCl/C.N.FREORESP/ 
5tN»N0SQKT2 S 

ALTER 160 S ADD LABEL f-GR ERRORS. 

LABEL ERRORS $ 

ALIER IbS, 166 $ REMOVE ERHORl ANO ERROR2. 


ALTER 168 $ EURCED VlbRATlCtM ERRORS 

LABEL ERkUHCI 5 CHECK NSEGSt KHAX AND OTHER CVCLIC DATA. 
PRTPARM //C.N.-7 /L ,N .C VC STAT 1 C S $ 

LABEL ERR0RC2 i COUPLED MASS NUT ALLOWED. 

PRTPARM //C.N.O /L . V .CU'JPMA SS S 
JUMP FINIS 5 

LAELL ERKOkCS $ SUPORT BULK DATA NUT ALLOWED, 

PKTPARM //C.N.-6 /C,N.UYCSTATICS S 

label ERRLRC'« S tPCINT BULK DATA NOT ALLOWED. 

PRTPAKM //C,N,0 /C.N.NOUE S 
JUMP FINIS J 


L ABLt. 
PRTPARM 
PRTPARM 
.iUMP . 
LABEL 
PHTPAKM 
PKTP Arm 
JUMP 


EHRORLB $ NtlTiTER FREw UK TSTEP WERE IN BULK DATA DECK.' 
//C,N,0 /C,N,NDFRL S 
//CeN,0 /C.N.NQTRL S 
MNIS A 

ERR0RL6 t BOTH FRtU AND TSTEP WtRE SELECTED IN CASE CONTROL. 
//CfNtO /C.NtNOFREQ $ 

//C.N.O /C.N.NOTIME S I 

F IN 1 S $ > : 


L'lf.ALTLR 


TIME 5 $ IBM 370/3031 

OIAG LA, 21 

CEND 
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rv'jR QUALl-nr 


CASS C 0 fi TT R 0 L DECK ECHO 


CARO 

COUNT 


1 

S 


i 

.'.-i 

2 

TITLE = FORCED VIBRATION ANALYSIS OF ROrATING CYCLIC STRUC’.URES 

3 

SUBTITLE = BLAOED OISC E5CANPLE 2 (CYC 

MODELS FREO LOAOScOi«YS£CAL 1/0 5 

4 

$ 


• 

5 

SPC * 30 


V 

6 

FREO = 1 


i 

7 

OUTPUT 


3. 

8 

SET 1 = 8al6. 

18 

■ ' ■ J 

9 

OLOAD - 1 



10 

DISPtSQRT2, PHASE) = ALL 


IL 

STRESS(50RT2» 

PHASE) = ALL 

1 

12 

SUBCASE 1 



13 

LABEL = 

SEGMENT 1' 

i 

i 

14 

OLOAD = 

1 $ FRED DEPENDENT 

LOADS i 

15 

SUBCASE 2 



16 

LABEL = 

SEGMENT 2 

} 

i 

17 

OLQAO a 

2 S FREO DEPENDENT 

LOADS 1 

18 

SUBCASE 3 


? 

19 

LABEL = 

SEG.'tCNT 3 

i 

20 

OLOAD = 

3 i FREQ DEPENDENT 

LOADS f 

21 

SUBCASE 4 


1 

\ 

22 

LABEL = 

SEGMENT 4 

] 

23 

OLOAD = 

4 $ FREO DEPENDENT 

LOADS ( 

24 

SUBCASE 5 


r 

« - 

25 

LABEL = 

SEGMENT 5 

! 

26 

DLOAO = 

5 $ FREO DEPENDENT 

LOADS ; 

27 

SUBCASE 6 


: 1 

28 

LABEL = 

SEGMENT 6 

1 

; 1 

29 

OLOAD = 

6 $ FREQ DEPENDENT 

LOADS ; 1 

30 

SUBCASE 7 


i 1 

' i 

31 

LABEL = 

SEGMENT 7 

: 1 

32 

OLQAO = 

7 S FREO DEPENDENT 

LOADS ; 

33 

SUBCASE 8 


! i 

34 

LABEL = 

SEGMENT 8 

i ! 

35 

OLOAD = 

8 * FREO DEPENDENT 

LOADS 

36 

SUBCASE 9 


• i 

37 

LABEL = 

SEGMENT 9 


38 

OLOAD = 

9 a FREO DEPENDENT 

LOADS 

39 

SUBCASE 10 


. .! 

40 

LABEL = 

SEGMENT 10 

■, i 

41 

OLOAD = 

10 i FREO DEPENDENT 

LOADS 

42 

SUBCASE 11 


‘ 1 

43 

LABEL = 

SEGMENT 11 

• . i 

• i 

44 

OLOAD = 

ll $ FREQ DEPENDENT 

LOADS 

45 

SUBCASE 12 



46 

LABEL = 

SEGMENT 12 


47 

OLOAD = 

12 S FREQ DEPENDENT 

LOADS 

48 

0UTPUT(XYPL0T3 



49 

PLOTTER NASTPLT, )50D£L 0,0 


50 

XPAPER.= 8,0 
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Or FOOn 


pf,2Z W 

Qij;a.rfV 



T 


■I? 

i 


i 

e3 


CARO 

COUNT 

51 

52 

53 

54 

55 

56 

57 
50 

59 

60 
61 
62 

63 

64 

65 

66 
67 
60 

69 

70 

71 


CASE CONTROL DECK ECHO 


YPAPER = 10.5 
KAXIS -- YES 
YAXIS = YES 
X6RI0 LIMES = YES 
YGRIO LINES « YES 
CUP.VELINE SYMBOL « L 
YLOG = YES 

XTITLE = FREQUENCY tHERTZ ) 

YTIYLE - GRID POINT DISPLACEMENTS { HAGNITUOE* INCH I 

TCUHVE = 14IT3RN) el8(T3RH» 

XYPLOT.XYPRI NT 01 SP RESPONSE I / 14 I T3RH I » 18 ( T 3RM ) 

TCURVE = 2IT3RMJ 

XYPLOTtXYPRINT OiSP RESPONSE 8 /2IT3RM) 

YTITLE = ELEMENT STRESSES I H&GNlTUDEiPSI 1 

TCURVE = 11 (35 ,11151 .11 (7J .liClOJ ,111 12J9 lit 14) 

XYPLOT.XYPRINT STRESS RESPONSE 1 Z 11 ( 3 ) . li( . lU 7) . 

imOt.lK 12). IK 14i 

TCURVE = U3» ,U5) ,117) ,U10) .1(12),1( 14) 

XYPLOT.XYPRINT STRESS RESPONSE 10 /1(3 ) , II 5) , 1( 71. 

1(10). 1112), 1(14) 

BEGIN BULK 


I'JFORHATi'JN MESSAGE 207, 


BULK DATA NOT SORTED, XSORT VHLL RE-ORDER DECK. 


I 

i 

I 

I 


'j 

ul. 
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OR POOR QUALiT/ 



. L 

a « 2 

« « 3 

• « 4- 

« ft 


C0RD2C 

L 


«3 

,0 

I 

&CUR12 

0.0 

1.0 

0.0 



C0UAD2 

4 

2 

2 

3 


C0UA02 

5 

2 

6 

7 

1 

C0UA02 

6 

2 

3 

4 


CQUA02 

7 

2 

7 

3 


CQUAD2 

8 

2 

4 

5 

T 

CQUA02 

10 

2 

a 

15 


CUUAi}2 

11 

3 

9 

16 


C QUA 02 

12 

3 

16 

IT 


CTRIA2 

1 

1 

1 

6 

•I 

CTRIA2 

2 

1 

1 

2 


CTRIA2 

3 

I 

10 

6 


CTRiA2 

9 

1 

8 

9 

T 

CYJOIN 

1 


1 

2 

i 

CYJOIN 

2 


10 

11 


OAREA 

L 

a 

3 

-1.0 

mm 

'J 

OAREA 

1 

16 

3 

1.0 

X 

□AREA 

1 

18 

3 

l.Q 


OAREA 

2 

8 

3 

-0.5 

n* 

OAREA 

2 

16 

3 

.5 


OAREA 

2 

18 

3 

.5 


OAREA 

3 

8 

3 

.5 


OAREA 

3 

16 

3 

-0.5 

mm 

( 

OAREA 

3 

18 

3 

-0.5 

4 

<• 

OAREA 

4 

8 

3 

1.0 


□AREA 

4 

16 

3 

-l.O 

mm 

.1 

OAREA 

4 

18 

3 

-1.0 

1 

OAREA 

5 

8 

3 

.5 

Ik* 

OAREA 

5 

16 

3 

-0.5 


OAREA 

5 

18 

3 

-0.5 

X 

'V 

OAREA 

6 

8 

3 

-0.5 

lift 

OAREA 

6 

16 

3 

.5 


OAREA 

6 

13 

3 

.5 

mm 

OAREA 

7 

3 

3 

-1.0 

1 

OAREA 

7 

16 

3 

l.O 


OAREA 

7 

18 

3 

1.0 


OAREA 

a 

3 

3 

-0.5 

1 

OAREA 

a 

16 

3 

.5 


OAREA 

8 

13 

3 

.5 


OAREA 

9 

8 

3 

.5 

T 

OAREA 

9 

16 

3 

-0.5 

I. 

OAREA 

9 

18 

3 

-0.5 


OAREA 

10 

8 

3 

l.O 


OAREA 

10 

16 

3 

-1.0 

1 ■ 

OAREA 

10 

18 

3 

-l.O 


OAREA 

11 

8 

3 

.5 


OAREA 

1 1 

16 

3 

-0.5 


□AREA 

l-l 

18 

3 

-0.5 

£3: 

OAREA 

12 

8 

3 

-0.5 


S 0 R 7 G 0 Q U L j ; D & 1 & E C Ihl 0 


8 

«0 



'.'■OF, POOR QUAUTY 


tiji 


1 



S 0 

R TED 

B U h 

R 0 

y 

. L 

.. 2 

• o 3 « 

. A 


. . 6 

IT 

I 

□ AREA 

12 

16 

3 

o 5 


□AREA 

12 

13 

3 

•of? 



FREQ 

1 

ITOOeO 

1750.0 

1777 

l795o 


CFRl 

1380.0 

1920.0 




f 

GROSET 


1 




1 

GRID 

1 


2.0 

30.0 

.0 


GRtO 

2 


3.1 

30.0 

.0 

If 

GRID 

3 


A.3 

30.0 

.0 

1 

GRID 

4 


5,2 

30.0 

.0 


GRID 

5 


7.1 

3C.0 

.0 


GRID 

6 


3.1 

A5.0 

.0 

r 

GRID 

7 


A.3 

A5.0 

.0 


GRID 

a 


5.2 

A5.0 

.0 


GRID 

9 


7.1 

AO 0 0 

.0 

r 

GRID 

10 


2.0 

60.0 

.0 

I. 

GRID 

ll 


3.1 

60.0 

.0 


GRID 

12 


A.3 

6G.0 

.0 

r 

GRID 

13 


5.2 

60.0 

.0 

1 

GRID 

lA 


7.1 

60.0 

.0 

»• 

GRID 

15 


7.1 

50,0 

.0 


GRID 

16 


8.^ 

AO.O 

-.25 

n 

GRID 

17 


9,7 

AO.O 

-.50 

\, 

GRID 

18 


0.5 

50.0 

.25 


GRID 

19 


9.7 

50,0 

.50 

f- 

HAFl 

1 

30.086 


.3 

7. A— A 

I ■ 

PARAM 

CYCIO 

81 





PARAM 

G 

.02 




r* 

PARAH 

GKAD 

FREORESP 



r 

PARAM 

KHAX 

2 




I. 

PARAM 

KMIN 

2 





PARAH 

LGKAD 

1 





PARAH 

NSEGS 

12 





PARAH 

RPS 

.0 





PQUA02 

2 

1 

.25 



“ 

PJUAD2 

3 

1 

.125 




PTRIA2 

1 

1 

.25 



1 » 

RLOADl 

1 

1 



100 


RLOAOl 

2 

2 



100 


RLOADl 

3 

3 



100 

ft 

RLOADl 

A 

A 



100 


RLOADl 

5 

5 



100 

r 

RLOAOl 

6 

6 



100 

!, 

RLOADl 

7 

7 



100 


RLOADl 

8 

8 



100 

p 

RLOAOl 

9 

9 



100 

[ 

RLOADl 

10 

10 



100 


RLOADl 

11 

11 



100 


RLOADl 

12 

12 



100 

r 

SPCl 

30 

6 

1 . 

THRU 

19 

L 

SPCl 

30 

123A56 

1 

iO 



0 A T a ECHO 


10 


i823«a5-ii032«O laSOil SFfti 
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i-i'.-JC iv1 1 

OF POOR QUALITY i 



« Loo 2 

TABLED! 100 
STBOi 0.0 
ENDOATA 


SORTED S U L K 0 A T A £ C H O' 


•• 3 o, o . . 3 o, £, oo 

1«0 lOOOoO l.O EL'DT 



S 


O O 


9 *o iO o 
OTSOl 









A. Description 

example uses the forced vibration capability with cyclic symmetry. 
Th'ii;puser input/output data pertain to harmonic representation. Frequency- 
dependent excitation is provided by both directly applied and base acceleration - 
loads. 

3. 1 nput 

1. Parameters: 

In addition to general input parameters, 

CYCIO = -1 harmonic cyclic input/output data 
KfllN = 0 minimum circumferential harmonic index 

KMAX = 2 maximum circumferential harmonic index 

NSE6S = 12 number of rotational ly cyclic sectors 
RPS = 600.0 revolutions per second 

BXTID, BYTID, BZTID \ Refer to TABLEDi bulk data cards to specify 
BXPTID, BYPTID, BZPTID/ magnitude and phase of base acceleration 

components. 

GKAD = FREQRESP\ Specify the form in which damping parameters are 
LGKAD = +1 / used. ! 

2. Constraints: 


Same as general input constraints. 

3. Loads: 

a) = A(f) specified on RLOADi bulk data cards. 

b) Base acceleration as shown in Figure 11. 

C. Results 

Results are shown ii Figures 12 through 20. 


Figures 12 and 13 present k = 0 results (subcase 1). The excit;ation consists 
of axial base acceleration and directly applied loads. The selected frequency 
band of excitation, 1700--1920 Hz, lies between the second out-of-plane d’se bending 
mode frequency (1577 Hz, k = 0, Table 2) and the first in-plane shear mode 
frequtr.cy (1994 Hz, k = 0, Table 2). Since the excitation is parallel to the 
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Esa 55?S3 ES3 ESJ EE3 EIS2 ES3 E23 ES3 CS3 EI3 EZ23 ■ ES3 


Figures 14 through 18 prasorit k ~ i results (subcases 2 (k = Ic) and 3 
(k = Is)). The excitation is due to lateral base acceleration only. Although 
the frequency band of input base acceleration is 1700-1920 Hz, the rotation of 
the bladed disc at 600 Hz (parameter RPS) splits the input bandwidth into tv/o 
effective bandv/idths: 

(1700 - 600) = 1100 to (1920 - 600) = 1320 Hz. and 

(1700 + 600) = 2^^ to (1920 + 600) = 2520 Hz. 

The only k = 1 mode in these effective bandwidths is the first torsional 
mode of the blade with the disc practically stationary (2460 Hz, k = 1, Table 2). 

This is shown by the out-of-plane displacement magnitudes of grid points 18 
(blade) and 3 (disc) respectively (Figures 14 (k = Ic) and 17 (k » Is)). The 
corresponding phase responses of these grid points are shown in Figure 16, 

Figures 19and 20 present k = 2 results (subcase 4 (k = 2c)). The excitation 

consists of directly applied k = 2c loads. The out-of-plane displacement magni- 
tude of grid point 18 (Figure 19) compares well with that obtained in example 2 
(Figure 7). Table 3 lists the out-of-plane displacement response of grid point 
18 as obtained in examples 2 and 3. The marginal difference in response in example 3 
is due to the Coriolis and centripetal acceleration effects at a rotational spped 
of 600 revolutions per second. 

No k = 2s loads are applied in this example (subcase 5). 
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D. Driver Decks and Bulk Dat a 

N A S T R A N EXECUTIVE C 0 N T .T 0 L DECK ECHO 


ID riASA,EXAllPLE3 

APP DISP 

SOL 8 

$ 

$ ALTER PACKAGE AS IN EXWiPLEZ 

$ 

TIME 12 $ IBM 370/3031 

DIAG 14,21 

CEND 



1 

i 

i 

1 

f 

i 

I 
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■•OffiGINAi. PAGE fS 
or f'COP. QUAUTV. 


CARO 

COUNT 

1 

2 


•C-a-.S £ 




0 

TITLE = f 
SUBTITLE 


iCEo mnnimcii &mvs£s- aafAiaos c\£C 5 .ss 

BLAOkO disc § eCT-C- EiSDSiS^GSOSeiiSS ££2B /LOS,D«HAGH 


t/m 



4 ■ 

fi- 


• ■ . - 


.... 


i 

r 

5 


SPC « 30 





i 

i 

6 


FREQ » i 
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This example uses the forced vibi'dtion capability vath cyclic syinnetry. 

The user input/output pertains to physical representation. Periodic loads are 
specified as functions of time on the secjmehts of the bladed disc corresponding 
to k = 2. For clarity of illustraticv! only* sinusoidal loads of varying 
amplitudes at a frequency of 1814 Hz are specified. The Fourier decomposition 
of these sine functions obviously contains contributions from first harmonic 
alone (£ = 1)-- the parameter LMAX accordingly has been set at 1 (£. = 0, Ic, 2l) 

B. Input 

1. Parameters: 

In addition to general input parameters, 

CYCIO = +1 physical cyclic input/output data 

KMIN = 2 minimum circumferential harmonic index 

KMAX = 2 maximum circumferential harmonic index 

LMAX = 1 maximum harmonic in the Fourier decomposition of periodic, 

time-dependent loads, 

NSEGS = 12 number of rotational ly cyclic sectors 
RPS = 600.0 revolutions per second 

GKAD = FREQRESPI Specify the form in which the damping parameters are 
LGKAD = +1 } used. 

2. Constraints: 

Same as general input constraints. 

3. Loads: 

,p"(t) = A(t) cos (ipr-(l>,^. , 

where n is the segment number, 

^represents k = 2, 

(f|) represents the total number of segments in the bladed 
disc, 

A(t) = A*sin (2Tr*1814*t) . 

P is specified on TLOADi bulk data cards. 



C. Results 

Results are presented -in Table 4 and ure In good agrae-Tient with those from 
example 3. 
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A. Description 

This example uses the forced vibration capability with, cyclic symmetry. 
The user input/output pertains to harmonic representation. Periodic loads are 
specified as functions of tim:; for the circumferential harmonic index k « 2. 

For clarity of illustration only, sinusoidal loads are selected. 

3. Input 

1. Parameters: 

In addition to general input parameters, 

CYCIO = -1 harmonic cyclic input/output data 
KHIN = 2 minimum circumferential harmonic index 

KHAX = 2 maximum circumferential harmonic index 

LItAX = 1 maximum harmonic in the Fourier decom; osition of periodic, 

time-dependent Toads. 

NSEGS =12 number of rotational ly cyclic sectors 
RPS = 600.0 revolutions per second 

GKAD = FREQRESP\ Specify the form in which the damping parameters 
LGKAD = +1 / are used. 

2. Constraints: 

Same as general input constraints. 

3. Loads: ■ 

P^^(t) =. A-fin (2T;*1314*t) , 

specified on TLOADi bulk data cards. 

C. Results 

Results are presented in Table 4 and agree well with those from example 3. 
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D. Driver Decks and Bulk Data 

K A S T R A N EXECUTIVE CONTROL DECK ECHO 
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Example 

rio. 

Finite 

Element 

Model 

of 

Applied loads specified as functions of 

Base 

Acceleration 

Rotational 

Speed 

Frequency (sinusoidal) 

; Time (periodic) 

Physical 
Componer ts 

Circuin. Harmonic 
Components 

Physical 

Components 

Circum, Harmonic 
Components 

1 

Complete 

Structure 

Yes 




No 

No 

2 

Cycl ic 
Sector 

Yes 




No 

. 

No 

’ 

1 
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Cyc 1 i c 
Sector 


Yes 

' 



Yes 

' 

Yes I 
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Cycl ic 
Sector 
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Yes 


No 

Ves . 1 

. 1 

1 
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Cycl ic 
Sector 


! i 


Yes 

No 

Ye': 1 
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*■ k is the circumferential harmonic index 


** Mode No. 4 for k = 0 at 1994 Hz represents an in-plane shear mode not excited 
by the applied forces. 






































TABLE 4: COF^PARISCM OF RESPONSE AT 1814 Hz 


Grid Pt.Disp. or 
Elem. Stresses 

Example 3 

Example 4 

Example 5 

k = 2c (subcase 4 ) 
Maq. (in)/Phase(deq) 

Segment 1 (subcase 1 ) 
iiaq.(in)/Phase(deq) 

k = 2c (subcase 4) 
I>1aQ.(in)/Phase(deq) 

8 (T3RH), 

■ . 

- 

5.4297 E-4/82.6 

5.4299 E-4/82.6 

5.4299 E-4/82.6 

18 (T3RiM), 

4.9177 E-4/266.8 

4.9180 E-4/266.8 

4.9180 E-4/266.8 


1.4841 E 3/84.7 

1.4842 E 3/84.7 

1.4842 E3/84.7 

” “yy.l 

2.0891 E 2/83.4. 

■■ ■ ■ 

2.0892 E 2/83.4 

2.0892 E2/83.4 

” <')• Vi 

1.0774 E 2/64.7 

1.0775 E 2/64.7 

1.0775 E2/64.7 

" ('<»• “xx,?* 

1.4677 E 3/263.3 

1.4678 E3/263.3 

1.4678 E3/263.3 

11 (12). 

2.2489 E 2/260.3 

2.2491 E 2/260.4 

2.2491 E2/260.4 

11 (14). 

1.85T0 E 2/253.0 

1.8511 E 2/253.0 

1.8512 E2/253.0 


* Fibre distances 1 and 2. 
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Segment 8 


Figure 1; HASTRAN Model of the 12-Bladed Disc 
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Figure 6 



1.68 . 1.72 1.7G 1.80 1.84 

FfiEOUENCT (HEHTZI 

m IT3RMI , 18 (T3RMI 

FORCED ViSRRTION RHOLTSIS OF RO T R II NG C TCL i C STRUCTURES 
BEADED DISC EXAMPLE 2 ICTC MCDEL.FREO LORDS , PHT S I C RE I/O) 
SEGMENT 1 SUSCnSE 1 


1.88 E3 


Figure 7 
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Figure 14 
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